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Abstract 

For nearly a century, more mass has been measured in galaxies than is contained in the luminous 
stars and gas. Through continual advances in observations and theory, it has become clear that the 
dark matter in galaxies is not comprised of known astronomical objects or baryonic matter, and 
that identification of it is certain to reveal a profound connection between astrophysics, cosmology, 
and fundamental physics. The best explanation for dark matter is that it is in the form of a yet 
undiscovered particle of nature, with experiments now gaining sensitivity to the most well-motivated 
particle dark matter candidates. In this article, I review measurements of dark matter in the Milky 
Way and its satellite galaxies and the status of Galactic searches for particle dark matter using 
a combination of terrestrial and space-based astroparticle detectors, and large scale astronomical 
surveys. I review the limits on the dark matter annihilation and scattering cross sections that can 
be extracted from both astroparticle experiments and astronomical observations, and explore the 
theoretical implications of these limits. I discuss methods to measure the properties of particle 
dark matter using future experiments, and conclude by highlighting the exciting potential for dark 
matter searches during the next decade, and beyond. 
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1. Dark Matter Searches: Historical Context 



1.1. Establishment of dark matter 

When it first became possible to measure the motions of stars in the Milky Way, its structure 
and dynamics was revealed to a level not previously possible. On examination of the brightest 
stars in the disk of the Milky Way, Oort found that the stars were only able to account for about 
one-third of the total mass contained within the disk (Oort, 1932). Oort deduced that the disk 
contained two-thirds "dark matter" that was predominantly in the form of stars less luminous than 
the Sun, as well as gas and dust in the interstellar medium. Nearly thirty years later, Kinman 



(1959) analyzed the velocities of Milky Way globular clusters, and deduced possible deviations 
from a pure bulge and disk mass model of the Galaxy, known at the time as "Schmidt's Law." 
Though there were large error bars associated with his measurements, the data was consistent with 
a linearly rising mass distribution beyond the disk. Kinman suggested that the data be held in 
reserve due to uncertainties in the globular cluster distances and the local circular velocity. 



Babcock (1939) was the first to measure the mass distribution of M31, the nearest galaxy with 



luminosity comparable to the Milky Way. From the measurement of radial velocities in the optical 
emission line region extending out to approximately 20 kpc from the center of M31, Babcock 
determined that the total mass-to- light ratio increases in the outer regions. He suggested that 



absorption of light was the cause of this increasing mass-to-light ratio. Two decades later, Kahn 



and Woltjer ( 1959 ) were the first to connect the dynamics of local galaxies to the universe on larger 



scales. These authors noted that if the Milky Way and M31 were once receeding due to Hubble 
expansion, and are now approaching at the observed velocity, there must be ten times more unseen 
material within both of the galaxies and within the inter-galactic medium to reverse their initial 
separation. They made the suggestion that this material was in the form of intergalactic gas and 
dust. 

The results on the mass distribution of M31 were further corroborated by the 21-cm observations 
of van de Hulst et al. (1957) and Roberts and Whitehurst (1975), which clearly indicated that the 
rotation curve of this galaxy is flat out to approximately 35 kpc from the center, and did not fall-off 
in the Keplerian fashion as expected if the stars and gas dominated the gravitational dynamics. 
Flat rotation curves were identified in larger samples of nearby spiral galaxies through a variety 
of techniques at several observatories ( Bosnia, 1978 Rubin et al. 1980). After several years of 



dealing with systematics, such as the effect of beam smearing, it became the widespread consensus 
amongst astronomers that the results were not due to observational limitations, and that rotation 
curves of many spiral galaxies remained flat in the outermost regions (Sofue and Rubin, 2001). 

The seemingly linear rise in the mass distribution of the Milky Way beyond the disk, combined 
with the measured rotation curves of M31 and nearby spiral galaxies, provided a fresh perspective 
on well-known observations obtained from larger and more distant systems. About a year after 
Oort published his results on the kinematics of the stars in the Milky Way, Fritz Zwicky famously 
deduced that there was hundreds of times more mass in the Coma cluster of galaxies than was 
contained within the luminous galaxies (Zwicky, 1933). Zwicky considered a few explanations for 
his data: the constituent galaxies contained much more mass than he ascribed to them, there was 
undetected intergalactic matter, or perhaps even the known physical laws were not appropriate for 
the mega-parsec scales he was considering. Upon noting that the cluster may be dominated by 
"dark matter," Zwicky both used Oort's phrase for the first time in the context of extragalactic 
astronomy, and also in a manner that most closely resembles its meaning in the modern scientific 
vernacular. 

As incontrovertible evidence built indicating that galaxies and clusters of galaxies require much 
more mass that is contained in the luminous material, the question became eminently focused: 
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What is the nature of the dark matter that dominates the mass of astronomical systems? While 
there has been, and continues to be, a variety of theoretical of explanations to this question, two 
particularly compelling explanations have lent themselves to rigorous investigation. On the one 
hand there is the plausible conjecture, entertained by many of the aforementioned astronomers, 
that our inventory of astronomical bodies is incomplete, in particular we do not have a good census 
of bodies that are faint sources of radiation that escape detection. On the other hand, a theoretical 
picture has been developed that predicts that dark matter is a new elementary particle that is 
distinct from the known particles of nature. In order to set the stage for the results that are 
presented in this article, it is informative to review the historical establishment of these two classes 
of dark matter candidates, the predictions that they make, and the searches for signatures of them. 

1.2. MACHOs: Review of searches and results 

We do not have a complete understanding of the population of low mass stars, stellar remnants, 
and planetary mass bodies in the Galaxy. Compact astronomical bodies that constitute a significant 
component of the mass of the Galaxy are historically referred to as Massive Astrophysical Compact 
Halo Objects (MACHOs). Because of their intrinsic faintness, or because they in some scenarios 
emit no radiation at all, MACHOs typically cannot be identified from direct imaging. This leaves 
gravitational interactions, and in particular gravitational lensing, as the best available technique 
to search for MACHOs. 

If a massive body, in this case a MACHO, intervenes between the Earth and a background 
source star, the massive lensing body creates multiple images of the background star. Though the 
images are too small to be resolved by modern telescopes, the source star brightens on a timescale 
that is related to the mass of the intervening object and its velocity. The amplification of the 
source is the ratio of the area of the images that are created to the area of the unlensed source 
star. For a stellar mass object moving at a typical velocity through our Galaxy, the characteristic 
brightening timescale of the star is in the range of a few weeks to a year. The above phenomenol- 



ogy of microlensing, initially predicted nearly a century ago (Einstein, 1936), is now a standard 



phenomenological tool in modern astrophysics (Mao 2012). If the dark matter in our Galaxy is 



comprised of faint stars, stellar remnants, planetary mass objects, black holes, or some combination 



of compact objects thereof, Paczynski (1986) suggested that a long term monitoring program of 



stars in the nearby Large Magellanic Cloud (LMC) would reveal dark matter in this form, even 
though the objects could not be identified directly. 

During the past two decades, the MACHO project, the Optical Gravitational Lensing Exper- 
iment (OGLE), and the Experience pour la Recherche DObjets Sombres survey (EROS) project 
have searched for MACHOs in the Milky Way via microlensing. Though microlensing events were 
indeed detected by these experiments, it is generally concluded that they do not have significant 
enough abundance to dominate the mass of the Galactic halo. As no conclusive signature for 
MACHOs was found, these experiments combined to place strong upper limits on the number of 
MACHOs in the Galaxy in the mass regime of approximately 1(T 7 - 30 M (|Alcock et ah] |1998 



Tisserand et al., 2007 Wyrzykowski et al. 2011). 



For MACHOs greater than about a hundred solar masses, the timescale of the typical microlens- 
ing event for a lensing object in the halo is a few years, which is near the lifetime of the surveys 
themselves. Therefore the efficiency for detecting objects in this mass range drops significantly, 
and other methods must be utilized to identify them. Two such methods rely on studies of the 



population of wide binaries in the Galaxy ( Bahcall et al. 1985 ) , and on the kinematics of the Galac- 



tic disk (Lacey and Ostriker, 1985). The former method utilizes the fact that wide binaries are 



broken up by encounters with MACHOs, either by weak and diffuse encounters that occur during 
the entire lifetime of the binary system, or by catastrophic encounters that break up the binary 



during a single interaction with a MACHO. The disruption timescale for binaries with a separation 
> 0.1 pc is less than approximately 10 Gyr for MACHOs with mass > 10 Mq, so the number of 
wide binaries with this separation or larger provides a constraint on the number of MACHOs in 
this mass regime. The most conservative constraints find that MACHOs of mass ~ 10 — 10 7 Mq 
comprise no more than ~ 50% of the Galactic halo (Yoo et al. 2004 Quinn et al. 2009). Above 



this mass scale, the indirect arguments from the velocity dispersion of the Galactic disk may be 



interpreted as placing strong constraints on MACHOs with mass > 10 7 Mq (Lacey and Ostriker 



1985) 



Though microlensing studies find that MACHOs do not dominate the mass of the Galaxy, there 
have been hints that significant populations of stellar remnants and low mass bodies do indeed exist 



in the Galactic halo. Oppenheimer et al. (2001) find a substantial population of old, cool white 
dwarfs with high proper motions, and that these objects are abundant enough to comprise > 2% 
of the Galactic halo. Subsequent analyses has however indicated that these objects are more likely 
members of the thick disk and the stellar halo, and not components of the more extended "dark" 



halo of the Galaxy (Reyle et al., 2001 Flynn et al., 2003) 



Though the analysis above is phrased in the context of constraints on low mass stars and stellar 
remnants, microlensing searches are sensitive to exotic objects such as primordial black holes if they 
constitute a component of the dark matter. Rather than forming as an end state of stellar evolution, 
or via growth from smaller seed black holes in the centers of galaxies, primordial black holes may 
have a wide range of masses. Primordial black holes may have formed in regions of the early 
universe that were sufficiently overdense that they halted expansion and recollapsed. They may 
also have formed during a phase transition in the early universe. Primordial black holes less massive 
than approximately 10 15 g ~ 5 x 10~ 20 Mq would have evaporated via Hawking radiation. Though 
reasonable physical processes predict their formation, there is no strong theoretical motivation 
for primordial black holes to be as abundant as the observed dark matter. General astrophysical 
limits on the abundance of primordial black holes of a given mass result from searching for their 



gravitational interactions or for particles that may be produced during their evaporation (Josan 



et al. 2009). During the next few years, microlensing observations are expected to further improve 
the limits on primordial black holes with mass < 10~ 7 Mq (Griest et al. 



2011). 



The above discussion indicates that Galactic searches for dark matter in the form of MACHOs 
have reached a significant level of maturity. These results have firmly established that dark matter 
does not dissipate energy to clump into objects as massive as stars or planets, indicating that the 
dominant matter component in Galactic halos is more diffusively distributed. This is probably the 
most direct empirical evidence that we currently have for the existence of particle dark matter. 

1.3. WIMPs: Theoretical motivation and initial search results 

The lack of MACHO detection was not surprising considering theoretical developments in cos- 
mology and particle physics. Going back over thirty years observations of the abundances of light 
elements have indicated a baryon (neutron and proton) density in the universe considerably less 
than the critical density, implying that some form of non-baryonic matter must be cosmologically- 
significant (Gott et al. 1974). Even at the time of these earliest calculations, this was consistent 



with the determination that the extended mass distributions around galaxies comprised a signif- 
icant fraction of the critical mass density of the universe (Ostriker et al. 1974). Non-baryonic 



matter was also found to be very appealing from the perspective of large-scale structure formation, 
because density fluctuations grow much earlier as compared to a universe dominated by baryonic 



matter (Peebles, 1982). This nicely accounts for fluctuations in the Cosmic Microwave Background 



(CMB) at the level of one part in 10 5 . These lines of evidence indicate that the dark component of 
spiral galaxies that was initially deduced through rotation curves may be connected to observations 
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on larger cosmological scales, and a profound connection was developing between the production 
of dark matter in the early universe, large scale structure formation, and the mass that dominates 
galaxies and clusters of galaxies. 

As these cosmological calculations of the matter abundance in the universe improved and be- 
came more precise, an interesting tension came to the forefront between these calculations and the 
developing standard model of particle physics. In particular it was not clear that the standard 
model of particle physics contained a particle with the necessary properties to be significant on 
cosmological scales. The most natural particle to consider for the non-baryonic dark matter was 
the neutrino. A neutrino with mass of a few electron volts (eV) was appealing because if it was 
in equilibrium in the early universe, the present mass density of neutrinos would be near the crit- 

However, under more detailed scrutiny neutrinos at 

On the one hand, 



ical density (iCowsik and McClellandl |1972h. 



this mass scale turn out to be disfavored as a candidate for the dark matter, 
phase space limits strongly constrain neutrino dark matter with mass less than a few hundred 
eV (Tremaine and Gunn, 1979). Further, numerical simulations of the large scale structure dis- 
tribution in the universe revealed that neutrinos with mass below tens of electron volts were too 
"hot" to explain the observed galaxy distribution; because they were relativistic when they decou- 



pled they produced significantly fewer low mass galaxies in comparison to the observations (White 



et al. 1983). 



Though neutrinos with mass of a few hundred eV or less became disfavored as the dominant 
component of non-baryonic dark matter, this analysis within a cosmological context motivated a 
general theoretical framework for determining the abundance of a stable particle species that was 
in thermal equilibrium in the early universe. Neutrinos with mass greater than about 1 GeV that 



fall out of equilibrium while non-relativistic were found to be cosmologically-significant (Lee and 



Weinberg, 1977), and that generically heavy leptons with similar interactions could be the dominant 
component of mass in galaxies and clusters of galaxies (Gunn et al. 1978). These results pointed 



to a preferred scale, the weak scale, to describe the interactions of a cosmologically-significant 
component of non-baryonic dark matter. For a standard thermal history in the early universe, the 
abundance of a particle is related to its thermally-averaged annihilation cross section times relative 
velocity as 

3 x KT^cmV 1 

<<w> * — — (1) 



flZeldovich[ |1965[ |Chiu[ |1966| |Steigman[ |1979[ |Scherrer and Turner} |1986[ |Steigman et al.[ |2012[ ). 
Recent cosmological measurements find that the fraction of non-baryonic dark matter relative to the 
critical density is f^DM^ 2 — 0.11 (Komatsu et al. , 2011 ). For this value of S^dm^- 2 ) for particles in the 
approximate GeV mass range, this scale for the annihilation cross section is characteristic of weak 
interactions. A particle with interactions at this scale is broadly referred to as a weakly-interacting 
massive particle (WIMP). 

Because of their weak scale interaction strength, WIMPs are detectable via non-gravitational 
methods. In addition to the prediction of the annihilation strength of WIMPs, Equation [T] leads to 
a basic prediction for the WIMP scattering cross section on ordinary matter (i.e. quarks). Com- 
bining this scattering cross section with the estimation for the local number density of WIMPs, the 
interaction rate is large enough for them to have been detected over two decades ago in the first 
generation of low temperature germanium experiments that were primarily designed to search for 
signatures of neutrinoless double beta decay. While it was certainly possible that these experiments 
would detect WIMPs and render the non-gravitational detection of particle dark matter relatively 
straightforward, the first two experiments to systematically search for WIMPs reported null re- 



sults (Ahlen et al. 1987, Caldwell et al. 1988). The reported constraints ruled out the simplest 
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model for WIMP interactions with ordinary matter, and in the process first showed that if WIMPs 
comprise the dominant component of dark matter in galaxies, physics must be invoked to suppress 
the scattering cross section relative to that derived from the most basic cosmological arguments in 
Equation [TJ 

1-4- Modern perspective on dark matter searches 

The results from MACHO searches and the seminal WIMP searches have formulated our modern 
picture of particle dark matter. While MACHO searches have largely ended, WIMP searches are 
continuing to pick up incredible momentum, largely because there are still very compelling regions 
of theoretical parameter space left to explore. The size and scope of searches for particle dark 
matter have expanded considerably since the first set of experiments that were designed to search 
for them. For these reasons, it is important to have a better classification of the types of WIMP 
searches going on around the world at present, and what each of them will ultimately be able to 
tell us about the properties of dark matter. The focus of this article is on two broad classes of 
observations that are used to probe the particle nature of dark matter: searches with astroparticle 
experiments and colliders, and searches using astronomical data sets. 

1.4-1- Astroparticle and collider searches 

There are three types of dedicated astroparticle and collider experiments actively searching for 
particle dark matter. Each type of experiment is sensitive to a different aspect of dark matter 
interaction with ordinary matter. 

• Indirect Searches: Indirect dark matter searches measure the annihilation and / or decay prod- 
ucts of dark matter from astrophysical systems. Schematically, they measure the rate for DM DM 
— > SM SM, where DM represents the dark matter particle and SM represents any standard model 
particle. In many instances, the particle represented by SM is unstable, and decays into particles 
(for example, photons or neutrinos) that are observable in detectors. Comparing to Equation [TJ 
the cross section probed by indirect searches is most closely related to the process that sets the 
abundance of the dark matter in the early universe, assuming that the dark matter was once in 
thermal equilibrium. In order to best interpret the results from indirect searches, we must have 
a good idea as to both how the dark matter is distributed in halos, and what standard model 
particles the dark matter preferentially annihilates or decays into. 

• Direct Searches: Direct dark matter searches, an example of which was described above, mea- 
sure the scattering of dark matter off of nuclei in low background underground detectors, DM SM 
—J- DM SM. In these experiments, the collision of dark matter is deduced through energy input 
into particles in the detector. At the most fundamental level, the SM particle is represented by up 
and down quarks that make up the protons and neutrons within nuclei, so direct detection experi- 
ments ultimate measure the strength of the interactions between WIMPs and these types of quarks. 
Though there is no specific cosmological prediction for the WIMP-nucleus scattering cross section 
similar to the prediction that leads to Equation [TJ as discussed throughout this article modern 
experiments are steadily improving limits and are probing well- motivated theoretical models. 

• Collider Searches: Collider searches measure the production of dark matter through the 
collision of high energy particles on Earth, SM SM — > DM DM. Hadron colliders, such as the 
Tevatron at Fermilab and the Large Hadron Collider (LHC) at CERN, collide protons to produce 
new stable particles either directly or through processes that are initiated by quarks and gluons. 
New stable dark matter particles produced directly cannot be observed, so these events must be 
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tagged by observable initial state radiation. However, the cross section for this process is typically 
small. On the other hand, the cross section is larger for the production of quarks and gluons, which 
eventually decay down into the lightest stable particle in the spectrum. In this process, much of the 
energy of the collision is not transferred to final state particles, so it is often not straightforward 
to reconstruct the kinematics of the collision and determine the mass of the stable particle that 
is produced. Lepton colliders, such as the Large Electron Positron (LEP) collider that operated 
from 1989-2000 at CERN, transfer all of the initial state energy into final state particles. LEP has 
produced a lower bound on the existence of stable charged particles of approximately 100 GeV, and 
through the non-observations of perturbations to standard model cross sections motivates existence 



of new stable particles near the electroweak scale (see the arguments discussed in Feng (2010a)). 



If dark matter is in the form of a new elementary particle, understanding its interactions requires 
measuring the cross section with standard model particles via each of the three methods described 
above. From the three sets of measurements, the theoretical challenge will be to determine whether 
the WIMP is a single stable particle, or it is embedded into a larger theory of high energy physics, 
perhaps one that contains a spectrum of particles with masses larger than that of the WIMP. 

Though it is probably most prudent at this stage to treat the respective cross sections described 
above as independent, under some simplifying theoretical assumptions it is possible to connect the 
cross section as measured by one experiment with that measured in another. A straightforward 
example is if the WIMP annihilates pre-dominantly to standard model fermions, and there are no 
resonant interactions that serve to enhance the annihilation cross section. In this case a measure- 
ment of the WIMP-fermion coupling from the strength of the annihilation cross section translates 
directly to the elastic scattering rate off of nucleons in underground detectors. Typically in these 
models, the WIMP is the only particle the constitutes the dark matter, with all other particles in 
the theory at a significantly larger mass scale so they do not affect the interaction rate of WIMPs 
with standard model particles. Models along these lines are discussed in more detail in Sections [6] 
and [3 

Given the compelling cosmological motivation for WIMPs, a vast number of theories have been 
developed that embed a particle (or particles) with these properties into a consistent framework. 
The most famous and well-studied example of WIMP dark matter is the neutralino in super- 



symmetry (Goldberg, 1983). Dark matter candidates also arise from theories of universal extra 



dimensions 1 


Cheng et al. 


2002 Hooper and Profumo 


2007), as well 


as many others ( 


Duffy and van 


Bibber 


2009; Feng, 2010b) (see the recent reviews of 


Bertone et al. 


(2005 


) and 


Bertone (2010) for 



a more complete list of dark matter candidates). In many cases these theories can be appealing not 
only from the perspective of having a well-motivated dark matter candidate, but also because they 
are often linked to various other unsolved problems in high energy physics, such as the stabilization 
of the Higgs mass, the unification of the forces at the GUT scale, or the explanation of the matter- 
antimatter asymmetry. A stable dark matter particle is thus embedded into a larger theoretical 
edifice, that almost always contains a larger spectrum of particles at higher energy scales. It is these 
theories that we hope to study over the next several years with the three separate astroparticle and 
particle experiments that are outlined above. 



1.4-2. Astrophysical searches 

In modern astrophysics, two types of observations are particularly sensitive to the nature of 
particle dark matter: 

• Abundance of low mass (dwarf) galaxies: Very roughly, the minimum mass dark matter 
halo that can form in the universe is related to the velocity of the dark matter at the time it 
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decoupled from the rest of the standard model particles in the early universe. In deriving Equation[T] 
the standard assumption was made that the dark matter mass is larger than the temperature 
at decoupling. If the opposite is true, and the dark matter mass is less than the decoupling 
temperature, then the dark matter decouples not as cold dark matter (CDM), but as so-called 
warm dark matter (WDM) (with hot dark matter in the form of standard model neutrinos already 
being ruled out). As the ratio of the dark matter mass to the decoupling temperature decreases, 
so too does the number of low mass dark matter halos. Though this argument provides a general 
guiding theoretical principle, numerical simulations are required to most accurately predict the 
abundance of dark matter halos as a function of halo mass given properties of the dark matter 
particle. 

In order to determine the minimum mass dark matter halo in the universe, a mapping between 
dark matter halo mass and galaxy luminosity is required for the faintest known galaxies. However, 
the uncertain theoretical mapping between galaxy luminosity and dark matter halo mass makes it 
difficult to interpret modern observations as a direct constraint on particle dark matter. Further, 
even though modern surveys now identify galaxies several magnitudes fainter than the Milky Way, 
these surveys are notoriously incomplete. In spite of these difficulties, steady progress is being 
made on both the observational and the theoretical fronts, and as discussed in significant detail in 
Sections [4] and [5j interesting parameter space for dark matter particle models is well within reach. 

• Central density of dark matter halos: The properties of particle dark matter are also strongly 
linked to the mass distribution in the central regions of galaxies. The two most straightforward 
properties of the dark matter that these observations probe are 1) the phase space density of 
dark matter, and 2) the self-interaction, or scattering rate, DM DM — >• DM DM. Because of 
its low velocity at decoupling, the phase space density of the WIMP is large, many orders of 
magnitude larger than the phase space density of a standard model neutrino that decouples while 
it is relativistic. Very broadly, due to the low phase space density neutrinos favor constant density 
cores in the centers of dark matter halos, while perfectly cold WIMPs favor more steeply rising 
central densities. However, because of limited observational resolution, and because the effect 
of baryons on the central densities of dark matter in galaxies is not well-understood, mapping the 
measurements of the central densities of dark matter halos into a constraint on particle dark matter 
is not straightforward. 

The self-interaction property is relevant because in order to form extended halos that are stable 
over cosmological timescales, there should be minimal interactions between dark matter particles 
in halos through non-gravitational methods. For a power law density profile, significant amount 
of interaction leads to the mixing of the "cold" (low velocity dispersion) central component of the 
galaxy with the "hot" (high velocity dispersion) component. It is straightforward to understand 
why, in the context of WIMP dark matter, interaction between dark matter is insignificant. From 
the weak scale interactions, the characteristic scale for DM-DM scattering is typically ~ 10 -40 cm 2 . 
However, this cross section is so small that WIMP- WIMP scattering will effectively never happen 
in a dark matter halo. Assuming a density of ~ GeV/cm 3 , and a WIMP of mass 100 GeV, for a 
characteristic WIMP velocity of ~ 100 km/s, the self-interaction rate of WIMPs is at least thirteen 
orders of magnitude longer than the Hubble time. 

The best modern limits on the self-interaction rate of dark matter do not arise from the obser- 
vations of the central densities of galaxies, but rather result from the interpretation of the collisions 



of clusters of galaxies, the most famous being the so-called bullet cluster (Clowe et al. 2006). In 
this system, three separate components are observed: the galaxies within the clusters (which are 
assumed to be collisionless) , the gas distribution as traced in X-rays, and the total mass as de- 
termined from a gravitational lensing map of the matter distribution. The collisionless galaxies 
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are observed to follow the total mass distribution, while the gas components are observed to be 
more strongly interacting. From these observations, the upper limit on the self-interaction rate 
of the dark matter for a 100 GeV WIMP is approximately 10~ 24 cm 2 , which is many orders of 
magnitude larger than predicted by the theory of weak interactions. Thus any model that predicts 
significant interaction between dark matter should have a much larger cross section than that which 
characterizes WIMPs. 



Within the cosmological context of the standard Lambda Cold Dark Matter (LCDM) model, 
particle dark matter is approximately five times more abundant than baryons, about 22% of the 
total energy density in the universe, and its energy density is a little less than a third of the dark 
energy density (Komatsu et al. 2011). Though the cold dark matter paradigm is fundamental 
to modern cosmology and crucial to the interpretation of many aspects of modern astrophysics, 
a clear validation of it has remained elusive. Detection of particle dark matter with the above 
properties will almost certainly be the best means to verify the cold dark matter paradigm, providing 
further validation of the LCDM cosmological model. Approaching a half of a century after the 
initial particle dark matter candidates were first postulated, and they were connected to the mass 
distribution around galaxies and to cosmology, within the next several years we are finally in a 
position to obtain a better understanding of the properties of WIMP dark matter, if it does indeed 
exist. 
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2. Scope of this article 



There are a number of excellent recent reviews on theoretical aspects of particle dark mat- 



ter (Bertone et al. 2005 Hooper and Baltz, 2008; Feng 2010a, Bergstrom, 2012) and experimental 



searches using astroparticle data (Gaitskell, 2004 Porter et al. 2011). The former articles in 



particular highlight how a large focus of the theoretical research has revolved around developing 
extensions to the standard model that contain either a single dark matter particle candidate or a 
spectrum of them. 

Now that we have new experimental data sets in hand, and these experiments are beginning to 
encroach into the regime of cosmologically well-motivated dark matter candidates, it is important 
to understand what information can be extracted from experimental data. The goal of this article 
is to review both the theoretical and experimental developments in the field of direct and indirect 
searches for Galactic dark matter, and to understand the implications of these results. As motivated 
in the introduction, Galactic searches for dark matter not only involve interpretation of data from 
astroparticle experiments, but also data from more classical astronomical surveys that are sensitive 
to the microscopic properties of the dark matter. For this reason, I devote the beginning sections of 
this article to understanding what we have learned in recent years from astronomical surveys, and 
examine both their direct implications on particle dark matter, and also the interplay that they 
have with astroparticle searches that are reviewed in more detail in subsequent sections. 

This article is specifically organized as follows. In Section [3] I review the observational evidence 
for dark matter in the Milky Way, and what is known about its distribution from a variety of 
different astronomical observations. In Section |4j I examine the population of Milky Way satellites, 
which are exceedingly important for several aspects of dark matter searches. In Section [5j I connect 
these observations to the theoretical predictions of the dark matter distribution in Milky Way-mass 
galaxies from modern numerical simulations. 

The astroparticle analysis begins in Section [6j where I review theoretical developments in indi- 
rect dark matter detection. In Section [7j I review experimental results in indirect detection, and 
review the implications for theoretical models. In Section |8j I review theoretical developments in 
direct dark matter detection. In Section [9j I review experimental results in direct detection, and 
review the implications for theoretical models. 

In Section [TO] , I outline future directions for direct, indirect, and astrophysical searches for 
Galactic dark matter. In this section I discuss some very recent results that connect aspects 
of direct and indirect detection experiments to searches for dark matter at the LHC. I close by 
examining future directions for astrophysical dark matter searches, and focus on what about the 
properties of dark matter can be extracted from these observations. These will become particularly 
important if WIMPs are not discovered during the next decade. Example observations I focus 
on are kinematic measurements of the stars in the Galaxy from wide-field surveys, searches for 
Galactic satellites, and gravitational lensing measurements in our Galaxy and beyond. I hope to 
bring to the forefront how astronomical observations can test, constrain, and possibly uncover the 
properties of dark matter. 
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3. Dark matter in the Milky Way 



The Milky Way (MW) is the best-studied of all galaxies. The three primary components of 
the MW in visible light are the bar /bulge, the disk which contains the spiral arms, and the stellar 
halo. The dominant component of the disk by mass is the thin disk, which has had a sustained 
star formation for nearly ten billions years, and thus has a diverse population of stars with a wide 
range of ages and metalliticies. The scaleheight of the thin disk perpendicular to the Galactic 
plane is approximately 300 pc, and the total mass of the thin disk is approximately 5 x 10 10 M . 
The second component is the thick disk, which has a scaleheight of approximately 1 kpc and is 
composed of older stars that are more metal-poor than the thin disk. By mass, about 90% of the 
visible material is contained in stars, while the remaining 10% is in the form of gas and dust. The 
total absolute visual magnitude of MW is estimated to be My = —20.9 (van den Bergh, 2000), 
corresponding to a total luminosity of approximately a few times 10 10 Lq. 

In the era of large scale surveys, extraordinary progress that has been made in understanding 
MW properties (Juric et al. 2008). Yet because our position viewing it from the interior, many 
aspects of the MW remain a mystery. For example, relatively little is known about the properties 
of the spiral arms in comparison to external galaxies of a similar type. Further, we do not yet have 
a firm understanding of the nature of the bulge. There is a growing line of evidence that the MW 
does not contain a so-called classical bulge, but rather a "pseudo-bulge" that may evolve in concert 
with the disk. On more theoretical grounds, we do not yet have a good understanding of how rare 
the MW is in terms of its number of satellites, globular clusters, merger history, and metallicity 
distribution in comparison to other galaxies of similar mass, luminosity, and morphological type. 

This section is primarily focused on using the distribution of stars and gas in the MW to extract 
information on the mass distribution, both visible and invisible. This information is extracted 
from a variety of different stellar populations, which act as kinematic tracers of the underlying 
potential. Though as most broadly classified, the MW stars can be grouped into bulge, disk, and 
halo components, it is possible to further divide these populations in sub-populations. These sub- 
populations can be classified according to such properties as their distance, velocity dispersion, 
age, or metallicity. The salient point is that the stars in each population had a similar formation 
history, so all, or a subset of these properties, are common to all stars within it. 

Because we are interested in uncovering the nature of the dark matter, there are, broadly, 
three particularly important components to separately consider: the Galactic center, the Solar 
neighborhood, and the distant outer part of the dark matter halo. The status of our knowledge of 
the dark matter distribution in each of these three regimes is discussed, paying specific attention 
to what quantities are best measured by the data, and the systematics that are incurred in the 
measurements. The dark matter distribution in each of these regimes is determined by utilizing 
different populations of stars and gas; it is important to review each of these observations in detail 
in order to gain an appreciation for the systematics that are propagated to particle dark matter 
searches discussed in the sections below. 



3.1. Dark matter in the Galactic center 

The Galactic center is a notoriously difficult region to study from our position within the MW 
disk. It is obscured by approximately 30 magnitudes in visual wavelengths, so, aside from a few 
windows of below average obscuration, the best method to observe it are through infrared and radio 
wavelengths. More recently it has becomes possible to resolve the properties of individual stars to 
determine the shape, ages, metallicities, and kinematics of bulge stars. 

Stars that belong to the bulge dominate the luminosity within approximately 1 kpc of the 
Galactic center. Though often referred to as a bulge, both integrated photometry maps and studies 
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of individual stars indicate that the bulge is actually shaped more like a bar. There is evidence 
from the kinematics of M giant stars that the bulge is cylindrically-rotating, with a peak rotation 
velocity of 

90 km/s at 



100 km/s at t ~ 10° ( |Howard et al.[ |2009[ ). At low Galactic latitudes, the velocity 

130 km/s at I 



dispersion of the bulge varies from 



0° to 



±6° (Kunder et al. 



2012). It is difficult to define a total mass for the bulge, because this requires a model-dependent 



separation of this component from the disk. The separation is even more difficult in light of recent 
observations that indicate the bulge is not a classical bulge that formed and evolves separately from 



the disk of the Galaxy (Shen et al. 2010). 



Closer in towards the Galactic center, the MW contains a nuclear star cluster (NSC) originally 
discovered with near-infrared imaging. The MW NSC has a spatial extent of approximately one 
parsec, with a total stellar mass of ~ 10 6 M . Like NSCs at the centers of many other galaxies, 
it contains a mix of old and young stars. Because of its proximity, it is possible to resolve the 
central region of the star cluster, where it is observed to have a cusp with a power law index of 
r -12 within the central 0.22 pc, and a break to a steeper index of r~ 1,7 for larger radii (Schodel 



et al. 2007). Kinematic studies of the MW NSC find that the total mass contained within it can 



be fully-accounted for by the constituent stars of the cluster (Schodel et al. 2009) 



From the astrometric and kinematic measurements of short period stars in the Galactic center, 
there is now conclusive evidence for a ~ 4 x 10 6 M supermassive black hole in the center of the 



Galaxy (Ghez et al. 2008, Gillessen et al. 2009). From measured values of the velocity dispersion 



of stars near the Galactic center, the sphere of influence of the central black hole is < 0.1 pc. 
Beyond this radius and out to a few pc, the population of stars within the central NSC dominate 
the potential near the Galactic center. 

Within the central few parsecs, the contribution from the dark matter and from the mass of 
the disk are sub-dominant. Assuming for the former a smooth dark matter density profile of the 
form predicted in Section [5j the expected integrated dark matter mass is ~ 1000 M Q within the 
central few parsecs. Explicitly including the disk component, the statistical constraints on the dark 
matter distribution in the Galactic center are extraordinarily weak (Iocco et al. 2011). From the 
current observations, it is not possible to tell if the dark matter profile is cuspy or is flat within the 
bulge region. Since the kinematics of this region are well-described by the observed populations of 
stars, and the kinematic data from which these results are derived is itself still relatively sparse, 
these results confirm that it is difficult to disentangle a subdominant dark matter component. At 
the stage, the conservative upper bound on the contribution of the dark matter is set by the upper 
limit deduced from the bulge contribution to the potential. This lack of constraint is somewhat 
unfortunate, because it is reasonable to expect that the brightest source of gamma-rays from dark 
matter annihilation comes from the Galactic center region (Section [6]). 

Theoretical calculations find that the central black hole affects the dark matter distribution 
right around the black hole (Gondolo and Silk 1999). In particular, the adiabatic growth of the 
central black hole may steepen the central cusp of dark matter due to phase space conservation. 
The calculation follows in a similar manner to the enhancement of the stellar density near a point 
mass. Assuming an initial dark matter profile in the Galactic center of the form p oc r~ 7 , it 
can be shown from both conservation of mass and conservation of angular momentum that the 
after the appearance of the black hole, the final mass profile scales as p cx r~ A , where A = 
(9 — 27)/(4 — 7) (Ullio et al. 2001). Using this formalism Gondolo and Silk (1999) derive a lower 
bound of ~ 0.24 GeV cm of dark matter near the Galactic center. Though the theoretical 
framework is compelling, it relies on assumptions that are not well understood, for example the 
adiabatic nature for the growth of the black hole itself (Ullio et al. 2001). Dark matter may also 



interact with the constituents of the NSC and create a shallower profile (Gnedin and Primack 
20041). 
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3.2. Local dark matter in the Solar neighborhood 

As discussed in a historical context in Section [TJ methods for determining the local mass distri- 
bution from the kinematics of disk stars were pioneered by Oort (1932). The methods of Oort were 



revisited and improved upon several decades later (Bahcall 



1984; Kuijken and Gilmore, 1991), and 



there was significant controversy in the interpretation of these results. In particular, it was not 
conclusive whether the determinations of the local potential pointed to Galactic disk that contains 
a significant amount of dark matter, or one that was entirely described by the observed stellar pop- 
ulations (Kuijken and Gilmore, 1991 Bahcall et al. 1992). All of these studies shared the primary 
motivation of searching for a dark component of the Galactic disk- a more spherically-distributed 
mass component was included in the analysis, but the motivation for this component was either 
matter in the form of an extended stellar halo or a dark matter halo as we typically refer to it today. 
The fact that this extended halo component could contribute significantly to the local inventory of 
matter was not the driver for the scientific investigation. 

More recently, there has been a revived interest in the topic of local dark matter, motivated 
by several factors. The most important motivation is that we now have improved experimental 
limits on the dark matter scattering cross section with nuclei (Section [9]), and also limits on the 
annihilation cross section from a variety of indirect detection experiments (Section l7j). Understand- 
ing the local distribution of dark matter is important for interpretation of these results. Further, 
there are now theoretical motivations for the presence of a dark matter disk (Section [5]), that is 
kinematically-distinct from the dark matter halo component. Establishing the presence of this com- 
ponent first requires an understanding of the measurements of the local inventory of dark matter, 
and the systematics, both theoretical and experimental, in its determination. 

Begin by understanding the generic theoretical framework which all of these investigations rely 
on in one form or another. Fundamentally, all determinations of the potential from the kinematics 
of local stars are based upon measurements of the stellar distribution function (DF). This quantity, 
which will be denoted as /, is defined so that there are f(x, v)d 3 xd 3 v stars in a volume of size d 3 x 
centered on position x and in a volume of size d?v centered on velocity v. Though a measurement 
of / is the fundamental goal of stellar dynamical studies, in practice only a fraction of the six phase 
space coordinates are typically accessible from observations. 

For a local population of stars, the accessible phase space coordinates are typically the position 
on the sky (two dimensions), the distance to the star, and the line-of-sight velocity. For the 
primary analysis in this section, assume that the distance to the star, d, along with the Galactic 
latitude and longitude (£, b) are determined. Define a standard coordinate system in which x 
points toward the Galactic center, y points in the direction of Galactic rotation, and z points 
towards the north Galactic pole. Define vi os as the line-of-sight velocity, and ^ as the proper 
motion of the star in the 6-direction. The z-component of the velocity of the star is thus given as 
v z (b,d,vi os ) = v ios s'mb + dn b . 

For the large sample of disk stars that have associated kinematical information, now number- 
ing in the thousands, there are several systematics that must be dealt with when attempting to 
reconstruct the local potential. First is the uncertainty in the kinematic properties that can be 
extracted from observations, because the measurement uncertainties on the line-of-sight velocities 
are typically in the range of 5 — 10 km/s ( jYanny et al. 2009). For cold stellar populations with a 
velocity dispersion of a z < 20 km/s, these uncertainties must be folded into the analysis. Further, 

accurate proper motion measurements are required; for local stars the typical uncertainty from the 

-l 



proper motion is 3 — 4 mas yr (Munn et al. 2004). It is straightforward to see that significant 



noise may be introduced from both the binning of the data and the observational uncertainties in 
the measurements of the positions and velocities. 
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Because the lack of full phase space information and observational uncertainties, it is important 
to consider independent analysis methods, and to understand the strengths and weaknesses of these 
analysis methods. Theoretically the methods can be broadly broken up into two main categories: 
those that directly model / based upon a set of assumptions for it, and those that rely on the 
moments of the DF and thus the more averaged properties of it. 



• Distribution Function-Based Methods 



Begin by considering the direct modeling of the distribution function. Assuming that the 
potential and the density distribution of the Galactic disk are axially-symmetric, the collisionless 
Boltzmann equation is 

d f m df x (w x V "A d f VRV * d f x p df n ^ 

VR M +Vz Tz + [ FR+ R)d^-nrd^ + Fz d^ = > (2) 

where Fr = —d$(R, z)/dR and F z = —d$>(R, z)/dz, and $(R,z) is the total potential for all of 
the mass components that contribute to it. 

As typically implemented, the key assumption of distribution function based methods is that / is 
separable in the coordinates that describe the motion out of the plane, so that / = fity(R, vr, v$)f z {z, v z 
For stars near the Galactic plane, this is a reasonable approximation, because disk stars have low 
orbital eccentricity, and the z-component of the force on a star is nearly constant with R for a 
stellar orbit. 

An implication of this separability assumption is that energy is conserved in the z-direction, so 
that Eg = v z /2 + $(z) is constant. From the definition of the distribution function, the number 
density of stars perpendicular to the plane is then 



oo 



p s {z) = I f z (z,v z )dv z = 2 / 7===, dE z . (3) 

J-oo J<5>{z) y/2(E z - $(») 

An Abel inversion of this equation provides a solution for the distribution function, 

ttJe z V 2 ($( z ) ~ E z) 

Equation [4] is the main equation that is utilized for the distribution function based method. For a 
given stellar population, p s (z) can be determined directly from star counts above the plane (for a 
realistic data set there will of course be uncertainties in the density profile incurred in binning of 
the data). With the stellar distribution in place, all that is required is a model for the potential 
perpendicular to the disk, &(z). As a reminder this is the total disk potential, which includes 
contributions from the stellar population and the dark matter component. 



Kuijken and Gilmore (1989c) introduced a convenient parameterization for the potential, 



2vrG 



K 



(Vz 2 + D 2 -D^ +Fz 2 , (5) 



where K, D, and F are constants. The first term in Equation [5] is designed to represent the 
contribution from a disk with a scaleheight D, while the second term represents the contribution 
from the more extended dark matter halo. These three parameters are determined from the local 
kinematic data. 
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With the model for 3>(z) in place, the total mass density, pr, can be determined from Poisson's 
equation: 

1 d ,„„ 4 . 8F g ' 



pr(R, z) -- 



4vrG 



+ a, 



(6) 



The first term in Equation [6] is related to the circular velocity through V c = Fr/R. In the plane of 



the Galaxy it is safe to neglect this term for plausible potential models of the disk and halo (Kuijken 



and Gilmore, 1989c), however extending up to approximately 4 kpc above the plane neglecting this 



term leads to an estimated 15% error in the determination of pt(R, z) (Bovy and Tremaine 2012). 
So for measurements in the disk midplane, the total volume mass density is directly related to the 
gradient of the potential perpendicular to the plane of the disk. 

The set of parameters K, D, and F can be determined through maximum likelihood or Bayesian 
methods. Equation [4] can be thought of as the conditional probability distribution for E z given 
the set of parameters that determine the potential, p(E z \F, D, K), or explicitly in terms of the 
observable position and velocity, p(z,v z \F,D,K). For the measured density distribution of stars, 
this can be re-expressed as p(z,v z \F,D,K) oc p(v z \z, F, D, K)p s (z). The total probability for the 
set of model parameters {F, D, K} from the measured set of n stars is then 



1 

P(F,D,K\z,v z ) = -P(F)P(D)P(K)l[p(v z , l \z l ,F,D,K)p s (z l ), (7) 

i=i 

where the normalization factor is, assuming that {F,D,K} are independent, 

/n 
dFdDdKP{F)P{D)P{K)]\p{v^\z l , F, D, K)p s (z, t ). (8) 
i=i 

Equation [7] can then be maximized for the parameters {F,D, K} to determine the best-fitting model 
of the potential. Then, by comparing pt derived from Equation [6] and p s as derived directly from 
the data, it is possible to measure the local density of dark matter contributed by the spherically- 
distributed halo component. 

• Moment-based methods 

Though distribution function-based methods are the most direct route to constructing the 
local potential, their accuracy depends on the validity of the assumption for the shape of /. In the 
above case, the critical assumption is that of separability. Complementing the distribution function 
methods, an additional widely used method relies on taking moments of Equation [2] This method 
has the benefit of not depending on the assumptions for the shape of the distribution function. 

Multiplying Equation [2] by v z and integrating over all velocities gives the jeans equation 

d(p s (vRVz}) d{p s {vp) Ps(vrv z ) d$ 

«5 1 a 1 5 h ^^~ = - ( 9 ) 

oR oz R oz 

Two additional independent equations can be obtained by similarly multiplying by vr and v^, 
though for the purposes of the discussion in this section we will be primarily interested in the 
information that is extracted from Equation [9| 

As in the distribution function-based method, in order to implement Equation |9j we must have 
a model for $ and a measurement of p s (z). An additional requirement is a measurement of the 
velocity dispersion perpendicular to the disk, (v z ). This is a key difference between the distribution 
function and moment based methods- in the distribution function based method all that is required 
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is a measurement of the stellar velocities. In addition, from examination of Equation [9j it is clear 
that we need an estimate of the scale of the "tilt" term, (vrv z ), or the correlation between the radial 
velocity in the plane and the velocity perpendicular to the plane. However, from the discussion of, 



for example, Binney and Tremaine (2008) this term is expected to be small. 



At a height z above the disk, the solution to Equation [9] for a single stellar population is 



Ps (z) 



1 d$ 



dz' 



(10) 



Through an iterative procedure, Equation [1Q| combined with Poisson's equation (Equation |6j), allow 
for a solution of the density in the stellar population p s , the total density, pr, and thus the density 
in any component that is not associated to the observed population. 

The calculations above show how the volume mass density is determined over the regime probed 
by the tracer stellar population. A related quantity, that is also derived from the poisson equation 
(Equation [6]) and is often used as a probe of the Galactic potential is the surface mass density. The 
integral surface mass density out to a distance z max from the Galactic plane is defined as 



T,(R,z r , 



Pt(R, z)dz. 



(11) 



Depending on the scaleheight of the stellar population that is used, the surface mass density probes 
the potential on different scales relative to those probed by the volume mass density, so it can be 
viewed as a measurement complementary to the volume density. Note that as defined Equation 11 
is the surface mass density from all components, and not just the component that is associated 
with the stellar population. 

• Application to data sets 

With the theoretical frameworks established we are now in position to understand the results 
obtained from their application to data. In order to reduce the systematics in measurements it is 
best to first define a population of disk stars, and from this population obtain accurate distance 
and kinematic measurements. There are several important properties that an observed population 
must satisfy. From the observational perspective, the population of stars clearly must have both 
accurate distance measurements and must be large enough to obtain a good measurement of the 
density fall-off p s (z). From the theoretical perspective, it is beneficial for the stars to be old enough 
to be well- mixed, which in the context of disk stars means that they have undergone several vertical 
oscillations. 

Different measurements of the local dark matter density have relied on different populations of 



stars. The original measurements of Oort (1932) utilized bright K giants as tracers, as did subse- 



quent measurements of Bahcall (1984) and Kuijken and Gilmore (1989a). Kuijken and Gilmore 



(1989b) utilized a sample of approximately 2000 K giant stars; the distances for this sample were 



re-calibrated by Garbari et al. (2012) in a more recent analysis. Holmberg and Flynn (2004) utilized 



a sample of K giants selected from the Hipparcos survey. Because of their typical age these stars 



have the benefit of being old and dynamically well-mixed. In an earlier study Holmberg and Flynn 



(2000) utilized a sample of A-F stars, which are typically not as old. The recent analysis of Bidin 



et al. (2012 ) utilized a sample of approximately 400 thick disk red giants towards the South Galactic 



Pole. 

The main results of these various analyses are summarized in Table [T] Also indicated in the 
final column is whether the method used by the respective authors most closely resembles the 
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Table 1 . Compilation of measurements of the local surface mass density and the local dark 
matter density. Above the horizontal line, the measurements of the surface mass density and dark 
matter density are derived directly from local measurements of different sets of stellar 
populations. The final column indicates the method used to obtain the results, either the 
distribution function method (DF) or the moment-based method (M). Below the horizontal line, 
the measurements are obtained from a combined analysis including additional measurements of 
Galactic model parameters. Note that 1 GeV cm -3 = 0.027 M Q pc~ 3 . 



Authors 


Total (visible) surface mass density 


Dark matter density Method 




(M pc- 2 ) 


(M pc- 3 ) 



Kuijken and Gilmore| (1991) 



Holmberg and Flyrm J2004 ) 



Bienayme et al.| ( 2006 



Garbari et al.| (|2012[ ) 
Bovy and Tremaine[(|2012|) 



Holmberg and Flynn (2000) 



Catena and Ullio (2010 



Weber and de Boer 



(|2010[) 



Salucci et al. 



McMillan (2011) 



(2010) 



(2010) 



71 ±6 (48 ±! 
74 ± 6 (56) 
64 ± 5 (53) 



0.007 

0.022i°-°J| (90%) 

o.oostffi 



0.011 



+0.0010 

-o.ooio 



[0.005 - 0.010] 
ni2 +0 003 

u,u — 0.003 

011 +o.oon 

U.U11_q oon 



DF 
M 
DF 
DF 
M 
M 



Estimated by subtracting the total mass density from the mass density in visible matter. 



distribution function based or moment based method. Measurements of the total surface mass 
density of the disk out to 1.1 kpc are indicated, along with the contribution to the surface mass 
density from the visible mass that is reported by the respective authors. Note that although the 
visible surface mass density is less than the total surface mass density, this does not mean that the 
difference is contributed by the component associated with the extended dark matter halo. In the 
parameterization of Equation [5j the F parameter most closely parameterizes the contribution from 
the extended dark matter halo. Some of the difference may be contributed by matter associated 
with the disk that is not included in the visible components. As discussed in more detail in Section[5j 
a dark matter disk is predicted to exist from cosmological simulations ( Read et al. 2008 ) , however 



using methods similar to those described above there is not yet kinematical evidence for such a 



structure in the Milky Way (Bidin et al. 2010). 



3.3. The mass of the Milky Way 



Dating back to the original work of Kinman (1959), the velocities of distant globular clusters 



and satellite galaxies have been used as kinematic tracers for the mass in the outer halo. Many 
authors have improved upon this analysis using both better theoretical models and observational 



es ( Little and Tremaine 


1987 


Zaritsky et al. 


1989, 


Kulessa and Lynden-Bell 


1992 


Kochanek 


Wilkinson and Evans 


1999 


Sakamoto et al. 2003 Watkins et al. , 2010 


. More recently, 



1996 



thanks to large scale surveys such as the Sloan Digital Sky Survey (SDSS), it has become possible to 
measure the dark matter halo mass profile using the kinematics of a large sample of Blue Horizontal 



Branch (BHB) stars (Xue et al. 2008 Deason et al. 2012). These BHB stars are particularly useful 
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because they are both intrinsically bright and have accurate distance measurements. 

Similar to the analysis of the dark matter distribution in the Solar neighborhood above, in mea- 
surements of the dark matter halo mass of the Galaxy there are two general theoretical methods 
to consider: those based on assumptions for the distribution function, and those based on the mo- 
ments of the distribution function. In each case, a stellar population is observed, and the kinematic 
information on this population is used to study the underlying properties of the dark matter halo. 
Though similar in this aspect, the analysis in this subsection differs from the Solar neighborhood 
analysis in a couple of different ways. First, in this case not only are stellar populations utilized, but 
also populations of systems of stars, in particular globular clusters and dwarf spheroidals (dSphs). 
Second, the populations in the halo analysis make a negligible contribution to the total potential. 
Although this is not a necessity of the analysis, it is a physical statement that the influence of the 
dark matter halo continually increases at larger radii from the Galactic center. 



• Distribution function based models 



Following the sequence above begin by considering models based on distribution functions. Ob- 
servations considered for this analysis typically provide at most four of six phase space coordinates; 
it is convenient to take these four coordinates as the latitude (b), longitude (I), distance (d), and 
line-of-sight velocity of the object (vi os ). The distribution function is defined as above, though in 
this case we will not make an assumption on the separability of the distribution function. Rather 
break up the distribution function-based models into two broad classes of sub-models: those that 
assume that the velocity distribution for the tracers is isotropic, and those that assume that the 
tracer velocity distribution is anisotropic. The degree of anisotropy is determined through the 
standard parameter, 

(3 = l-a?/a 2 r , (12) 

where a spherical polar coordinate system (r, 6, <p) about the Galactic center is assumed, a r is the 
radial component of the velocity dispersion, and of = Og + cr^ is the sum of the components of the 
velocity dispersion orthogonal to the radial component. The components of the velocity dispersion 
orthogonal to the radial component are taken to be equal, <7g = a^. For (3 < 0, the orbits are more 
circular, while for f3 > 0, the orbits are more radial. For f3 = the orbits are isotropic. 

For the isotropic model, the distribution function depends on a single integral of motion, which 
is the total energy. The mapping between a given density /potential model and the distribution 
function, /(e), is given by the Eddington formula, 

ft \ 1 d [ ^ /I Q\ 

m = 7^deJ m7^> (13) 

where as above p s represents the density of the tracer population. Here \£ is the total potential 
from all of the components, and e = ^/ — v 2 /2, where v is the velocity. 



Equation [13] is used to determine the velocity distribution of a population described by a set of 
theoretical parameters, and it is then to compared to the observed population. For the analysis in 
this sub-section, consider a simplified model for the Galactic potential, which includes contributions 
from the bulge, the disk, and the dark matter halo, so that ^ = + Vt^ + f,j m . It will be most 
straightforward to assume that the potential components are spherically-symmetric. The bulge 
potential can be modeled as 

•.-7^, (14) 
{r + c ) 
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where M b ~ 10 10 M is the total mass of the bulge and cq = 0.6 kpc is its scalelength. Though 
the potential from the bulge near the Galactic center deviates from a spherically-symmetric model, 



Equation 14 will be appropriate at much larger distances. The disk potential can be modeled as 



GM d [l-exp(-r/6 d )] 
r 

where ~ 5 x 10 11 M is the mass of the disk, and bd — 3 kpc is the disk scale length. For 
the dark matter halo, there are a variety of well-motivated models. In this case, we can start from 
the spherically-symmetric density profile. A popular modern model is the Navarro-Frenk- White 
(NFW) profile, 

PNFw{r) = P ° — (16) 

(r/ro)(l + r/r ) 2 

which is motivated by the results from modern numerical simulations (See Section [5] for more of a 
discussion on this profile). Then from poisson's equation, the potential corresponding to the NFW 
profile is 

9 NFW (r) = 4vrGp r 3 ln(1+ r r/r0) . (17) 

At large radius, the sum of the potentials above can be approximated by a pure power law that 
scales as oc r~ 7 . This is primarily because beyond the scale radius the potential is dominated by 
the NFW mass profile. For the tracer density p s it is standard to consider two limiting models. One 
model is described by a pure power law in density, p s oc r~ a , which is appropriate for a stellar halo 



population (Bell et al. 2008) or probably for the dSph population (Evans and Wilkinson, 2000). 
In these cases, a ~ 2 — 4. The second model assumes that the observed stellar population traces 
the underlying dark matter halo, p s oc pdm- 

The final step involved in comparing the theory to the observational data requires transforming 



the distribution function in Equation 13 to a function of observed line-of-sight velocities. This is 
done as 

f(v r] r)= f d 3 v'f(e)S(v' r -v r ), (18) 



where v' 2 = v' 2 + v'q + v' 2 , and r denotes the position in the halo. A similar relation can be derived 
for the tangential motion, though these are typically subject to much larger uncertainties. 

In a manner similar to Equation [7] we can perform a maximum likelihood analysis to determine 
the best fitting parameters that describe the MW dark matter halo. Assuming that the anisotropy 
is fixed, (3 = 0, the likelihood can be maximized for the parameters (0,7), so that 

m 

C(a,j) oc Y[f(v rjl ,f t \a,i), (19) 

i=l 

where the product is over the positions and the velocities of the components of the population. 

The analysis becomes more model-dependent when considering anisotropic distribution function 
models with /3 7^ 0, for the tracer population. In this case, even when a single population contributes 
to the potential, the distribution function no longer depends on a single integral of motion, and there 
is also no longer a one-to-one mapping between the density and the potential. To make the problem 
tractable, it is standard to assume a distribution function that is separable in the energy and the 
angular momentum per unit mass, F(e,L) = g(L)f±(e). If it is assumed that g(L) = L~ 2/3 , then 
the anisotropy (3 for the tracer population is constant as a function of radius, and the distribution 



function splits up into even and odd components as (Watkins et al. 2010) 

F even ocL- 2/ Vi(e), (20) 
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with 

/ 1 ( e )oc £ ^('r- 2 )]/'H-(«/7)-3/2 j (21) 

and 

Fodd oc (1 - 77) tanh(L z /A)F even , (22) 

where L z is the ^-component of the angular momentum, and A is a smoothing parameter. The 
value of Tj determines the amount of rotation: for 77 = there is maximum prograde rotation, for 77 
= 2 there is maximum retrograde rotation. Similar to the likelihood analysis in the case of isotropic 
models, parameters of the distribution may be determined from the data sets. These models form 
a representative set that are typically used to measure the mass profile of the Galaxy from discrete 
populations. 

• Moment-based methods 



Now consider measurements of the Milky Way mass distribution that are based not on param- 
eterizations of the distribution function, but rather on its moments. As above, both isotropic and 
anisotropic models can be considered, and the potential of the Galaxy is assumed to be spherically- 
symmetric, with (r, 9, (/)) the coordinates describing the system. The spherical jeans equation is 



GM(r) 



(23) 



which is derived by taking moments of the Boltzmann Equation, similar to the manner in which 
Equation [9] is derived. 



The velocity dispersions that appear in Equation 23 are with respect to a coordinate system 



centered on the Galaxy. To perform the analysis in this subsection we need to convert the measured 
velocity dispersion, ctgsr, hrto the velocity dispersion in the Galactic coordinate system, a r (as in 
e.g. Dehnen et al. ([2006)). First assume that the Sun is on the z-axis, and 9 is the polar angle. 
The distance to a star is 



d 1 



+ Ri- 2rR G) cos 9. 



The time derivative of this quantity is the velocity in the direction of the star, 



d 



Rq cos 9 
~d 



-R0 sin 9 

v r H v e . 

d 



(24) 



(25) 



The velocity dispersion is then the RMS velocity over all solid angles, (Jgsr 
written as 

<r r y/l-P(r)A{r), 



where 



A(r) 



&GSR 



" 2 + R% 



In 



r + R, 



© 



R, 



© 



(d), and can be 
(26) 



(27) 



8r s R Q 

These equations show that for stars local to the Sun, there is an important correction factor to the 
measured line-of-sight velocity that must be accounted for. For large distances, which in this case 
means r 3> Rq, the correction factor is insignificant, so that cfgsr — °V- 

From the above formalism, Figure [l] shows the velocity dispersion versus radius in comparison 
to three recent observational samples. Shown is the velocity dispersion due to the three primary 
dynamical components of the Galaxy, the bulge, disk, and dark matter halo. The dark matter halo 
clearly becomes the most relevant in the outer regions of the Galaxy. The potential is determined 
assuming that p s oc r -3 ' 5 and that the orbits of the tracers are isotropic. 
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Figure 1 Velocity dispersion versus Galactocentric radius from the data of Xue et al. (2008), Brown 



et al. (2010), and Battaglia et al. (2005). The dot-dashed curve is the velocity dispersion from the 



bulge potential, the dashed curve is from the disk potential, and the dotted curve is from the dark 
matter halo. The solid curve is the sum of all three components. The dark matter halo potential 
is calculated from an NFW profile with po = 6.5 x 10 6 M kpc -3 and ro = 22 kpc in Equation 16 
The velocity dispersion is assumed to be isotropic. 
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Table 2. Compilation of the measurements of the mass of the Milky Way dark matter halo. The 
third column gives the direct measurements within the radius R in the second column using 
kinematic data. The fourth column gives measurements of the total Milky Way mass; in this 



column, all except the Watkins et al. (2010) results rely on an extrapolation using data from 
smaller radii. The fifth column gives the types of tracers that are used: halo stars (S), globular 

clusters (GC), and/or dwarf spheroidals (dSph). The sixth column gives the statistics of the 

orbits of the tracers that are implied from the best-fitting masses, either isotropic (I), circular 
(C), or radial (R). The final column indicates the most closely- resembled method used to obtain 

the results, either the distribution function method (DF) or the moment-based method (M). 
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The Dehnen et al. (2006) results are from a reanalysis of the Battaglia et al. (2005) data. The 



quoted result from Battaglia et al. (2005) is the truncated flat model that gives the best fit. 



Table [2] shows a compilation of recent measurements of the Milky Way mass using both distri- 
bution function and moment-based methods described above. In nearly all of the analyses reported, 
there is a "direct" measurement of the mass on a scale that is directly probed by the tracer popu- 
lation. Nearly all of the measurements of the total halo mass require an extrapolation beyond the 
radius of the visible tracer population. 

For the measurements that use the dSphs as the tracers, there are two systematics that are 
particularly important: first is the unknown density law p s due to the small number of objects that 
are measured, and second is the issue of whether or not Leo I is bound to the Galaxy. Because 
Leo I is moving at a high Galactocentric velocity and is at a Galactocentric radius of 260 kpc, it is 
not yet clear whether this is an outlier that is bound to the halo or if it is unbound and is on its 
first pass through the Galaxy. Most recent proper motions indicate that it may be bound to the 



Galaxy (Sohn et al., 2012 Boylan-Kolchin et al., 2012). 



• Escape Velocity 

The escape velocity is directly related to the mass of the Milky Way dark matter halo. It is 
derived independently from a local population of high velocity stars, providing an important cross 
check on the mass calculations above. For this measurement, it is convenient to parameterize the 
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high velocity tail of the stellar distribution function as 



f(v) oc (v esc - v) 1 



(28) 



This formula can be theoretically-motivated by noting that at high energies the distribution function 
can be parameterized as a power law, /(e) oc e k for isotropic orbits. Once the escape velocity is 
determined, the total potential of the Galaxy at the Solar radius can be derived as 'P total = ~ v esc/^- 
The potential here is the sum of the respective components from the bulge, disk, and halo as defined 



above in Equations 14 



and 15 



Smith et al. (2007) have recently determined the escape velocity by combining this theoretical 
formalism with a maximum likelihood analysis on the sample of stars from the RAVE survey. From 
a sample of a couple of dozen stars they find that the escape velocity lies within the range 498 km/s 
< v esc < 608 km/s at the 90% c.l. There is a strong correlation between v esc and k, so that only a 
lower bound can be determined from the data on these parameters, which is at about v esc = 450 
km/s. Negative values of k are ruled out in the analysis. 

• Full Milky Way Mass Model 

In order to develop a full mass model of the Milky Way, the above measurements of the local 
dark matter distribution and the distribution in the outer halo are supplemented by several other 
sets of observations. 

The first of these parameters is the distance from the Sun to the center of the Galaxy. This 
has been measured recently to high precision using the orbits of short period stars around the 
central black hole, giving R = 8.33 ± 0.35 kpc (IGillessen et al.l [20091) or R = 8.4 ± 0.4 kpc phiz 



et al. 2008). A second measurement is the so-called terminal velocity, which is defined as the peak 



velocity of the interstellar medium along a line-of-sight in the Galactic plane. Relative to the Local 
Standard of Rest (LSR), the terminal velocity is defined as 



v c (Ros'm£) — v c (Rq) sim 



(29) 



where v c is the circular velocity. Peak velocities for several lines-of-sight from surveys have been 



compiled in Malhotra (1994, 1995). A third observation often utilized are Galactic maser sources. 



High precision trigonometric parallaxes and proper motions of masers in high-mass star-forming 
regions across the Milky Way have been determined, providing three-dimensional location and 



velocity vectors relative to the Sun. Reid et al. (2009) provide the measurements of 18 proper 



motions and line-of-sight velocities for masers. A fourth observation is the proper motion of Sgr 
A*, the source at the center of the Galaxy. Since Sgr A* is fixed at the Galactic center to within 
approximately 1 km/s, its measured proper motion stems primarily from the motion of the Sun 
around the Galactic center, which is due both to the motion of the LSR and the peculiar motion 
of the Sun relative to the LSR. 

By combining all of these observations, several authors have recently reported measurements 
of the local dark matter density from the ensemble of data. Since the model parameterizations 
used in each of the analyses are somewhat different, it is important to clarify the differences in the 



assumptions. Catena and Ullio (2010) use two different parameterizations of the dark matter halo, 



the first is the NFW profile and the second is the Burkert model, defined by 



p(r) 



Pb 



{l+r/r b )(l + r/r b ) 2 ' 



(30) 



Their motivation was to determine whether the shape of the density profile affects measurement 



of the local dark matter density. From their Bayesian modeling, Catena and Ullio (2010) conclude 
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that for an NFW profile the local dark matter density is p@ = 0.389 ± 0.025 GeV cm -3 . The mean 
of this result is confirmed by McMillan (2011), who finds Pq = 0.40 ± 0.040 GeV cm -3 also using 
an NFW dark matter profile. 

Though several authors have converged on similar values for the mean dark matter density 
with the NFW halo based analysis, it is important to dissect this further to see if there are any 
systematics that may affect these results. The first obvious one is the assumption of spherical 
symmetry; breaking from this assumption makes the problem theoretically much less tractable. As 
discussed at the end of this section there is some evidence that the dark matter halo is aspherical. 



An additional systematic that was brought to light in Weber and de Boer (2010) is that there is 



a degeneracy between the scale radius of the Milky Way disk and the dark matter halo. As an 
example, using a generalized parameterization of the NFW dark matter density profile 



pir) 



Po 



c a (l + x b )( c - a )/ fe ' 



(31) 



where r/ro, there is a strong correlation between the local dark matter density and ro- 

In sum, all of the results presented above and summarized in Table [l] and Table [2] show that 
astronomical measurements are providing interesting constraints on the local dark matter density 
and the Milky Way halo mass, though significant systematics still remain in these measurements. 



3.4- Shape of the Milky Way halo 

The discussion to this point has focused on spherical models for the Milky Way potential. The 
motivation for this focus was twofold. On the one hand, spherical models are the most straight- 
forward to understand and interpret theoretically. Extending to non-spherical potentials certainly 
enlarges the model parameter space, and it makes the comparison between different theoretical 
models often less straightforward. A second motivation for restricting to spherical models is that 
the data sets considered to this point are well-described by them, meaning that they typically do 
not require more complexity. This section on the observational properties of the Milky Way dark 
matter halo closes by examining sets of observations that are, in principle at least, able to determine 
if the Milky Way potential deviates from a purely spherical model, and examines the state of these 
observations. 

There are three separate sets of observations that can determine the shape of the Milky Way 
dark matter halo. The first relies on the measurements of the velocity ellipsoid from different sets 
of stellar populations in the Galaxy. The second relies on the properties of the gaseous HI disk. 
The third relies on the spatial and velocity distribution of disrupted satellite galaxies. 

In order to understand the analysis discussed below, it will be first necessary to define a set 
of ellipsoidal coordinates, and relate them to the spherical coordinates utilized above. For a given 
Galactic model, the elongation of the potential will generally be different than the elongation of 
the corresponding density distribution. The ellipsoidal coordinates, either for the potential or 
for the density, can be defined by transforming the spherical coordinate r — > m, where m 2 = 
x 2 /a 2 + y 2 /b 2 + z 2 /c 2 . The ellipsoid is thus characterized by the axes a, b, c, and from these define 
the axis ratios q = c/a and s = b/a. From these ratios we can define the three classes of ellipsoids, 
which are prolate, a > b ~ c, oblate, a ~ b > c, and triaxial, a > b > c. A spherical halo has 
a = b = c. 

For the purposes of the discussion below it will also be necessary to define the standard spherical 
polar coordinate system, in which the radial velocity coordinate v r points away form the Galactic 
center, vg is the velocity in the direction of the zenith angle measured from the North Galactic 
pole, and v<j> is the velocity in the azimuthal direction opposite of Galactic rotation. 
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Shape of the Velocity Ellipsoid 



Now consider in more detail measurements of the shape of the halo from the velocity ellipsoid. 
For this measurement, begin by considering a sample of stellar velocities from an observed popula- 
tion of stars. From the sample of stellar velocities, the velocity dispersion tensor along any set of 
orthogonal coordinates is defined as 



i j 



((«» - i^))( v j ~ i v j)))- 



(32) 



As above the average is over the phase space distribution function, i.e. (vi) = f v t f(x, v)d?v. 
Since the velocity dispersion tensor is a second-rank tensor that is symmetric, it can always be 
diagonalized, and the resulting set of principal axes form the velocity ellipsoid. 
For any pair of coordinates, the tilt of the velocity ellipsoid is defined as 



tan(2a l:/ ) 



of 



(33) 



.1,1 



Defined in this manner, the tilt represents the angle between the i th the major axis of the ellipse 
formed by projecting the velocity ellipsoid onto the zj-plane. 

The tilt angles in Equation [33] have important theoretical properties. It has been long known 



that if all of the tilt angles are zero, then the potential is spherically-symmetric (Eddington, 1915 



Lynden-Bell, 1962). More generally, if the velocity ellipsoid is radially aligned everywhere, then 
the underlying potential of the system is spherically-symmetric (Smith et al. 2009). Even if the 



velocity ellipsoid is aligned along the spherical polar coordinates and the potential of the dark 
matter halo is spherical, the density distribution of the population that is used to determine the 
velocity ellipsoid itself need not be spherical. 

Thus the measurements of the velocity ellipsoid provides a powerful probe of the potential of 
the Galaxy, and thereby the shape of the dark matter halo. Clearly the observational challenge is to 
obtain a large population of stars that have full six-dimensional phase space information associated 
with them, so that the shape of the potential may be determined by finding the tilt angles via 
Equation 33 Further, since we are particularly interested in understanding the properties of the 



outer region of the Galactic halo, it is most beneficial to obtain a sample of stars that are outside 
of the disk, where the potential is dominated by the dark matter halo. 



Several recent studies have measured the tilt of the velocity ellipsoid. Siebert et al. ( 2008 ) used 



the phase space information from a sample of 763 red clump giant stars from the RAVE survey that 
are distributed between distances of 500 to 1500 pc below the Galactic plane. From their sample 
they measured tilt angle of auw = 7.3 ± 1.8° at ~ 1 kpc below the plane. Because of the sensitivity 
to the disk potential at this radius, the fact that the tilt is different than zero does not imply that 
the potential of the dark matter halo is non-spherical. Using an axisymmetric NFW dark matter 



halo model, Siebert et al. (2008) show that the observed tilt may be explained by either a prolate 
dark matter halo with a disk scale length of ~ 3 kpc, or an oblate dark matter halo with a disk scale 
length of ~ 2 kpc. These results highlight the degeneracies between the halo and disk potential 
even for a fixed NFW model. 



Smith et al. (2009) utilize a sample of approximately 1800 sub dwarfs in SDSS stripe 82 in 250 



deg^ to measure the tilt of the velocity ellipsoid in all three of the spherical coordinates (r,v,(p) 
defined above. Since these stars extend out to a distance of 5 kpc from the disk, they are a more 
ideal sample to probe the dark matter potential than the above RAVE sample. Interestingly, they 
show that each of the angles a r( p and ag^ are consistent with zero, and that the third angle a r g 
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is consistent with zero at the three-sigma level, and its deviation from zero can be explained as a 
perturbation from the Galactic disk potential. 

Ideally an even more extended sample of stars would provide the most robust measurement of 
the shape of the halo; the sample of halo stars used above to determine the Milky Way circular 
velocity out to 60 kpc particularly interesting. A measurement of the kinematics of this population 
will be possible with the GAIA satellite, which will measure the proper motions of bright halo 
stars with an uncertainty < 1 mas. It is only possible at present to determine the shape of 
the distribution of the halo stars, and because they most likely formed via a different physical 
process than did the underlying dark matter halo, a direct measure of their density distribution 
cannot reliability constrain the shape of the dark halo. Studies of halo stars out to approximately 
35 kpc using the Canada-France-Hawaii Telescope indicate the oblateness of the stellar halo is 
c/ a = 0.70±0.01 ( |Sesaret al.}|2011[ ). The best-fitting Einasto fit to the stellar density is n = 2.2±0.2 
and r e = 22.2 ± 0.4 kpc, which is less steep than the dark matter as is predicted in the LCDM 
theory of structure formation (Section [5]). 

• Thickness of HI gas disk 

The thickness of the HI gaseous disk of the Milky Way is observed to increase with radius, 



with an extent of 0.1 kpc at approximately 5 kpc to 1 kpc at approximately 25 kpc (Wouterloot 



et al. 1990). Since the HI gas disk is more extended than the stellar disk, it is more sensitive to 



the potential from the dark matter in the outer halo than the stellar disk, and it is interesting to 
determine whether dark matter halo models can reproduce the above observed "flaring" of the HI 
disk. Several studies have shown that it is difficult to reproduce the flaring within the context of 



spherical halo models (Narayan et al. 2005; Kalberla et al. 2007). Banerjee and Jog (2011) find 



that the flaring can be reproduced by a prolate dark matter halo that becomes the most prolate at 
~ 25 kpc, q ~ 2. 

• Dynamics of Stellar Streams 



The dynamics of stellar streams also provide a probe of the shape of the Galactic dark matter 
halo. In particular, the positions and the velocities of many of the stellar components of the 
Sagittarius stream are now well-determined from a combination of photometry from SDSS, 2MASS 
and high resolution spectroscopy. Assuming the parameterizations of the bulge, disk, and halo 
discussed above, however, it has proven to be difficult to obtain a fit to the detailed properties 
of the Sagittarius stream. Two particular difficult aspects to fit are first what appears to be a 
bifurcation of the leading stream, and second the position of the stream on the sky. Assuming that 



the bifurcation is not related to the physical dynamics of the stream, Law and Majewski (2010) 



find that the best-fitting dark matter halo is triaxial, which for the potential give q = c/a = 0.72 
and s = b/a = 0.99. For the halo density the respective axis ratios are q = c/a = 0.44 and 
s = b/a = 0.97. This corresponds to a nearly oblate ellipsoid with the minor axis contained 
in the Galactic plane. The Sun does not lie along a principal axis, but rather is rotated about 
approximately 10 degrees from the minor axis. 

Because of the difficulty in matching the Sagittarius stream to a halo model, and because in 
future data sets fainter streams are likely to be more prevalent, it is also possible to utilize fainter 



and thinner streams to study the shape of the halo. As an example, Lux et al. (2012) show that 



the globular cluster tidal stream NGC 5466 favors an oblate or triaxial halo due to a feature on 
the westward side of the stream. Orbital integrations and theoretical modeling designed to match 
observations of faint streams and globular clusters will become more important in the coming years. 
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In sum, through several complementary observational techniques it is now possible to determine 
whether the Milky Way dark matter halo is non-spherical, and even if the halo is prolate, oblate, 
or triaxial. As discussed each of these methods have their respective strengths and weaknesses, 
and there are different lines of evidence that point towards spherical, prolate, oblate, and traixial 
models. In addition to the potential for improving these results with future large scale surveys, 
a variety of other interesting methods may become available to complement these methods, for 



example via the use of precision astrometry of hypervelocity stars (Gnedin et al. 2005). 
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4. Dark matter in Milky Way satellite galaxies 



The first galaxy ultimately identified as a dwarf satellite of the Milky Way was discovered 
by Harlow Shapley in the constellation Sculptor (Shapley 1938a). Initially Shapley classified the 
system in Sculptor as a cluster of galaxies, however follow-up observations revealed that the con- 
stituents of the system were not galaxies, but were in fact individual stars, and the distance to the 
system was approximately 80 kpc. At the time, the extent of this system, and the fact that the 
stars were individually-resolvable, made the system more akin to a globular cluster than a known 
type of galaxy. Further observations by Shapley revealed a similar organization of stars in the 
constellation Fornax (Shapley, 1938b). 

Fast forwarding up to the end of the 20th century, seven more systems with similar properties 
to Sculptor and Fornax were discovered. Those discovered prior to 1990 (Draco, Ursa Minor, Leo 
I, Leo II, Carina) were located by visual inspection of photographic plate data. Sextans, which was 
discovered in 1990, was the first satellite discovered via automated plate machine scans. In 1994, the 
Sagittarius galaxy was serendipitously discovered during the course of a survey of the outer regions 
of the Galactic bulge (Ibata et al. , 1994). Though this system is more luminous than Fornax, it was 
not located earlier because it lies in a region of high foreground density behind the Galactic center 
(For a more detailed review of the discovery of the Milky Way satellites, see Willman (2010)). 

These nine systems have come to be classified as dwarf Spheroidal (dSph) galaxies, and for 
reference from this point forward these nine dSphs will be called "classical" satellites . The classical 
dSphs range in Galactocentric distance from approximately 15 — 250 kpc, and their overall distance 
distribution in the Galactic halo is much more extended than the more centrally-concentrated 
Globular cluster population. The two brightest of the Milky Way satellite galaxies, the Large and 
the Small Magellanic Clouds, are of course visible to the naked eye and were known long before 
the discoveries of Shapley. These galaxies are also distinct from dSphs in that they contain HI gas, 
whereas the dSphs are devoid of gas up to the present observational limits. 

At the turn of the century, it was certainly reasonable to assume that our census of Milky 
Way satellites was incomplete. However, the nature of the incompleteness was unknown. From 
a systematic search of the Southern sky, Kleyna et al. (1997) reported that the census of Milky 
Way satellites is likely complete down to approximately 1/8 of the luminosity of Sculptor, or 
approximately a few times 10 5 L . Willman et al. (2004) estimated that there may be a few times 
more dSphs at the luminosity of Fornax that may have gone undetected in previous surveys due 
to systematic effects such as disk obscuration. Whatever the nature of the incompleteness, it was 
clear that large scale galaxy surveys coming online would be able to find more faint satellites, if 
they indeed were to exist. 

The Sloan Digital Sky Survey (SDSS) is a photometric and spectroscopic survey of ~ 1/5 of the 
northern Galactic sky. The SDSS is spectroscopically-complete from objects brighter than a mag- 
nitude r = 17.7; since this is brighter than the most luminous stars in any of the classical dSphs, it 
was expected that SDSS would be sensitive to the detection of new satellites as resolved overdensi- 
ties of stars through imaging data. As of the publication of this article, searches of the SDSS data 



for overdensities of stars in color space have revealed a dozen new Milky Way satellites (Willman 



et al. 2005 Belokurov et al. 2007b). The distances to these objects are determined by comparing 



the population of stars to old stellar populations in globular clusters, and determining the position 
of the turn-off of the stellar main sequence. 

From deeper photometric study of these objects, it was revealed that the satellites discovered 
in the SDSS are unique in many ways from the classical satellites. For starters, they all have 
total luminosities less than the faintest known classical satellite, so it indeed turns out that, at 
least in the part of the sky surveyed in detail by SDSS and for objects with surface brightnesses 



30 



-15 



Dwarf spheroidals 



Globular clusters 



-10 



QJ 

3 
2 

bo 




10° 10 1 10 2 

Half-light radius (pc) 



10 3 



Figure 2 Relationship between the magnitude and the half-light radius for Milky Way globular 
clusters (red squares, left) and dwarf spheroidals (blue triangles, right). The ultra-faint satellites 
appear here as the blue triangles in the lower middle portion of this figure. Globular cluster data 



from Harris (1996). 



characteristic of the classical satellites, the census of classical satellites is complete. Examination of 
their metallicities indicates internal spreads that are more characteristic of galaxies that have had 
multiple episodes of star formation, as opposed to globular clusters, which largely formed its stellar 
population in a single burst. Further, several of these have half-light radii in the range ~ 30 — 100 
pc, i.e. in between the respective half-light radii of typical globular clusters. This exceedingly 
faint population with intermediate half-light radii can only be detected out to distances of ~ 50 
kpc, a small volume of the Galactic halo. The population of satellites discovered in the SDSS 
are typically referred to as "ultra-faint" satellites. Figure [2] shows the magnitude of the known 
population of globular clusters and dSphs as a function of their half-light radii. In this figure the 
ultra-faint satellites appear towards the bottom middle and right of the figure, while the classical 
dSphs appear largely in the upper right. Figure [3] shows a sky distribution in Galactic coordinates 
of the known Milky Way satellites; the satellites discovered in the SDSS are clearly clustered in the 
northern Galactic cap. 

The most relevant question about both the classical and ultra-faint satellites for this article is, 
do they contain dark matter, and if so, how much? About forty years passed from the discovery 
of the first satellites by Shapley before spectrographs were sensitive enough to measure velocities 
for constituent stars and study their kinematic properties. The first measurements of the velocities 
of stars in Draco and Ursa Minor were taken by Aaronson (1983). Even from his measurements 



of the velocities of only a few stars the Aaronson (1983) observations hinted that these objects 



have total mass-to-light ratios of a hundred or greater, which if true would make them the most 
dark matter dominated objects known. Improved statistics on the sample of stellar velocities in the 
dSphs was clearly needed. In the following decade, the sample of velocities was improved upon by 



Mateo et al. (1993); these observations also indicated that the dark matter mass-to- light ratios of 



the dSphs was of order a hundred or greater, and there were additional intriguing hints that these 
objects had very similar masses. 
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Figure 3 Aitoff projection in Galactic coordinates of the globular clusters (squares, red) and dwarf 
spheroidals (triangles, blue). 



Concurrent with the early determinations of the dSph masses, it was recognized that they 



strongly constrain the nature of particle dark matter (Gerhard and Spergel, 1992). Because of 
the extreme total mass-to-light ratios, it is unlikely that the luminous material has altered the 
distribution of dark matter in these systems, making them important testbeds for dark matter 



models (Hogan and Dalcanton, 2000). In the following years, significant theoretical attention has 
been paid to the dSphs as perhaps the best astrophysical probe into the nature of particular dark 
matter. 

Modern kinematic data sets, which utilize multi-slit spectroscopy, now measure the velocities 
of hundreds to thousands of stars in both the classical and ultra-faint satellites; theoretical in- 
terpretation of these modern results is a main focus of this section. The theoretical framework 
in this section provides the basis for the analysis of satellites with indirect dark matter detection 
experiments, as discussed in Section [7} 

4-1- Theoretical modeling 

The dark matter mass distributions of the dSphs are determined via methods similar to those 
discussed in Section [3| In this case, the observations of bright stars in these systems provide 
three out of six phase space coordinates; these three coordinates are the two-dimensional projected 
position of the stars and the line-of-sight velocity from the stellar spectra. Similar to the analysis 
in Section [3j reconstruction of the full stellar distribution function is difficult because of the lack 
of complete phase space information. It is therefore standard to use moments of the distribution 
function to analyze these observations. 

• Moment-based methods 
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Assuming that the stellar distribution and the dark matter distribution are spherically-symmetric, 
the analysis begins with the jeans equation, Equation 23 In this case, 
dimensional density profile of the stellar distribution. 



Ps 



represents the three- 
The mass profile, M(r), is the sum of the 



contribution from the dark matter and the stars. The three-dimensional density profile of the stellar 
distribution is related to the observed projected distribution, I S (R), via 



Ps{r) 



1 f°° dI s (R) dR 



7T 



dR 



(34) 



The line-of-sight velocity of a star in a dSph is a mixture of the radial and tangential velocity 
components, and therefore depends on the intrinsic velocity anisotropy of the stars (Equation 12). 
At a projected position R, the line-of-sight velocity dispersion is 
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(35) 



The three-dimensional radial velocity dispersion a r is determined from Equation 23. From the 
examination of Equation [35] it is manifest that there is a degeneracy between the mass profile and 
the velocity anisotropy, /3(r), if the mass profile is determined only from line-of-sight velocities. 

Because of the assumption of spherical symmetry, there is no equivalent separability condition 
that led to Equation [4j so the stellar distribution function in dSphs does not simplify in the 
same manner as the distribution function of local disk stars. Perhaps the simplest theoretical 
approximation that can be made for the intrinsic line-of-sight velocity distribution is that it is 
gaussian with a variance given by ai os . Define the true line-of-sight velocity of a star as vt rU e, the 
measured velocity as vi os , the mean line-of-sight velocity as v, and the measurement uncertainty 
on a star to be a m . Take the measurement uncertainties to be gaussian-distributed. In this model, 
the line-of-sight velocity distribution for a star at a projected position R is 



f(vi os ,R) 
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(36) 



In Equation 36, af reflects the intrinsic variance, as derived from the distribution function. Here 



no rotational component is assumed; the mean velocity v simply reflects the systemic motion of the 
dSph. The assumption of a small organized rotational component is valid because there are strong 
upper limits on the rotational velocity of dSphs, < 1 km/s, which implies that they are dominated 
by disorganized motion at all radii. 

Similar to the calculation presented in Section [3j the best-fitting parameters for the mass dis- 
tribution are determined by maximizing the likelihood function, 



C{9) = \\f(v los ^R l 



(37) 



where n is the number of stars with line-of-sight velocities, and 6 represents the set of model 
parameters that describe the mass distribution. 

The approach outlined above is actually a hybrid method that utilizes both the moment-based 
modeling through the spherical jeans equation combined with an assumption for the gaussian nature 
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line-of-sight velocity distribution. The observed velocity distribution function is a convolution of 
the intrinsic distribution with the distribution of measurement uncertainties. If the latter are 
gaussian and are of similar order to the intrinsic velocity dispersion then even if the true form 
of the distribution function differs from the gaussian model, the method outlined above is still 
likely to provide an accurate model for the distribution function, because it is dominated by the 
measurement uncertainties. This is the appropriate regime of parameter space for the ultra-faint 
satellites, as discussed below. However, if the measurement errors are small relative to the intrinsic 
velocity dispersion, as is the case with the classical dSphs, the line-of-sight velocity distribution 



may deviate from the gaussian model in Equation 36, so the approximation made in constructing 



Equation [36] must be more rigorously tested. 

A related point is that for an assumed parameterization of the velocity anisotropy profile, the 
potential, and the stellar density distribution, there is no guarantee that a unique physical stellar 
distribution function exists. If a physical distribution function does not exist, then there may be a 
bias that is incurred when determining the mass, and other physical properties of the dark matter 
distribution. This provides an additional important motivation for testing the method outlined 
above with self-consisent distribution function models. 

It is possible to relax the assumption of spherical symmetry in the context of moment-based 
models. Indeed, in many ways non-spherical models are more appropriate, because despite their 
name, the dSphs are not spherical; typical ratios of the minor-to- major axes from the photometry 



are approximately 0.3 (Irwin and Hatzidimitriou, 1995). However non-spherical models do introduce 



a larger deal of complexity, and it is easy to lose sight of physical motivation for the parameters that 
describe them. For this reason, at this stage it is best to provide a basic outline for the ingredients 
of a non-spherical model for comparison to the spherical case. 

Take (R, (ft, z) to be the standard cylindrical coordinates, and assume that the distribution 
function depends on the energy and the z-component of the angular momentum. With these 
assumptions (v^) = (v^). From the axisymmetric jeans equations, the solutions for the velocity 
dispersions are 

1 f°° <9<E> 

^ = mW 2 dz ' p ^ ] ^ (38) 

{v *> ~ {Vr) + Ps (R,z) dR +R dR (39) 

In order to compare to the measured dispersion profile, the (R, (ft, z) coordinates must be converted 
to coordinates along the line-of-sight. The differences between the dark matter halo masses in 



spherical and certain non-spherical models using this modeling have been discussed in Hayashi and 



Chiba (2012). 



Distribution function-based models 



The first modern set of self-consistent models for the dSphs were developed by Wilkinson et al 



(2002). These authors started with the premise that the stellar density distribution was described 
by a Plummer model, which is generally in good agreement with the measured stellar density 
profiles. The surface brightness for a Plummer model is 

4 PpiflVpi 
3 [1 + {R/r pl f 



T * w# (40) 
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and the corresponding density is 
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PAr) = - , /2 - (41) 



The Plummer model can then be fit to a dSph to obtain the Plummer radius, r p \. With the 



parameters of the Plummer profile fixed for a given dSph, Wilkinson et al. (2002) developed a 
model for the distribution function specified by two parameters: one that controls the shape of the 
potential, a, which is defined via 

„ > J(l + r 2 ) Q/2 if a^O 

Mr) <* S ; h ! 2 N n ( 42 ) 

[log (1 + r z ) if a = 0, 

and a second parameter that controls the anisotropy of the stellar velocity distribution. The param- 
eters that describe the potential were chosen to vary between a pure mass-follows-light model (i.e. 
no dark matter) and a dark matter-dominated model with a harmonic potential that corresponds to 



a constant density core. From the solution to Eddington's formula (Equation 13), a self-consistent 
model for the velocity distribution function is obtained. 

Though these models are interesting because they allow for an analytic calculation of the distri- 
bution function, there is room for theoretical improvement. First, if anisotropic models are consid- 
ered, the radial dependence of the velocity anisotropy is restricted in its functional dependence, so 
that the velocity anisotropy is an increasing function of radius. From numerical simulations how- 
ever it is more probable that the velocity distribution of both the stars and dark matter are more 
circular because of the stripping of radial orbits due to interactions with the Galactic tidal field. 
Second, the Plummer profile is probably too restrictive of a model for the light profile, because it 
only allows for a core in the stellar distribution in three dimensions, and precludes the existence of 
cusped stellar density models. 

Though there is clear motivation for expanding the parameter space of distribution function- 
based models, because many are non-analytic it is difficult to gain a physical intuition for more 
general models of the potential and the stellar density. For example there is clear motivation to 
analyze both the photometry and the kinematics of the dSphs independently, in particular in order 
to allow for cusped central profiles of the stellar density. Though the distribution functions are 
typically non-analytic, it is straightforward to obtain it from a numerical solution to Eddington's 



formula (Equation 13) 



For both the moment-based and distribution function based modeling above, a functional pa- 
rameterization for the velocity anisotropy is assumed, which ultimately leads to a degree of model 
dependence. Recently methods have been developed that bypass this issue and do not assume a 
form for the velocity anisotropy, but rather determine it directly in the analysis. This numerical 



method, originally developed by Schwarzschild (1979), which has long been used to determine the 



masses of central black holes in more massive galaxies, is now being applied to discrete velocity 
samples of dSphs. 

In the context of dSphs, very broadly Schwarzschild modeling begins by assuming a form for the 
dark matter potential that is axially-symmetric to match the observed stellar distribution. With 
the assumed potential, orbits are launched and integrated for a large number of crossing times. 
Each orbit is then assigned a weight, and a set of weights are chosen so that the photometry 
and kinematics are matched. Because the orbits are directly determined in this method, the 
velocity anisotropy for the best-fitting model is automatically determined. On the down side, 
the computational time required to implement these models for the large modern data samples is 
significant. 

In sum, there are now several theoretical methods that have been developed to analyze the 
kinematic data from dSphs. These models are very broadly classified into the following general 
categories: 1) moment-based methods in combination with a likelihood analysis described by Equa- 
tion 36 , and 2) distribution function models with either empirically or theoretically-motivated forms 
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for the stellar density, potential, and velocity anisotropy. Below, these methods are compared to 
the most basic estimates for the mass obtained from the virial theorem. It is of course true that 
the dark matter masses as determined from each of these methods must be mutually consistent. 



4-2. Classical satellites 

Classical Milky Way satellites are more extended and luminous than their ultra-faint counter- 
parts discussed below. The half-light radii for these systems are typically of the order few hundreds 
of parsecs, and their luminosities spread over a range of nearly two orders of magnitude, approxi- 
mately 10 5 — 10 7 Lq. The photometric profiles as derived from star counts are typically consistent 
with a cored model in projection, followed by a turnover into an exponential fall-off in the outer re- 
gion of the galaxy where the stellar density blends into the background star counts. For the classical 
dSphs there are hundreds, and in some cases thousands, of bright giant stars that have measured 
velocities to a precision of a few km/s or less (Walker et al. 2009). The distribution of stars with 
line-of-sight velocity measurements does not necessarily follow the same density distribution as the 
integrated star counts. 

An observed sample of velocities must be cleaned of interloping stars that are in the direction 
of the dSph, but are not physically associated with it. Several authors have discussed methods for 
cleaning interloping stars. Klimentowski et al. (2007) utilize a method, calibrated to simulations, 
that rejects stars with a velocity greater than a maximum velocity available to a star at a given 
projected radius. The maximum velocity is determined either by assuming that the star is on 
a circular orbit or is infalling freely into the galaxy potential. Martinez et al. (2011) develop a 
Bayesian method that may be applied to classical satellites that assigned a membership probability 
based on metallicity, distance, and velocity. 

Applying the above selection criteria to the data samples, a maximum likelihood analysis that 
utilizes Equation|36|provides the integrated dark matter mass within a fixed physical radius r, for an 



assumed parameterization for the dark matter density profile (Strigari et al. 2008; Martinez et al. 



2009[ I Walker et al.[ |2009[ ). Several analysis methods have shown that the line-of-sight velocity 



data sets constrain the integrated dark matter mass within the approximate half-light radius, 



independent of the parameterization of the velocity anisotropy (Strigari et al. 



2009, 


Wolf et al. 


2010) 
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where here (<x 2 



lost 



is the luminosity weighted velocity dispersion. Because each of the dSphs have 
different half-light radii, the mass distribution is best constrained at different fixed physical radii. 
Table [3] shows some properties of many important satellites whose masses can be readily derived. 

Though the mass at the half-light radius does clearly indicate that the dSphs are dark matter- 
dominated and provides a measurement of the central dark mater density, it does not give a handle 
on the total mass of their respective dark matter halos. Unfortunately it is difficult to determine 
these total halo masses, primarily because it is not certain how far the halos extend beyond the 
observed stellar distributions. Therefore, any estimate of the halo mass necessarily introduces a 
significant degree of model dependence. Strigari et al. (2010) examined isotropic velocity models, 
utilizing the light profiles, kinematic data, and dark matter potentials that are predicted in LCDM 
models (Section]!]) to show that the total mass distributions of the brightest satellites are consistent 
with being in the range 10 8 — 10 9 M . Figure [4] shows the velocity dispersion profiles and the 
best-fitting photometric and kinematic models that correspond to these dark matter halo masses. 
These results are broadly consistent with other determinations of the mass distribution that utilize 
moment based methods (Lokas 2009). 
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Table 3. Compilation of properties of Milky Way satellites. Classical satellites are above the 
horizontal line, while ultra-faint satellites are below the line. This table only shows dSphs with 
well-measured kinematic data, in which cases the mass can be determined using Equation 43 In 



the fifth column the half-light radii are reproduced from McConnachie (2012). In the sixth 
column the luminosity weighted velocity dispersions are reproduced from Wolf et al. (2010). 



Name 


Distance [kpc] 


Magnitude 


Half-light radius (r^) (pc) 


K s ) (km/s) 


Sculptor 


80 


-9.8 


283 


9.0 


Fornax 


138 


-13.1 


710 


10.7 


Leo I 


250 


-11.9 


251 


9.0 


Leo II 


205 


-10.1 


176 


6.6 


Draco 


80 


-9.4 


221 


10.1 


Ursa Minor 


66 


-8.9 


181 


11.5 


Carina 


101 


-9.4 


250 


6.4 


Sextans 


86 


-9.5 


695 


7.1 


Willman 1 


38 


-2.7 


25 


4.0 


Ursa Major I 


106 


-5.5 


319 


7.6 


Ursa Major II 


32 


-4.2 


149 


6.7 


Hercules 


138 


-6.6 


330 


5.1 


Leo IV 


158 


-5.0 


206 


3.3 


Canes Venatici I 


224 


-8.6 


564 


7.6 


Canes Venatici II 


151 


-4.9 


74 


4.6 


Coma Berenices 


44 


-4.1 


77 


4.6 


Segue 1 


2006 


-1.5 


29 


4.3 


Bootes I 


2004 


-6.3 


242 


9.0 
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Figure 4 Photometric profiles (left) and velocity dispersion profiles (right) for five classical dSphs, 



using LCDM-based models for the dark matter potentials. From Strigari et al. (2010) 



Schwarszchild mass estimates have now been published for three dSphs, Fornax, Sculptor, and 



Draco. Using a cored model for the stellar light profile, Jardel and Gebhardt (2012) find a mass 



within the half-light radius that is consistent with those deduced from moment and distribution 



function-based methods. Breddels et al. (2012) determine that the mass of Sculptor within 1 kpc 
10 8 Mq, which is again in agreement with the above methods. Jardel et al. (2012) determine 



is 



a lower bound to the mass of Draco of a few times 10 8 Mq within a physical radius of about 500 
pc where kinematics of stars are measured. 

4-3. Ultra-faint satellites 

Measuring the velocity dispersion, and thus the mass, of ultra-faint satellites poses different 
sets of challenges in comparison to measuring the velocity dispersions of classical satellites. First, 
by their very nature, the constituent stars are fainter, with typical target stars having a magnitude 
of r = 20 — 21. For a realistic exposure level, the Keck/DEIMOS spectrograph provides a signal- 
to-noise on a star of this magnitude of approximately 15 (Simon and Geha, 2007). Second, the 
measured uncertainties derived from the stellar spectra are approximately 2 — 3 km/s, which, 
because it is within about a factor of two of the intrinsic velocity dispersions of the systems, 
complicates the extraction of the intrinsic velocity dispersion that arises from the distribution 



function (Low velocity dispersions at this level have been measured in globular clusters (Cote 
2002)). Third, the measured line-of-sight velocities may contain a component that is due 



et al. 



to the motion of the star around a binary companion. Early studies of classical dSphs indicated 
that this binary contamination was ~ 1 — 2 km/s, so it is a small systematic to classical satellite 
mass measurements (Kleyna et al. 20021). However, as measurements of low velocity dispersion 



globular clusters clearly indicate that binaries do significantly contaminate the velocity dispersion 
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of bound stellar systems (Bradford et al. , 2011 ), detailed understanding of this effect is required in 
ultra-faint satellites. 

Another interesting difference between the ultra-faint and the classical satellites relates to how 
their total luminosities are derived. Determination of the total luminosity is important because, in 
addition to the kinematics, it provides a measurement of the total dark matter-to-luminous matter 
ratio in the system. The stellar density in ultra-faint satellites is low enough that it is not possible 



to obtain photometry from integrated unresolved starlight (Martin et al. 2008 Muhoz et al 



2010). Further, because of the low numbers of stars, removal of interloping stars is particularly 



important, because a single bright star not belonging to the satellite may significantly inflate the 
total luminosity (recall that the faintest of the ultra-faint satellites are about as luminous as a single 
red giant star). With the interloper population removed, photometric density profiles are derived 
assuming that the measured stars are an unbiased sampling of an underlying density profile. It 
is then assumed that the stars in the object are comprised of a single population, in this case an 
old, > 8 Gyr, and metal-poor, [Fe/H] < —2, population. If the sample of stars in the system is 
complete brighter than a given magnitude, from the stellar population model the total number of 
stars at fainter magnitude can be determined, and from this the total integrated luminosity of the 
system. 

Similar to the classical dSphs, the ultra-faint satellites are not spherical systems. From the 
isodensity contours, the ellipticities of the systems, again measured in terms of the ratio of the 
minor to the major axis, are in the range ~ 0.3 — 0.7 (Martin et al. 2008), and as a mean population 
are more elongated than the classical dSphs. Given the small samples of stars that belong to these 
systems, and the uncertainty in the precise properties of the dark matter halos that they inhibit, it 
is difficult to determine if they are elongated by tidal effects, in particular because their direction 
of elongation is random with respect to the Galactic center (Muhoz et al. 2010). 

Before applying the theoretical models discussed above to the ultra-faint satellites, and an- 
swering the important question of their dark-to-luminous mass ratio, it will be useful to obtain an 
estimate of the velocity dispersion for a system with the total luminosity of an ultra-faint satellite 
that is bound and in equilibrium. For this example consider the system Segue 1, which is the least- 
luminous of the ultra-faint satellites. For a simple estimate, we can appeal to the virial theorem, 
which states that the total mass, M v i r , of a system with a half-light radius of is related to the 
velocity dispersion as a 2 ~ IMr^jG (note that here we have approximated the total extent of the 
system as the half-light radius). If we assume a mass-to-light ratio for the stars as = 3, as 

is appropriate for a stellar population as observed in the ultra- faint satellites (Maraston, 2005), 
from the total luminosity the velocity dispersion of this system would be approximately 0.5 km/s. 
For the brightest of the ultra- faint satellites, the corresponding velocity dispersion is < 1 km/s. 
A measured velocity dispersion that is significantly larger than this reveals that a dark matter 
component dominates the dynamics of the system. 

The first measurements of stellar line-of-sight velocities in ultra-faint satellites were obtained 

[2009] ). All of these 



a few years ago (Martin et al. 2007 Simon and Geha, 2007 Geha et al. 



measurements determined that the velocity dispersions are in the range of 4-5 km/s; the velocity 
dispersion for each satellite, averaged over the entire radius of the satellite, is shown in Table |3| 
From the estimates above, these velocity dispersions indicated that the masses of these systems are 
significantly greater than what is contributed by the luminous population of stars. Interpreting 
these results at face value, the ultra-faint satellites contain dark matter mass-to-light ratios in the 
range of hundreds to thousands, making them the most dark matter dominated objects known. 
Combined with the measurements of the internal spreads in metallicities, these systems are clearly 
more akin to galaxies than they are to the less extended globular cluster systems that do not contain 
appreciable dark matter. 
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Figure 5 Histogram of line-of-sight velocities of Segue 1 (left) and Coma (right), compared to a 
model with an NFW density profile that provides a fit to the data. 



Figure [5] shows an example velocity distribution for the ultra- faint satellites Segue 1 and Coma 
Berenices. Also shown are fits to the velocity distributions assuming cusped NFW models for the 
dark matter distribution and using Eddington's formula in Equation [l3j Given the relatively small 
data samples, it is not surprising that good fits can be obtained to the data for an input potential. 
Segue 1 has a luminosity of a few hundred solar luminosities and is at a distance of 23 kpc, and 
Coma Berenices has a luminosity of a few thousand solar luminosities and is at a distance of 
approximately 44 kpc. Geha et al. (2009) originally derived in detail the mass distribution of Segue 
1 from a sample of 24 line-of-sight velocities, concluding that it was the most dark matter dominated 
system known. Even though the sample of stars was small, amongst ultra-faint satellites Segue 1 
is a unique system given that it has a systemic velocity of ~ 205 km/s, which nicely separates it 
from the velocity distribution of Milky Way foreground stars. Folllow up studies of Segue 1 with 



three times the sample of stars has shown that this result holds (Simon et al. 2011), and that it is 



robust to systematic effects such as binary stars (Martinez et al 



2011). 



Though as a population there is now strong evidence that the ultra-faint satellites are dark 
matter-dominated, there are outlier systems with unique observational properties, so it is still 
necessary to examine each system on an individual basis. An example of such an intriguing system 
is Willman 1, which presents a unique set of challenges from a few perspectives. Though Willman 
1 was the first ultra-faint satellite discovered, and was hinted to be a dark matter dominated object 
from early studies, the newest results present a less clear picture (Willman et al. 2011). With a 
systemic velocity of -11 km/s, there is significant contamination from Milky Way foreground stars 
in the direction of Willman 1, so determining a clean membership sample is difficult. While this 
object does appear to have properties more in common with that of a galaxy than a globular cluster, 
the velocity distribution shows a double peaked structure that is difficult to interpret. 

Given the relatively small kinematic samples from the ultra- faint satellites, it is not possible to 
determine whether they have cored or cusped central dark matter distributions. However, similar 
to the classical satellites, the uncertainty on the mass profiles is minimized at about the half- 
light radius. The best constrained integrated half-light radius mass for an ultra-faint satellite is 
determined to approximately 20% (Canes Venatici I) (Wolf et al. 2010). This can be compared 
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Figure 6 Mass within 300 parsecs versus luminosity for Milky Way satellite galaxies. From Strigari 



et al. (2008). 



to the best constrained classical dSph, Fornax, which has a half-light radius mass determined to 
approximately 5%. 

The above discussion details how the dark matter distributions in the classical and ultra-faint 
satellites are determined; it is appropriate to put these determinations into a broader context. As 
outlined above, it is not straightforward to compare the classical and ultra-faint satellite popula- 
tions because of their wide range of observed properties, in particular their approximate order of 
magnitude range in half-light radii. Since the total dark matter halo masses are unknown from the 



kinematic data samples, a better observational proxy is the mass within a fixed physical radius ( Stri- 



gari et al. 2007b, 2008); a convenient radius may be chosen as the approximate average half-light 



radius for the satellite population. Phrased in this context, Strigari et al. (2008) highlighted how 



there exists a common mass scale for all of the Milky Way satellites of IO* 5 Mq within 100 pc, and 
10 7 M within 300 pc (Fi gure ^ . The intriguing aspect of this result is that it highlights how the 
Milky Way satellites share a mean integrated central dark matter density, which is independent of 
the luminosity of the satellite. This is particular striking considering that the luminosities of the 
Milky Way satellites span a range of nearly five orders of magnitude. 

Though there have been subsequent revisions to the velocity dispersions of the ultra-faint satel- 



lites (Aden et al. 2009), the salient point still remains from the Strigari et al. (2008) analysis. 



It is certainly plausible that this scale could be set by dark matter physics that differs from the 



CDM framework (Lee and Lim, 2010), though it can also be viewed as a prediction of CDM-based 



models of galaxy formation (Maccio' et al. , 2009 Li et al. 2009; Stringer et al. 2010; Busha et al 



2010). In addition, in the context of CDM, it has important implications for indirect dark matter 



detection, as will be discussed in more detail in Sections [6] and [7[ In particular, if the dark matter 
annihilates to provide an appreciable gamma-ray flux, the luminosity of a galaxy in gamma-rays is 
independent of its optical luminosity. 

4-4- Core versus cusp issue 

The discussion above has focused on determining integrated masses of dSphs. Given the sig- 
nificant improvement in the kinematic data sets, it is also prudent to now ask if the data are able 
to distinguish between cored or cusped central dark matter density profiles. This question is of 
direct relevance for Galactic dark matter searches, because different dark matter particle models 
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make different predictions for the central density of dark matter halos. As an illustration, it is 
possible to define a phase space density for the dark matter as Q = p/& 5 , where here a is the three- 
dimensional velocity dispersion and p is the density (Hogan and Dalcanton, 2000). This quantity 
can only decrease from its primordial value as a result of collisional shocks or coarse graining. For 
a neutrino-like thermal relic particle, the phase density is given in terms of the particle mass, m, 
as 

m \ 2 M /pc" 3 
lkeV/ (kms 



-1A3- 



(44) 



The key assumption in the derivation of this Equation is that the particle freezes out while it is 
still relativistic. It is also possible to obtain phase space densities of this order with non-thermal 
particles produced from decays (Kaplinghat, 2005 Strigari et al. , 2007c). Cold dark matter particles 
have Q values that are orders of magnitude larger, because their velocity dispersion is negligible at 
freeze-out. 

In order to probe the phase-space properties of the dark matter, the ideal systems to study are 
those that are the most dark matter dominated, so that central dark matter density is set by phase 
space properties, and not by interactions with baryons. The most dark matter-dominated systems 
known, whose core sizes are measurable given the scale of modern observations, are the dSphs and 



the low surface brightness galaxies (Simon et al. 2005 Oh et al. 2011). It is certainly true that a 



measured deviation from a cuspy profile in these galaxies would be the most significant clue from 
astrophysics that dark matter is different than is assumed in the LCDM model. 

This motivates a more detailed examination of the central density profiles of dSphs. Before 
examining them from the perspective of the line-of-sight velocity data, it will first be worthwhile 
to review related observations that test for the existence of a dark matter core or cusp. Indeed, 
several authors have suggested that observations not directly related to the analysis of the line-of- 
sight velocity data indicate that the dSphs may favor cored dark matter density profiles. 

As a first example, in Fornax dSph there are five surviving globular clusters, making it, along 



with Sagittarius, the only dSphs with surviving globular cluster systems. Goerdt et al. (2006) 



determined that, based on dynamical friction arguments, the globular clusters would not have 
survived to the present day if the dark matter density profile rises as 1 jr towards the center of the 
galaxy. These authors surmise that the only explanation for their existence is a dark matter core in 
Fornax, of the scale < 1 kpc. Though this result is intriguing, it is based on the projected positions 
of the globular clusters; three-dimensional measurements of the globular cluster distances from the 
center of Fornax are required for more precise calculations. 

There are also arguments in favor of a core in the Ursa Minor dSph based upon observation 
of kinematically-cold stellar substructure (Kleyna et al. 2003, Gilmore et al. 2007). A stellar 
substructure that is this old, > 10 Gyr, is not expected to have survived in a galaxy of this size 
if the dark matter distribution is significantly cuspy, though may survive in a system with a dark 
matter core < 0.85 kpc. Pace et al. (2012) corroborate the presence of a kinematic substructure in 
Ursa Minor, though these latter authors find a different magnitude for the velocity dispersion of the 
substructure, and speculate that this substructure may have its own dark matter halo associated 
with it. 

In order to test these observations, it is best to turn to the line-of-sight velocity dispersion 
data. The likelihood analysis discussed above provides a means for determining whether a cored or 
a cusped dark matter distribution provides a better description of the data. However, significant 
degeneracies preclude robust conclusions. The most significant degeneracies are between the central 
slope of the dark matter density, the central slope of the stellar density, and the velocity anisotropy. 
Before addressing the degeneracies between these quantities in a more rigorous numerical manner, 
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it is first worthwhile to set up an analytic description, for which we again appeal to the spherical 
jeans equation. 

For isotropic models, the spherical jeans equation reduces to a form of hydrostatic equilibrium, 



dr 



d<5> 

' dr ' 



(45) 



where the quantity Ps&r plays the role of the pressure. As a further simplification take the velocity 
dispersion to be constant, which is a good approximation to the dSph data, and a r = ai os . Taking 



the potential to be dominated by the dark matter, combining Equation [45] with poisson's equation 
then gives 



dr 
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pr 



(46) 



los 



where p represents the density of the dark matter, and 7 S = — d(lnp s )/dlnr is the log-slope of the 
stellar density profile. In the limit that r — ¥ 0, assuming that 7 S = Equation [46] implies that the 
central dark matter density must also be cored. If j s > 0, then the central density of the dark 
matter must be cusped as p oc r~ a , with a < 2. In other words, assuming isotropy, if the stellar 
density is cusped then the dark matter density is cusped, while if the stellar density is cored then 



the dark matter density is cored (Evans et al. 2009). 



With this background we can now go into more detailed modeling with the spherical jeans 
equation. In order to implement the jeans equation we need an assumption for the dark matter 
density profile; take this to be of the form Equation 31 The log-slope of this profile as a function 
of radius is 

dlogp ( 
= a — [a — c) - 
\ r o, 



dlogr 



1 + 



To 
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Note that this is a five parameter model for the density profile. Also including a model for the 
velocity anisotropy, /3, implies at least six parameters must be input when fitting to the line-of-sight 
velocity data. Despite the quality of the data sets described above, theoretical analysis shows that 
the degeneracies amongst parameters are so significant that it is not possible to distinguish between 
central cores and cusps. This can be stated on more statistical grounds for hypothetical data sets 



that are even larger than the current sets (Strigari et al. 2007a). The addition of stellar proper 



motions hold promise for breaking these degeneracies (Wilkinson et al. 2002; Strigari et al. 2007a), 



though the required accuracy of approximately a few km/s per star will likely be challenging for 
modern observatories. 

The core versus cusp issue can also be studied via distribution function and Schwarzschild based 
methods outlined above. As an example, Strigari et al. (2010) study five classical dSphs (Fornax, 
Sculptor, Carina, Sextans, and Leo I) with high-quality kinematic and photometry data. They 
assume an isotropic stellar velocity dispersion, allow for cusps in the stellar density p s , and utilize 
cusped central dark matter density models that are predicted in LCDM simulations. The line-of- 
sight velocity distribution is determined from the isotropic Eddington formula. Though the models 
presented are in general non-analytic, numerical solutions find that the cuspy LCDM-based profiles 
may fit the full line-of-sight velocity distributions, even under the assumptions of isotropic stellar 



orbits. This result generalizes the simple analytic formula presented in Equation 46 for line-of-sight 
velocity profiles that are not strictly constant. 



More recently, Schwarzschild modeling has been applied to Fornax ( Jardel and Gebhardt , 2012 ), 



Draco (Jardel et al. 2012), and Sculptor. For Sculptor, 



47 and confirm the results from moment- 



zation of the dark matter profile presented in Equation 
based methods that the kinematic data set from Sculptor is unable to constrain the central slope 



Breddels et al. (2012) use the parameter!- 
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a. Interestingly these authors do find that the velocity anisotropy parameter (3 transitions from 
positive, implying radially-biased orbits, to negative, implying circularly-biased orbits, outside of 
approximately the half-light radius. For Draco, Jardel and Gebhardt (2012) find that the data are 
consistent with a cusped central density profile with a = 1.0 ± 0.2 for r > 20 pc. 

The results above seemingly indicate that even with high quality modern data sets, the line-of- 
sight velocity data is unable to distinguish between cored and cusped dark matter central density 
profiles, possibly a fundamental limitation of the information that is extractable from the data. 
However, more recently a new approach has come to light that looks at the data from a different 
perspective. It does not revolve around larger data sets, but rather around a reinterpretation of 
the data sets that already exist. The idea is straightforward, and goes back to the discussion of 
the previous section on multiple stellar populations. Because stellar populations can be separated 
based on properties such as their metallicity, velocity, and spatial distribution, they may be viewed 
as distinct kinematic data sets. This is of practical interest for the dSph data because a given data 



set can be reinterpreted as multiple data sets that probe the same underlying potential ( Walker and 



Penarrubia 3 2011). Depending on the nature of the data sets, degeneracies may be broken between 



the stellar density, velocity anisotropy, and the dark matter density, possibly pinning down more 
precisely the central dark matter distribution. 

The Sculptor dSph has received the most attention, because it is observed to have two distinct 
populations of stars: a metal rich component with [Fe/H] > —1.5, and a metal poor component 
with [Fe/H] < —1.7. The half-light radius for metal rich component is ~ 230 kpc, and the half-light 
radius for the metal-poor population is 350 kpc (Battaglia et al. , 2008). From jeans-based modeling 
with the assumption of cored Plummer profiles in each case, Walker and Penarrubia (2011 ) indicate 
that a cored model is favored for Sculptor, ruling out a cusped NFW model at > 96% c.l. Because 
of the implications of this result, it is imperative to scrutinize the assumptions in more detail, in 
particular the assumptions of the slope of the central stellar density profile, which was shown above 
to be important in determining the central density slope (Strigari et al. 2010). 

A simplified analytical model has recently been presented to expand upon the results obtained 



m 



Walker and Penarrubia (2011), that derives from the projected virial theorem (Agnello and 



Evans, 2012). In general, it is difficult to implement the projected virial theorem to dSph data sets, 



because it involves an integral over the line-of-sight velocity dispersion, weighted by the spatial 
distribution of the stars. Performing this integration requires a fit to the line-of-sight velocity 
dispersion, so it is a weighted average of the second moment of the distribution. It is clear to see 
how this method loses information relative to using the information that is directly contained in 
the line-of-sight velocity distribution. Because of these issues there are still significant degeneracies 
that preclude the measurement of the mass distribution from a single population, however, for the 
case of multiple populations there is a twist on the argument. 

For a spherically-symmetric system, the projected virial theorem provides a simple relation 
between the line-of-sight kinetic energy tensor, K[ os , and the potential projected along the line-of- 
sight, Wi os , 2K[ os + Wios = 0. In terms of the observable properties of a galaxy, these quantities 
are given by, 
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The key idea brought to light in Agnello and Evans (2012) is that if there are two populations of 
stars in a dSph, each with a distinct surface brightness profile and velocity dispersion then the ratio 
of the virial quantities may be taken. For two populations, labeled 1 and 2, the virial theorem then 
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implies that if/ os ,i/^os,2 = Wios,l/Wios,2- So from the surface brightness profile and the velocity 
dispersion for each population, a mass distribution that generates the underlying potential can be 
assumed, and it is then possible to determine if this mass distribution satisfies the ratio of the 
virial quantities. Though it is a more simplified theoretical argument than the jeans-based analysis 
presented above, perhaps the most intriguing aspect of this method is that it depends only on the 
projected light profile, and is not subject to the degeneracies incurred from deprojection of the light 
profile. 

Though theoretically-interesting, the projected virial theorem is probably still not the final 
word on the issue of multiple stellar populations in dSphs. Systematic issues that pertain to the 
use of the projected virial theorem include testing the assumption of spherical symmetry, testing the 
robustness to the variation in the projected light profile, and understanding if significant information 



is lost when performing the integrals in Equations 48 and 49 

Though there are some intriguing hints of dark matter cores in dSphs, a more definitive answer 
will likely have to wait until proper motions for individual stars are obtained, as well as higher 
quality photometric data sets in the central regions. Both of these observations will help in breaking 
the degeneracy between the velocity anisotropy and the central dark matter slope. 

4-5. Streams and disrupted satellites 

To this point this section has focused on understanding bound, self-gravitating populations of 
stars and dark matter around the Milky Way and the extraction of their dark matter distribu- 
tions. Though we now have measurements of over two dozen dark matter-dominated satellites, 
and approximately six times more globular clusters, it is very likely that many of these objects are 
transients on cosmological timescales. Interactions with the Galactic potential may induce many 
of them to dissipate their stars and dark matter into tidal streams and debris that will eventu- 
ally phase mix into the Galaxy. Understanding the dynamics of these debris streams is not only 
important as a point of comparison to the predicted distribution of disrupted satellites in LCDM 
simulations (Section [5]), it is of direct relevance for direct dark matter searches, for which event 
rates are significantly affected by dark matter streams passing through the Solar neighborhood 
(Section [8}. 

There is now significant evidence that tidal disruption is a common phenomenon that is an 
integral part of the formation of the Galaxy. Several streams have now been identified as distinct 
features in phase space. The most spectacular example is the Sagittarius dSph; this object was 
previously utilized in Section [3] as a probe of the shape of the Milky Way dark matter halo. The 
major axis of stars associated with Sagittarius extends well over 10°, and a nearly constant velocity 



dispersion of ~ 10 km/s is now measured in the core region of ~ 2 kpc (Frinchaboy et al. 2012). 
The core of Sagittarius has a luminosity that is similar to Fornax, and also like Fornax, Sagittarius 
has several globular clusters associated with it. Modeling indicates that the dark matter mass of 



Sagittarius before it fell into the Milky Way halo was in the range 10 — 10 M (Law et al. , 2005 ) 



Another intriguing stream in the Milky Way is the Magellanic stream, which is a neutral hy- 



drogen structure that appears to trail behind the Magellanic Clouds (Mathewson et al. 1974). 
Previous modeling has associated this stream with the Magellanic clouds, relying on the hypothesis 
that the Magellanic clouds have undergone multiple close passages with the Galaxy, and the stream 
is a result of either ram pressure stripping or tidal stripping. However, recent proper motions have 
questioned this multiple passage model, showing instead that the clouds are most likely on their 



first passage through the Galaxy (Kallivayalil et al. , 2006). More recent modeling has shown that 



it is possible to connect the Magellanic stream to the clouds under the first passage scenario (Besla 



et al., 2010 Diaz and Bekki, 2012) 
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Table 4. Prominent stellar streams in the Milky Way halo that may be associated with 

tidally-disrupted dwarf galaxies. 



Name 



Distance [kpc] Reference 



Magellanic Stream 
Sagittarius 
Orphan Stream 
Styx 

Bootes III 
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Mathewson et al. 



Ibata et al. 



(1994) 



Belokurov et al. 
Grillmair| ( |2UM[ 
Grillmair| 



(1974) 



(2007a) 



Carlin et al. 



(2009) 



In addition to the Sagittarius and the Magellanic streams, several other streams have recently 
been identified in the SDSS ( [Belokurov et al.[ |2007a[ |Duffau et akj |2006[ |Grillmair[ |2009[ ). These 
streams may have originated either from a dwarf galaxy or a globular cluster. Table [4] lists promi- 
nent streams in the Milky Way halo that may be plausibly associated with a tidally-disrupted 
dwarf galaxy (several other known streams exist in the Milky Way, but likely originated from the 
disruption of a globular cluster). The detection of a multiple streams, either from disrupted dwarf 
galaxies or globular clusters, lends more support to the theory that the stellar halo of the Galaxy 
was formed via the chaotic accretion model (Searle and Zinn, 1978), as opposed to the theory of 



formation via monolithic collapse (Eggen et al. 1962) 



Probably the simplest means to determine whether a stream originated from a dark matter- 
dominated dwarf galaxy or from a globular cluster or star cluster is through the width of the stream. 
For example, Grillmair (2009) identified four streams in the SDSS data ranging in distances from 3 
to 15 kpc (Acheron, Cocytos, Lethe, and Styx), and indicates that three of them have widths < 100 
pc, while the fourth, Styx, has a width of approximately 2.6 kpc. This conclusion is strengthened 
by measurements of the velocity dispersion of ~ 14 km/s of the Bootes III object, which is a clump 
of stars that appears to be in the initial phases of tidal disruption and may be the progenitor of 
the Styx stream (Carlin et al. 2009). An additional diagnostic of the origin of a stream comes 



from understanding the metallicity of the constituent stars. Boer et al. (2012) show that dSphs are 
typically deficient in Ni and Na relative to globular clusters, and use this diagnostic to conclude 
that the progenitor of the recently-discovered Aquarius stream is a globular cluster. 

In addition to their role in understanding the population of accreted dwarf galaxies, streams 



are also excellent probes of the population of dark matter substructures in the Galaxy. Carlberg 



(2012) showed that gaps in streams are created by dark matter substructure, and that the high 
density of gaps in narrow streams may be indicative of a population of > 10 5 substructures greater 
than approximately 10 6 M . If the orbit of a stream can be reconstructed, it is also possible to 
constrain the shape of the Milky Way potential that gives rise to that stream. This was evident 
in the discussion of the Sagittarius stream in Section [3j However, due to the typically sparse data 
associated with the streams, it is difficult to reconstruct the orbital trajectory and obtain unique 
fits to the data for a given potential model (Eyre and Binney, 2009). 

Finally, a discussion of the population of Galactic satellites and stellar streams cannot be 
complete without mentioning the unique spatial distribution of the satellites. It was noted many 
years ago that the observed Milky Way satellites are distributed to a good approximation along 
an inclined plane relative to the Milky Way disk (Lynden-Bell, 1976). As discussed in the next 
section, this may be a reflection of the fact that dark matter satellites accrete in an anisotropic 
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manner along filaments into the Milky Way (Libeskind et al. 2005, Zentner et al. 2005). It has 



also been hypothesized that the observed population of satellites are so-called tidal dwarf galaxies, 
which form when a gas- rich galaxy merges into a larger host halo, and fragments to produce distinct 
populations of stars (Metz and Kroupa, 2007). 

The spatial distribution of satellites, as well as the other dark matter properties of Milky Way 
satellites and streams that were discussed in this section, comprise a series of compelling issues 
that must be understood in a larger cosmological context. This now leads into the subject of the 
next section, which lays out the predictions of the structure, and substructure, in the Milky Way 
within a cosmological context. 



47 



5. Simulations of Galactic halos: Predictions and observational connections 



The observations discussed in the previous sections have provided a wealth of information on 
the structure of the Galaxy and its satellites. They provide measurements of the local dark matter 
density, the mass of the Galactic dark matter halo, and the masses of the satellite galaxies. These 
quantities will be important for interpreting results from direct and indirect dark matter searches 
(Sections [6] and [8]). It is important to recognize, however, that there is still some limitation to 
the scope of these observations, and that dark matter searches require information not directly 
measured from modern astronomical data. Two particularly important examples are the dark 
matter velocity distribution, and the properties of dark matter substructure. In order to model 
properties such as these that are not directly determined from modern observations, we must appeal 
to theory, which now mainly revolves around the use of numerical simulations. The primary goal 
of this section is to advance the theoretical description of the dark matter properties in the Galaxy, 
using both analytic models and N-body simulations, focusing on those that are of most practical 
use for dark matter searches. 

The first analytic models for the formation of dark matter halos were based on the theory of 



self-similar accretion of material around a seed perturbation ( Gunn and Gott , 1972 ; Fillmore and 



Goldreich, 1984; Bertschinger 1985). Though a great deal of theoretical insight has been gained 
from this theory, it is becoming clear that this picture of self-similar accretion must be considered 
together with a more violent and chaotic picture for the formation of dark matter halos that involves 
hierarchical merging of smaller mass halos. In the hierarchical merging scenario, dark matter halos 
of a given mass can have significantly different formation histories; for example some halos may 
experience a violent late merger, while others may experience a quiescent merger history. 

In recent years, it has become widely appreciated that the "smooth" component of a dark matter 
halo is only part of its phase space structure. One of the more profound theoretical insights in dark 
matter over the past couple of decades is that halos formed in a LCDM universe are not smooth, but 
rather contain an abundance of substructure. This substructure contains a wealth of information on 
the nature of particle dark matter. Indeed, if dark matter has an annihilation or an elastic scattering 
cross section that is too low to be detectable by indirect or direct methods, probably the best 
method for studying the nature of particle dark matter will be through the properties of satellites 
and substructure. As a result there is clear motivation for a detailed theoretical understanding 
of the mapping between a particle dark matter model and its substructure properties, and how 
different models can be studied observationally. 

This section begins with a description of numerical simulations that are used to study the 
structure of dark matter halos. It then goes into more detail on the output of the simulations, in 
particular the important quantities such as the mass and the velocity distributions. The theoretical 
predictions for dark matter substructure in the Galaxy are then discussed, and the connection is 
made to the observations of satellite galaxies discussed in Section [4} 

5.1. Simulations of Milky Way-like galaxies 

Since this section will frequently quote the results from numerical simulations, it is beneficial 
to have a general understanding of the most important properties of the simulations. Standard 
N-body (i.e. dark matter-only) simulations assume that dark matter is the dominant component 
of mass in the universe, and that the particles responsible for setting the formation of structure 
interact only gravitationally, so their direct encounters can be ignored on cosmological scales. In 
order to run the computations in a reasonable amount of time on modern computers, a single dark 
matter particle is made to represent a much larger number of smaller mass particles. 
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In a simulation the unperturbed universe is represented by a uniform grid of particles. From seed 
linear density perturbations on this smooth model, structure in the universe forms as overdense 
regions grow and eventually break away from the expansion of the universe to form self-bound 
gravitating dark matter halos. The number and the structure of the halos depends sensitivity on 
both the initial conditions for the density perturbations and the nature of the dark matter. In 
the cosmological volume within which modern numerical simulations are run, the mass function of 
dark matter halos is well-determined in the range of approximately 10 8 Mq to scales as massive as 
clusters of galaxies, approximately 10 14 M Q . 

Halos at the mass scale of the Milky Way, ~ 10 12 Mq, are identified in a larger volume N-body 
cosmological simulation at redshift z = according on their mass and a density criteria on the 
matter distribution in the local vicinity of the halo. Once the halo is identified, a Lagrangian 
region within which the halo formed is determined. Within this region, a larger number of lower 
mass particles are generated according to a new set of initial conditions. Regions of larger and larger 
distance outside this Lagrangian volume are populated with more and more massive particles. Since 
the sampling is much higher in the Lagrangian region than in the initial larger scale simulation, 
not all of the particles within this region wind up in the final halo, and a negligible fraction of more 



massive particles in the outer regions end up in the final halo (Springel et al. 2008a). 



The above procedure for simulating the formation and the subsequent identification of Milky 
Way mass dark matter halos in an expanding universe works well for particles that travel a small 
distance compared to the grid size at the redshift when the simulation is started. This is a good 
approximation for the cold and warm dark matter models considered here. However, in warm dark 
matter models, the sharp cut-off in the power spectrum gives rise to artificial halos that form along 



filaments on scales that are characteristic of the inter-particle separation (Lovell et al. 2012). This 



turns out to be an important effect in the regime of interest for both observations and theory. 
5.2. Mass distribution 

The mass profiles of dark matter halos in N-body simulations are of great interest for many 
aspects of this article. In combination with the observational results discussed in Section [3j the 
results from N-body simulations are used for predicting signals in indirect dark matter experiments 
(Section [7]). Because the formation of halos and the resulting mass profiles are directly related to 
the properties of dark matter, before delving into the results from the modern simulations, it is 
worthwhile to take a step back and understand different aspects of theoretical predictions for mass 
profiles of halos in an expanding universe. 

The first predictions for the mass profiles of galaxies that form in an expanding universe were 
based on the self-similar accretion model. The self-similar model begins by defining a spherical 
region, and assumes a perturbation in this region of the form 5M/M ~ M~ e , where the case 
e = 1 corresponds to accretion about a point mass (Gunn and Gott, 1972). Assuming spherical 



symmetry and entirely radial orbits, self-similar models predict single power-law models for the 



density profile, of the form p oc r 2 for e < 2/3, and p oc r 9e /( 1 + 3e ) for e > 2/3 (Fillmore and 



Goldreich, 1984). These results predict that density profiles depend on the halo mass, and the 
density profiles are always steeper than r~ 2 . Because the density at a given position is set by a 
sum over mass shells that contribute to the density, the interior of the density profile is smooth as 
a result of the summation of a large number of mass shells, whereas in the outer region of the halo, 
only a small number of shells contribute to the density profile. The power law shape of the density 
profile breaks at the transition to the infall regime, where the density profile is characterized by 



sharp caustic features (Bertschinger, 1985). 



Self-similar halo models have now been extended to include non-radial motions and deviations 
from spherical symmetry, and solutions for the density profile have been obtained. The addition 
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of angular momentum increases the range of slopes allowed, predicting a pure power law that 
scales as r -1,3 for Milky Way-mass halos (Nusser, 2001). When torques are introduced into the 



spherical infall model, the central slope is mass dependent, with larger mass halos being described 



by steeper density profiles (Zukin and Bertschinger 2010). It is also interesting to note that even in 



the simplest case of pure radial orbits, numerical simulations with spherical initial conditions show 



that the radial orbit instability turns the resulting halo into a non-spherical structure (Vogelsberger 



et al. 2011) 



Understanding self-similar halo models is not just important because it provides physical in- 
tuition for halo formation, it also brings up the important question of the dependence of density 
profile on halo mass. A related question is the scatter in the density profile at a fixed halo mass. 
If density profiles are found to not depend on mass, or if there is significant scatter in the density 
profile at a fixed mass, it would provide direct evidence that the physics of halo formation is not 
entirely captured in the self-similar model. 

In numerical simulations, the formation of halos does not occur through a smooth and nearly- 
spherical accretion, but rather through a combination of smooth particle accretion and the merging 
of lower mass halos. It is still a matter of debate what fraction of the mass of a halo comes 
from mergers and what fraction comes from smooth accretion; some studies find that the merger 



component constitutes nearly the entire mass of the halo (Angulo and White, 2010), while others 



find that approximately 50% is due to mergers (Stewart et al. 2008. Genel et al. 2010) 



Whatever the origin of the mass is that comprises dark matter halos, it is clear that halos 
that form via hierarchical clustering are not spherical. Analysis of the shapes of halos over a 



mass range of approximately 10 10 — 10 i4 M indicate that halos become less spherical as the mass 
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increases (Allgood et al. 2006; Schneider et al. 2012). Halos are found to be more preferentially 



prolate, and the major axis is aligned with the large scale structure distribution, indicating an 
intriguing interplay between the properties of the halo and the local environment. 

In order to best understand whether the density profiles of halos are universal, it is most 
convenient to work in the context of spherically- averaged profiles. Initial studies of the mass 
distribution of dark matter halos in cosmological simulations indicated that the central slope is 
shallower than is predicted by the secondary infall model, with the central slope reaching a constant 
asymptotic inner value p ~ r _1 (Dubinski and Carlberg, 1991 ). Further, it was found that the outer 
regions of the density profile fall-off more steeply, as r _a , and that this double power law model, 
the Navarro-Frenk- White (NFW) profile, is nearly universal in halos over a wide range of mass 



scales (Navarro et al. 1996| ) 



Though these results put forth the profound idea that density profiles are the same for any 
mass halo, they still suffered from insufficient resolution of the central regions. The recent sample 
of Aquarius simulations, which resolve down to less than one percent of the halo virial radius for 
approximately 10 12 M Q halos, now find that this asymptotic slope is not yet achieved, and that the 



profile becomes progressively more shallow towards the center of the halo (Navarro et al. 2010) 
These recent results favor the Einasto profile, 



p e (r)=p- 2 exp[-((r/r_ 2 r-l))], 



(50) 



where r_2 is the radius at which a log slope of 
radius. For Milky Way-mass halos, a 



-2 is achieved and p_2 is the density at that 



for subhalos (Springel et al. 2008a) 



= 0.17 ( jNavarro et al. 2010). This value of a also works well 
though recent results suggest that a modification may work 
better (see discussion of substructure below). In addition to the fact that the density profile does 
not achieve a pure power law, there is a small mass trend that is seen in the density profiles as 
a function of halo mass. However, neither of these trends are as significant as is predicted in the 
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self-similar halo formation model as described above. These simulations also show that different 
halos at a fixed mass cannot be scaled to look identical. 

In addition to the density profiles, simulations predict additional important trends for the 
structure of dark matter halos. It is found that lower mass halos are more centrally-concentrated 
than halos of larger mass, where the concentration can be defined as the ratio of the radius at 
which the density reaches 200 times the critical density, defined as r2oo> to the characteristic radius 
of the halo. For the Einasto profile, the concentration is c = 7"2oo/ r -2> wnue f° r t ne NFW profile 
the concentration is c = ^oo/^O; where ro is defined in Equation 16 For the updated values of 



cosmological parameters and for halo masses above approximately 10^ M , the concentration- mass 



relation for the Einasto profile is (Duffy et al. 2008) 



8.82 



l + z 



,0.87 



M„ 



2 x lO^/i-iMQ 



-0.106 



(51) 



It is important to ask whether the halo profile in Equation 50 and its corresponding properties 



are specific to the LCDM cosmological model. Maccio' et al. (2012) have recently examined the 



density profiles in universes in which a component of the dark matter is warm. Due to the cut 
off in the matter power spectrum in warm dark matter models, halos typically form later, so their 
concentrations are lower than in pure cold dark matter universes. Though the concentrations are 
lower, the NFW model still is a good description of the density profile, which points to the intriguing 
idea that the NFW-like profile results from a general dynamical processes in halo formation. 

All of the results discussed to this point concern pure dark matter-only simulations, and do 
not include the effect of baryonic physics. Inclusion of baryons leads to a great deal of challenges 
from the perspective of numerical simulations, so understanding halo structure in more realistic 
universes will likely remain elusive for some time. Higher resolution simulations of single halos, in 
combination with analytic models, have recently indicated that the central densities of dark matter 
halos become less steep than those found in pure N-body simulations because the baryons induce 
repeated epochs of feedback due to star formation activity. Cusps are turned into cores because 
energy is repeatedly input into the collisionless particles due to fluctuations in the gravitational 



potential (Pontzen and Governato, 2012). This model provides an appealing explanation for low 



surface brightness galaxies, for which there is now good indication that the central densities are 



flatter than the Einasto model (Oh et al. 2011). 



5.3. Velocity distribution 

The density profiles of dark matter halos are only a part of the phase space information that is 
of interest for dark matter searches. A second important quantity is the dark matter velocity distri- 
bution function (VDF) at the Solar radius. While there has been a significant amount of literature 
devoted to understanding the universality of halo mass profiles, by comparison relatively little work 
has been dedicated to understanding the global properties of the VDF in N-body simulations. The 
VDF is of particular importance for direct dark matter searches (Sections [8] and [9]), which are most 
sensitive to the particles in the high velocity tail of the distribution due to the kinematics of the 
collision between the dark matter and nucleus. Indirect searches for high energy neutrinos from 
dark matter particles that get captured in the Sun (Sections [6] and [7]) are sensitive to the lowest 
velocity component of the VDF. 

Following the discussion above on halo density profiles, it is first informative to gain some 
theoretical intuition to the nature of the VDF. For this we can compare to the observations of 
collisionless stars in elliptical galaxies. For stars in elliptical galaxies, it is straightforward to show 



that the two-body relaxation scale is much longer than the Hubble time (Binney and Tremaine 
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2008). However, since their density profiles are near isothermal, a physical mechanism beyond 



two-body relaxation must set the dynamics of galaxies, unless for some reason galaxies form in 
a state near equilibrium. The above paradox led to the development of the theory of 'violent 



relaxation' by Lynden-Bell (1967). The main premise of this theory is that particles in a galaxy do 



not respond to interactions with individual objects, but instead to changes in the form of groups 
of objects, which causes the potential of the system to vary rapidly during its formation. The stars 
continually respond to different groups in a statistical manner; for completely random interactions 
the central limit theorem dictates that the velocity distribution approaches a maxwell-boltzmann 
distribution and thus an isothermal density distribution. It is straightforward to see why this logic 
may translate over to dark matter in halos, which form via a combination of mergers and accretion. 

Though the broad arguments above suggest that the VDF should be close to a maxwell-boltzman 
distribution, it is important to recognize that different VDFs are predicted in other theoretical 
models. For example the self-similar model described above predicts that the density profile at each 
location is comprised of distinct flows of particles. These flows manifest themselves in the VDF 
as a series of delta-functions that are normalized by the density of particles in each flow relative 
to the total density at that location. Also the disruption of a dark matter satellite generates a 
VDF that is represented by a single flow of particles. A significant amount of modern theory is 
devoted to understanding the relative contribution between these non-thermalized flows and a more 
thermalized component 

In modern N-body simulations, the VDF is studied in two separate statistical regimes. On the 
one hand, the resimulation of dark matter halos in a full cosmological volume via the methods 
described above provides a large number of particles in a relatively small number of halos. In 
this case, although single halos are well-sampled, it is difficult to obtain a measure of the scatter 
in the VDF between halos of a given mass. On the other hand, the VDF is analyzed in a large 
sample of halos taken from a large cosmological volume. While this later method does give a better 
estimate of the scatter in the VDF between halos of a fixed mass, accurately measuring the VDF 
in the central region of the halo is difficult due to lack of numerical resolution. In both cases some 
statistical information is lost, either in the total number of halos analyzed or in the number of 
particles per halo. 

To understand this matter more concretely, consider the relevant scales of modern simulations. 
The two largest volume cosmological simulations to date that resolve a statistical sample of Milky 
Way- mass halos identify approximately 10 4 halos with mass near 10 12 M (|Boylan-Kolchin et al. 



2009 Klypin et al. 2011 ). Due to the particle resolution of a few times 10 M Q , halos at this mass 
contain 10 3 — 10 4 particles per halo. For comparison, the Aquarius (Springel et al. 2008a) and 



Via Lactea-II (Diemand et al. 2008) simulations of Milky Way- mass halos contain nearly a billion 



particles within the virial radius of their respective halos. 

In order to interpret the VDF from these simulations, it is best to establish a point of compar- 
ison. With this motivation define the maxwell-boltzmann velocity distribution, 



f(v) 



1 



V27TCJ 2 



exp[-i; 2 /2cT 2 ], 



(52) 



where a is the velocity dispersion, which as discussed below in section [8] can be related to the 
circular velocity of the halo. The maxwell-boltzmann velocity distribution is isotropic in that 
it only depends on the magnitude of the velocity, and may be self-consistently derived from an 
isothermal density profile. Though it is convenient and motivated in part by theory, there is no a 
priori reason why we should expect halos to have this exact distribution. 

As shown in Figure [7j the VDF at the Solar radius differs from a pure maxwell-boltzmann 
distribution in the highest resolution dark matter only simulations ( Vogelsberger et al. , 2009 Kuhlen 
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Figure 7 Velocity distribution function in dark matter only simulations. The left panel shows the 



distribution from the dark matter-only Aquarius simulations (Figure reproduced from | Vogelsberger 
et al. ( 2009| )). The right panel shows the stacked velocity distribution function for halos from Milky 
Way mass to cluster mass scales (Figure reproduced from Mao et al. (2012)). In the right panel, 



from left to right, the VDF is shown at increasing radius relative to the scale radius of the halo. The 
solid dashed- line in the left, and the solid line in the right give the best-fitting maxwell-boltzmann 
distributions. 



et al. 2010). In both cases the deviation from the corresponding best fitting maxwell-boltzmann 
model from Equation |52| is similar. First, the distributions have a larger fraction of particles at 
velocities less than its peak velocity. Second, the peak of the distribution is broader than the 
maxwell-boltzmann model. Third, there are fewer particles in the extreme high energy tail of the 
distribution. For the individual halos studied, these results incorporate the scatter in the VDF at 
the Solar radius by measuring it at different locations within an ellipsoidal radius that is determined 
from the halo density distribution. 

The Aquarius and Via Lactea-II simulations also reveal the presence of sharp features that are 
due to substructure. These features are only recognizable when a halo has a significant fraction of 
its mass, > 50%, in an individual substructure. To a good approximation, the VDF is dominated 
by the "smooth" component, which reflects particles that have been both smoothly accreted and 
those that were once part of a subhalo that has been tidally-destroyed. The VDF also reveals 



broader and more extended features that are related to the energy distribution ( Vogelsberger et al. 



2009). These features are not seen in the individual velocity components, which are well-described 



by gaussians of different widths; the triaxiality of the VDF is clearly indicated in Vogelsberger et al 



(2009) and Abel et al. (2011) 



Turning now to the other regime of identifying a larger sample of halos with lower resolu- 



tion, Mao et al. (2012) have recently studied the properties of the VDF in halos from Milky Way 



mass up to the cluster mass scale. From a stacked analysis of the VDF from thousands of halos, 



the general results from the individual high-resolution halos are corroborated. Mao et al. (2012) 
find that the VDF is well described by the form 



f(v) = exp(-\v\/v )(v - v esc ) 



v 

esc I i 



(53) 
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where vq and p are correlated fit parameters that depend on the scale radius of the halo. The 
key aspect of this VDF is that it is broader at the peak and steeper at the tails than the maxwell- 
boltzmann model, with the power-law index p controlling how steep the distribution falls off beyond 
the peak. The right hand panel of Figure [7] shows the VDF at different radii in a large sample of 
halos from cosmological simulations. 

The above methods rely on direct analysis of the output from cosmological simulations. Com- 
plementary methods for studying the VDF combine these results with various analytical methods. 
One such method involves sampling halos from cosmological simulations, merging them with a halo 
that is similar in mass, and evolving the merged system until it is relaxed. Here relaxed means that 
the system is in equilibrium, so that evolving it further in time does not change its phase space 
distribution. This methodology was utilized in Hansen et al. (2006), who started from an initial 



condition of two halos with an NFW profile, and used them to generate a distribution function 



according to the isotropic Eddington model (Equation 13). Hansen et al. (2006) et al. suggested 



that the best-fitting VDF of the merged structure after many dynamical timescales is the Tsallis 
distribution, 

2" 

/ r 

!-(!-<?) 



f(v) 



V 



1 

1-9 



(54) 



where q ~ 1 is the Tsallis index and k is a free parameter. The maxwell-boltzmann distribution is 
recovered in the limit that q — > 1. 

Numerical simulations also show a correlation between j3, velocity anisotropy parameter (Equa- 



tion 12), and the log-slope of the density profile such that /3 increases with radius from the center 



of the halo. This behavior of the anisotropy can be qualitatively understood in terms of the in- 
fall properties of particles in the halo, because infalling particles with radially-biased orbits tend 
to dominate the outer region of the halo. The Ospikov-Merritt model predicts behavior of the 
anisotropy of this form, 



P{r) 



9 i 9 > 

r z + r p 



(55) 



where rp is the scale radius of the velocity anisotropy profile. Though it does qualitatively provide 
the correct behavior for /3, the slope of the rise of j3 is much steeper than what is seen on average 



in simulated dark matter halos. Hansen and Moore (2006) more accurately quantify the relation 



between the anisotropy and the log-slope, 7, as 



/3 = 1- 1.15 1 + 



(56) 



Recent analysis of cluster mass galaxies corroborate this relation, although they find that there is 



a significant amount of scatter (Lemze et al. 2012). 



5.4- Galactic substructure 

The discussion to this point has focused on the smooth component of dark matter halos. One of 
the great insights into dark matter that has been deduced from N-body simulations is that the mass 
distribution within Galactic dark matter halos is in fact not entirely smooth, but that a significant 
fraction of it is bound up in substructure, or subhalos. Modern simulations of Milky Way-mass 
galaxies are now complete in their measurement of subhalos down to a mass of approximately 10 6 
Mq, corresponding to approximately 10 _6 % of the total mass of the halo. For masses > 10 6 Mq, 
the mass function of subhalos is measured to be a power law that scales as dN/dM oc M~ a , with 
a = 1.9. Because 1 < a < 2, subhalos at the low mass end of the mass function dominate the 
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distribution by number, while the subhalos at the high mass end of the mass function dominate the 
total mass in substructure. The fraction of the total mass of the Milky Way halo in substructure 
is in the range 10% - 50% (IDiemand et al.l [2008] ISpringel et al.l |2008a|). 



The most massive subhalos in a typical Milky Way-mass dark matter halo are approximately 
10 10 Mq. Subhalos near this upper end of the mass range are preferentially found in the outer 
regions of the halo, because dynamical friction effects destroy massive subhlaos as they approach 
the central regions of the halo. Less massive subhalos are found to more faithfully track the smooth 



dark matter density profile ( Springel et al. , 2008a ) . The predicted radial distribution of subhalos 



in primary halos is of particular importance for indirect dark matter searches (section [7]) . 

The formation of substructure down to the mass scale resolved in N-body simulations is a unique 
prediction of cold dark matter theory. For standard cold dark matter particles with mass around 
the weak scale and freeze-out temperatures > GeV, the minimum mass subhalo is predicted to be 



somewhere in the range of an Earth mass (Green et al. , 2004 Loeb and Zaldarriaga, 2005 Profumo 



et al. 2006; Martinez et al. 2009). It is not yet clear if subhalos at this low mass scale would 



survive tidal disruption in the Milky Way halo though (see Diemand et al. (2005)). In warm dark 
matter models or its variants, free streaming induces a cut-off in the power spectrum, and there is 
expected to be a corresponding cut off in the mass function of substructure at a significantly larger 
mass scale. 



5.4-1- Connection to Galactic satellites 

The bound subhalos that survive within Milky Way-mass dark matter halos are expected to 
be related to the dSph galaxy population; understanding this mapping provides a unique probe 
of cosmology and the properties of dark matter. The correspondence between subhalos and dSph 
satellite galaxies is by no means straightforward; it has been known for two decades that there are 
orders of magnitude more dark matter subhalos predicted in simulations of Milky Way-mass halos 



than there are observed dSph satellite galaxies (Kauffmann et al. 1993 Klypin et al. 1999 Moore 



et al. 1999 ). This discrepancy, often known as the "missing satellites problem," has received a great 



deal of attention in the literature (for recent reviews of different aspects of this issue see Kravtsov 



(2010) and Bullock (2010)), and it has inspired a wide range of theoretical variations of the LCDM 



model (e.g. Kamionkowski and Liddle (2000); Zentner and Bullock (2003); Hooper et al. (2007)). 

At the time when it first received a great deal of attention in the astronomical community, 
it was unclear what to make of the missing satellites issue. For example, it was well-appreciated 
that the census of Milky Way satellites was incomplete, even for satellites that are as luminous 
as the brightest dSphs (Willman et al. 2002). In addition, kinematic data sets from dSphs were 



not as complete as they are today, so there was a great deal of uncertainty on their dark matter 
masses. Because of these uncertainties, there were well-motivated theoretical solutions to the 
satellite discrepancy. First, and perhaps most simply, it could be that the observed dSphs reside 
in the most massive dark matter subhalos. This would imply a sharp cut off in the formation of 
galaxies below a fixed present day mass scale, which could result from, for example, supernova 



feedback processes ( Stoehr et al. , 2002 ) . It was also hypothesized that the dSphs do not reside in 



the most massive subhalos at present, but rather in the most massive subhalos that formed before 
a redshift near reionization. This solution directly connects the luminous fraction of satellites to 



the properties of the intergalactic medium at the time when the halos formed (Bullock et al. , 2000 



Benson et al. 2002). Combinations of the above solutions were developed, which hypothesized that 



the dSphs reside in subhalos that were the most massive at the time of accretion into the Milky 



Way (Kravtsov et al., 2004). 



Within the past decade, advances in the census of Milky Way satellites from the SDSS and in 
their kinematic data (see Section [4]) have served to better focus the understanding of the missing 
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satellites problem. Assuming that the subhalos in the simulations host the dSphs, data and the- 
oretical modeling indicates that the peak circular velocities for the dSphs are in the range 10-30 
km s _1 (Strigari et al. 2010), confirming earlier estimates based on sparser data samples (Zentner 



and Bullock 



2003 



Kazantzidis et al. , 2004). This now clearly implies that there are dozens of dark 



subhalos that are devoid of stars entirely, or host stellar populations that are too faint to resolve 



in modern surveys ( Boylan-Kolchin et al. 2011). This is explicitly shown in the left hand panel 



of Figure [9j which compares the measured circular velocities of the dSphs to the circular velocities 
of the most massive subhalos in the Aquarius simulation. Given that the satellite population is 



complete down to at least a luminosity of the classical satellites (Walsh et al. , 2009), in the regime 



surveyed by SDSS it is no longer possible to postulate that a large population of undetected satel- 
lites exist at the same luminosity as the classical satellites, so long as they have surface brightnesses 
similar to the classical satellites. 

Now that the abundances of the Milky Way satellites and subhalos are more secure, does the 
new experimental data point to a solution of the missing satellites problem? Combining both the 
theoretical and the observational advances, the possible solutions can be broadly classified into 
three main categories. First, the solution may be that our knowledge of galaxy formation at the 
scales of the faintest known galaxies in the universe is still incomplete. It is certainly true that 
modern simulations, which are unable to appropriately account for baryonic physics, do not include 
all of the relevant physical processes in galaxy formation. Second, the solution may be that our 
observational knowledge is still incomplete, either in the form of our census in the Milky Way of 
ultra-faint satellites, or in our understanding of the distribution of satellites around a typical galaxy 
like the Milky Way. In other words, the Milky Way may be an outlier in terms of its population of 
luminous satellites. Third, as mentioned above, the solution may be that the cosmological model 
differs from the standard LCDM model. This could result in either a reduction of the number of 
subhalos, a reduction in their central dark matter densities, or both. Any one of these three effects, 
or a combination of them, may be a part of the solution. Each of them are now considered in turn. 

5.4-1 -1- Inclusion of baryonic physics. To better address the first hypothesis, there are now numer- 
ical simulations that have examined the effect of baryonic physics on the formation and structure 
of Milky Way-mass dark matter halos and satellite galaxies (Wadepuhl and Springell 2011; Parry 



et al. , 2012 ). The physical processes included in these simulations are supernova feedback, radiative 



cooling, cosmic ray pressure, and energy input from the supermassive black hole at the center of 
the host galaxy. Because of the low density of baryons in the satellites and the scales that the 
simulations resolve, these cosmological simulations do not point to a significant change in the dark 
matter properties of the satellites. This picture has been challenged recently in simulations of 
isolated satellite galaxies, which indicate that baryonic feedback can significantly affect the central 



densities of dark matter halos (Brooks and Zolotov, 2012). The reduced central densities affect the 



population of satellites in two ways. First, it alters the mass distribution of the most massive satel- 
lites, putting them seemingly more in line with the kinematics of the observed satellites. Second, 
because of the reduced central densities of the satellites, it lends them more susceptible to tidal 
stripping, reducing the overall number of satellites. 

5. 4-1-2. Satellites of external galaxies. Now addressing the second hypothesis, a complete assess- 
ment of the satellites discrepancy requires a better understanding and identification of faint satellites 
around the Milky Way and external galaxies. Observations of the Milky Way satellite population 
were covered in section [4] and will be discussed again in section [TUJ Here two modern methods 
are discussed for examining the satellite population around external galaxies: the first is the direct 
counting of satellites from large scale surveys, and the second utilizes the method of gravitational 
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lensing. 



Direct counting methods 



Recently a great deal of progress has been made on searches for faint satellite galaxies around 
systems that have physical properties similar to the Milky Way. These properties include absolute 
magnitude, color, and isolation. From modern galaxy surveys, in particular the SDSS, it is now 
possible to build up a statistical sample of Milky Way-like galaxies along with their respective 
satellite populations. To understand better how these studies work, begin by considering satellites 
brighter in absolute magnitude than the Magellanic Clouds. Modern kinematic measurements 



indicate that the Magellanic Clouds have maximum circular velocities > 60 km/s (van der Marel 



et al. 


2002 


Stanimirovic et al. 


2004) 



Way-mass halos in cosmological simulations contain two satellites with maximum circular velocities 



greater than the accepted maximum circular velocities of the Magellanic Clouds (Boylan-Kolchin 



et al. 


2010 


Liu et al. 


2011 



with masses greater than the mass of the Magellanic Clouds, but probably not uncomfortably so. 

In order to better probe the satellites discrepancy, it is necessary to extend this measurement 
down to fainter satellite magnitudes. However, because it becomes increasingly more difficult to 
reliably assign redshifts to satellites that are several orders of magnitude fainter than the Milky 
Way, the sample of host galaxies that can be used in the measurement dramatically decreases. 



Using SDSS data release 7, Guo et al. (2011) measure the satellite luminosity function down to the 
magnitude scale of Fornax, which is about eight magnitudes fainter than the Milky Way. These 
authors show that Milky Way-like galaxies have on average about a factor of two fewer satellites 
than our Galaxy or M31. Using SDSS data release 8, and a combination of photometric and 
spectroscopic redshifts, Strigari and Wechsler (2012) translate this into an upper limit of ~ 13 



satellites brighter than the luminosity of Fornax around Milky Way-like galaxies (Figure [8]) . At 
this point, these results broadly indicate that, down to at least the magnitude of Sagittarius, and 
probably down to the magnitude of Fornax (see Table|3]), the Milky Way and M31 are not significant 
outliers in terms of their respective numbers of satellites. As discussed again in section 10, in the 
next several years these measurements will significantly improve because of new larger scale galaxy 
surveys, thereby allowing us to study the missing satellites issue in a complete statistical framework. 

• Gravitational Lensing 

Strong gravitational lensing is an interesting technique that may be utilized to measure the 
distribution of substructure in galaxies and to study the nature of dark matter. In a strong lensing 
system, there are two basic observables, the positions of the images that are created, and the 
ratios of the fluxes between the images. If the lensing potential is smooth, then relations can be 
established between the fluxes of the separate images. A deviation from a smooth model for the 
potential may then reveal the presence of a dark matter subhalo, even though the galaxy associated 
with it is not identified in optical surveys. 

The largest sample of strong lensed systems comes from the Cosmic Lens All-Sky Survey 



(CLASS) (Myers et al. 2003). CLASS identifies lenses where multiple images of background radio 
sources are formed; radio sources are beneficial because they are free from dust extinction and 
they do not suffer from stellar microlensing. Interestingly, almost all of these lensing systems have 
images that violate the expected flux ratios for a smooth potential model. Several authors have 
argued that substructure, as is expected in LCDM models, causes the flux ratios to differ than 



expected in the smooth model (Mao and Schneider, 1998; Metcalf and Zhao, 2002) 
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Figure 8 Mean number of satellites brighter than an apparent magnitude around Milky Way-like 
galaxies, versus the satellite apparent magnitude. Daimonds (blue) and squares (black) are from 
SDSS data, and triangles (red) are Milky Way data. From Strigari and Wechsler (2012). 



Utilizing simulations of dark matter halos that have similar properties to the lensing galaxies 
observed, several authors studied whether it is possible to reproduce the flux ratios in the LCDM 



cosmological model (Xu et al. 2009 2010). The general conclusion of these studies is that there is 
too little substructure in LCDM halos to account for the observed flux anomalies. These authors 
also find that, since most of the lensing signal is due to subhalos greater than approximately 10 6 
Mq, extrapolating the mass function to lower subhalo masses is unlikely to affect the conclusion. 
The addition of baryons in the subhalos, and the corresponding contraction in the dark matter 
distribution that makes the subhalos more dense, also does not significantly affect the results. 
Current studies are now focused on the contamination of low mass galaxies along the line-of-sight 
towards the lensing system. 

Detailed modeling of individual lensing systems now uncover galaxies that are too faint to be 
seen in optical light. Vegetti et al. (2012) report the detection of a ~ 10 8 M@ mass satellite in 
the Einstein ring lens system JVASB1938+666, which is an elliptical galaxy at redshift z = 0.88 
that is lensing a bright infrared background source galaxy at redshift 2.059. This follows up on the 
eariler detection of a ~ 3 x 10 9 M & substructure associated with the lens system using a Bayesian 
reconstruction method (Vegetti et al. , 2010). In the former case, the authors utilize this detection to 
make the first measurement of the substructure mass function beyond the Local Group, finding that 
dN/dM cx M~ a , with a = 1.1±^;|. They determine the fraction of the total mass in substructure 
to be ~ 3.3%. Given the uncertain mapping between mass and luminosity at this mass scale, it is 
difficult at this stage to predict the luminosity of this galaxy, though it does sit in the same mass 
range as the classical dSphs of the Milky Way. 

Assimilating the results presented above, direct satellite counting experiments seem to be indi- 
cating that the Milky Way is relatively normal amongst galaxies of its kind in terms of its population 
of satellites. While promsing, the gravitational lensing observations are still probably too sparse 
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Figure 9 The circular velocities for the subhalos with the most massive progenitors, in a CDM 
simulation (left) and a WDM simulation (right). The data points are the masses within the half- 
light radius for the brightest dwarf spheroidals. Figure from Lovell et al. (2012); see larger sets of 
related figures in Boylan-Kolchin et al. (2011). 



at this stage to deduce more general conclusions about the population of dark matter substructure 
around galaxies. Taking the former results at face value, the evidence is that the satellite abundance 
issue points to an inadequacy of modern theory, either in the understanding of galaxy formation, 
the properties dark matter, the standard cosmological model, or some combination. 

5. 4-1-3. Variation of the LCDM Model. Finally, it is interesting to ask in more detail what mod- 
ifications to the physics of dark matter provide the best solution to this problem. As mentioned 
above, the most basic solution involves giving the dark matter larger velocities at freeze-out than 
in the CDM framework. In this case, the formation of dark matter halos with mass scale below the 
corresponding free-streaming scale of the particle is significantly suppressed. These particles may 
be produced thermally or non-thermally, and well-motivated models typically require the particle 



mass to be approximately the keV scale (Abazajian, 2006 Abazajian et al. 2012). Recent N-body 



simulations of Milky Way-mass galaxies find that a particle of this mass is effective in reducing 
both the total population of Milky Way subhalos and their central densities. As a result, the most 
massive dark matter subhalos are now able to host the observed dSph satellite galaxies, as indi- 
cated in the right hand panel of Figure [9] ( |Lovell et al. 2012). A second class of models involves 
introducing a larger-than-standard self-interaction cross section, as discussed in Section [T] (Spcrgcl 
and Steinhardtl |2000[ [Loeb and Weiner[ |2011| ). Several groups have shown via N-body simulations 
that subhalos in this model are also less dense, and it is also more likely for the dSphs to populate 



the most massive subhalos of the Milky Way (Vogelsberger et al 



et al. , 2012 ) 



2012 


Rocha et al. 
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5. 5. The mass of the Milky Way 

Very recently, an additional aspect of the missing satellites issue has come to the forefront that 
is related to our uncertainty on the mass of the Milky Way. Specifically, if the Milky Way mass is 
approximately 10 12 Mq, the number of subhalos with maximum circular velocity > 30 km/s can 
be made consistent with observations (Wang et al. 2012). Though this "low Milky Way mass" 



solution is appealing, there are still some aspects to it that require further investigation. One of 
which is that reducing the mass of the Milky Way to this level seems to make the existence of 
the Magellanic Clouds even less probable than the ~ 5% chance discussed above. Further direct 
measurements are clearly required. 

Cosmological simulations discussed in this section have been utilized to estimate the mass of 
the Milky Way. These estimates rely on identifying a Milky Way or Local Group-like system in a 



large volume cosmological simulation, and then applying a version of the "timing argument" (Kahn 



and Woltjer 1959). The mass of the Local Group is determined by identifying Milky Way-M31 



like systems in the simulation, while the mass of the Milky Way is determined by identifying Milky 



Way-Leo I like systems. Using the millenimum simulation, Li and White (2008) use this method 



to determine a total Local Group mass of approximately 5 x 10 12 Mq, and a 95% c.l. lower limit 
on the Milky Way mass of 8 x 10 11 Mq, which is consistent with the low Milky Way mass solution 
to the satellite abundance problem. This result is compared to the mass measurements from the 
kinematics of halo stars, globular clusters, and dSphs in section [3j 

The mass of the Milky Way may also be estimated from properties of satellites other than 



Leo I. Busha et al. (2011) recently construct the probability distribution for the Milky Way mass 



by locating the number of dark matter halos that have two satellites that are as massive as the 
Magellanic Clouds. Their results indicate a 68% c.l. lower limit on the mass of approximately 
5 x 10 11 Mq, which is also both consistent with many other mass measurements and the low-mass 
solution. While it would be ideal to extend this analysis to less massive satellites, both the lack of 
resolution in large scale N-body simulations and the uncertainty in the statistical mapping between 
the luminosity of a satellite and its dark matter mass preclude an analysis at this time. 

5.6. Caustics 

This section closes out with a brief discussion of an additional property of dark matter halos that 
is intriguing from both the perspective of their phase space properties and for the detection of par- 
ticle dark matter. For cold and collisionless dark matter, since the phase space density is conserved 
along the particle trajectory, the collapse of halos is expected to give rise to caustics. For com- 
pletely cold particles with zero velocity dispersion, particles effectively occupy a three-dimensional 
sheet in the total six-dimensional phase space, and the projection of the sheet into configuration 
space results in infinite density singularities. In reality these singularities are smoothed out because 
particles are cold and do not have precisely zero velocity dispersion. 



Caustics have been studied a great deal in the context of the spherical infall model (Fillmore 



and Goldreich, 1984 Bertschinger 1985). Within the context of hierarchical structure formation, 



caustics are also expected to form as a tidal stream folds back upon itself. It is not yet possible to 



resolve caustics in modern N-body simulations due to numerical resolution (Diemand and Kuhlen 



2008 ) , however analytical schemes have been developed that allow for the prediction of annihilation 



radiation from caustics (Hogan, 2001), including the smoothing of them through the finite velocity 



dispersion (White and Vogelsberger 2009). The issue of detection of annihilation radiation from 



caustics will be addressed again in the next section, in particular comparing the expected signal to 
the more often considered signals from the smooth halo and from halo substructure. 
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6. Indirect Detection: Theoretical Developments 

As discussed in Section [TJ particles with roughly weak scale interactions that were in thermal 
equilibrium in the early universe may decouple with an abundance to constitute much if not all 
of the observed dark matter in the universe today. Via the same processes that set their relic 
abundance, dark matter particles annihilate today in regions of high density, for example in the 
centers of galactic halos over a wide range of mass scales. The particles that are produced in 
the annihilation processes may be detectable through a variety of modern ground and space-based 
observatories (Section [7]). 

Discussion of detection of standard model particles produced in the annihilation of dark matter 



dates back more than three decades (Gunn et al. 1978). Since then, the idea that particle dark 



matter may be detected through astrophysical observations has inspired an immense amount of 
theoretical work. From the theoretical perspective, the most pertinent questions are: Where is the 
dark matter? What standard model particles does it produce, and what is the rate of standard 
model particle production? The answer to the former question may be addressed through the 
observations and analysis methods discuss in Sections [3} [5] The latter question is addressed by 
appealing to a specific particle physics framework, several of which are introduced for the first time 
in this article in this section. 

On the theoretical side there has also been interest in models beyond the standard WIMP frame- 
work. For example, dark matter particles that are unstable, such as sterile neutrinos or lightest 
particles in supersymmetric theories that violate R-parity, provide unique phenomenological signa- 
tures both from the perspective of astrophysics and particle physics. It has also been understood 
that dark matter may be detected indirectly in searches over a wider spectrum of energies than are 
considered in searches for standard WIMPs. 

This section reviews recent theoretical developments in the area of indirect dark matter detec- 
tion, bringing together the astrophysical and particle physics frameworks that are the most relevant 
for interpretation of modern experimental results. In introducing a representative set of particle 
models, the goal is to make the most straightforward possible connection between the cross section 
that is relevant for dark matter production in the early universe and the cross section that is most 
relevant for detecting dark matter via indirect searches. 

6.1. Dark matter annihilation 

Along a direction £, the observed flux of photons from dark matter annihilation above an energy 
Emin per solid angle is 

^ = ^%^f<-»- < 57 > 

where ^ represents the angular separation from the center of the galaxy, D is the distance to the 
center of the galaxy, so that r 2 = £ 2 + D 2 — 2£D cos\I/ . The upper and lower boundaries to the 



integral are £± = D cos ± y r 2 — D 2 sin 2 ^, where rt is the tidal radius of the dark matter halo. 
The number of photons produced is the integral over the energy spectrum of decay photons, 

rM dm dN 

»■> = L ^ dE - (58) 



Equation 59 neatly separates into a component that depends on the properties of the dark matter 
particle, A^y(o" ann v)/(2Mj m ), and a component that depends on the distribution of the dark matter 
in the halo. Integrating over a solid angle Af2 gives the total flux within the solid angle 

£(Afi) = N, 2 [ -j0^-27r / sinM / p 2 [r(£)]d£. (59) 

dm JO J t— 
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Explicitly separating out the component that depends on the dark matter distribution we have 



J (AO.) 



1 



An 



sin'M I p 2 [r(£)]d£. 



(60) 



Equation 60 will be referred to here as the "J-value." 

Annihilation radiation is expected to arise from three primary components of a halo: the smooth 
component, subhalos, and caustics. The smooth component is straightforward to model under the 
assumption of a dark matter density profile (e.g. Equations 16 and 31). The contribution from 



halo substructure and caustics requires a more detailed theoretical treatment. 
• Substructure 

If there is significant substructure in the dark matter distribution, as is predicted by N-body 
simulations, there is a corresponding "boost" to the predicted flux over this smooth component []] 
There are several papers in the literature that provide a theoretical treatment of the substruc- 



ture boost factor (Strigari et al. 2007d Pieri et al. 2008 Martinez et al. 2009 Charbonnier 



et al. , 2011 ), which complement the results derived from N-body simulations above their resolution 
limit (Springel et al. 2008b Kuhlen et al. 2008). The following discussion closely follows the 



methods outlined in Strigari et al. (2007d) and Martinez et al. (2009). 



Start from the definition of the luminosity of a source in Equation 59 
from a galaxy, including the substructure boost factor, as 

C{M h ) = [1 + B(M h ,M mm )]C(M h ), 



Define the total flux 



(61) 



where B is defined as the boost factor, which is zero for a totally smooth dark matter distribution, 
and C is the contribution from the smooth halo. Here the boost is defined to depend on the total 
halo mass, Mh, and the minimum mass in the substructure hierarchy, M m j n . While the former 
quantity can be estimated from observations of a given source, the later is of course only weakly 
theoretically-constrained. 

Due to the mass-dependence of the radial distribution of substructure discussed in Section [5j 
there is predicted to be a radial dependence of the substructure boost factor, which depends strongly 
on how efficiently subhalos of a given mass are destroyed by tidal interactions within the host halo. 
Simulations show that the substructure boost factor increases towards the virial radius of the 



halo (Springel et al., 2008b), reflecting the distribution of the most massive subhalos in a host 



galaxy. Smaller mass subhalos contribute progressively more to the annihilation luminosity in the 
inner regions, though because of the resolution of modern simulations at this time the boost factor 
within the approximate scale radius of a galaxy can only be determined from extrapolation of the 
results in the outer region of the host galaxy. Because of this theoretical uncertainty in the radial 
dependence of the boost factor for subhalos below the mass limit of modern simulations, and also 
because the angular resolution of modern experiments implies that they are likely not sensitive 
to the radial dependence of the boost factor (except perhaps for a couple of targets such as M31 
and nearby clusters), the primary calculations in this section focus on the total enhancement from 
subhalos integrated over the entire primary halo. 

To determine the value of the boost, annihilation luminosity factors C{m) for subhalos within 
mass m of the host halo are integrated, BC(M) l ) = J (dN/dm)C(m)dm. With the above assump- 



tions for no radial dependence of the boost factor, from the definition in Equation 61 this boost 



1 Note that this factor is to be distinguished from increases in the annihilation cross section, such as from Som- 
merfeld enhancements 
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factor can be explicitly written out as 



B(M h 



qM h m 

C{M h ) J Mm . n dm 



C{m)dm 



C{M h ) 



qM h 



dN 
dm 



1 + B(m)]C(m)dm 



AM a' f qM h 



C{M h ) 



[1 + B(m)\£(m) m 



-l-a' 



dm. 



(62) 



where a subhalo mass function of dN/dlnm = A{Mj l /m) a is assumed, and q reflects the fact that 
the subhalo mass function continues up only to a fraction of the mass of the host galaxy, typically 
q ~ 0.1. From the results in Section [5j a 1 = 0.9. For a fixed a', the normalization factor A is 
determined by demanding that a fraction of the total mass in substructure is within a given mass 
range. 



To proceed in the evaluation of Equation 62 we need a model to relate the luminosity of a 
subhalo to its mass. For the purposes of this analysis approximate this C — m relation as an 
analytic relation, though its probably more realistic that there is scatter in this relation that can 



ultimately only be probed through higher resolution numerical simulations (Springel et al. 2008b 



Kuhlen et al.[ 2008 ) or Monte Carlo analysis ( Martinez et al. , 2009 ) . Motivated by the concentration 
mass relation (e.g. Bullock et al. (2001)), the subhalo luminosity-mass relation has been shown to 
be well- approximated by a power- law relation, £ oc M 87 ( [Strigari et al. 2007d Kuhlen et al. 



2008) 



It is clear to see that the substructure boost factor is a strong function of the halo mass M^. For 
example, for a typical galaxy cluster, the boost factor may be in the range B ~ 10 2 — 10 3 depending 
on the minimum mass cut. This significantly improves the detection prospects of clusters. For 
galaxies about as massive as the Milky Way, such as M31, for an extrapolation of M m i n down 
to Earth mass, B ~ 10 2 . On the other hand, for dSphs, B is likely only a factor of few at the 
most, which in the context of the above calculation is a reflection of the small host halo mass for 
dSphs. If the radial distribution of the boost is also considered, it is also likely that the boost 
is insignificant for dSphs, because it is plausible that subhalos in the outskirts of the dSphs are 
efficiently stripped off by the Milky Way halo. However at this stage these the arguments are 
heuristic because theoretical knowledge of dark matter substructure in Milky Way-mass galaxies, 
much less sub-substructure distributions that may be present in dSphs, is subject to significant 
uncertainty. Though the above arguments do not improve the detection prospects of dSphs, they 
do motivate more of a clean interpretation of the annihilation cross section limits that are derived 
from the data, which are discussed in Section [7j 



Caustics 



As discussed in Section [5] in the self-similar halo model caustics will form in regions where the 
mapping between phase space and physical space is singular. Caustics are also expected to form 
when a satellite is tidally-destroyed upon accretion into a parent halo. Because of the high density 
achieved at caustics, the dark matter annihilation rate is expected to be enhanced relative to the 
smooth dark matter halo contribution. This enhancement can be quantified by noting that, on one 
side of the caustic, the density diverges as p ~ A/ ^fx, where x is the distance from the center of 
the caustic and A is the fold coefficient. 

Though the above relation shows that the density diverges for purely cold (i.e. zero velocity dis- 
persion) particles, in practice the dark matter is not purely cold but has a finite velocity dispersion. 
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For a WIMP of mass 100 GeV, conservation of phase space density gives the velocity dispersion 
today for particles that are just infalling into the halo as ~ 0.03 cm/s. This finite velocity dis- 
persion cuts off the divergence of the caustic, so that the dark matter annihilation rate per unit 
surface area is J p(£)d£ ~ f^A 2 /£ oc ln(R/5£), where R is the radius of the halo. This predicts 
that the increase in the annihilation rate over the smooth halo due to the caustic is a factor of 



several (Hogan, 2001; Natarajan and Sikivie, 2008). 



6.2. Annihilation cross section models 

With the input from astrophysics in place for calculating the gamma-ray flux from dark matter 
annihilation, it is next necessary to evaluate the annihilation cross section in some representative 
particle models. On general theoretical grounds, particle dark matter candidates can be classified 
according to their intrinsic spin. Scalar particles have spin and fermions have spin 1 /2 (Higher spin 



dark matter particle models have also been studied (Barger et al. 2008 Yu et al. 2012)). Fermion 
dark matter candidates can be further classified as either Dirac or Majorana particles, which differ 
from one another in the interactions that are allowed. With these definitions interaction operators 
can be constructed between the dark matter and standard model particles that transform as odd 
or even under charge, parity, and time reversal transformations. This approach is valid under some 
simplifying theoretical assumptions, most importantly that the dark matter and standard model 
interactions are mediated by new and unknown physics above a certain energy scale. Models along 
these lines are sometimes called effective theory models, and have recently been examined in the 



context of WIMP direct and indirect detection by several authors (Beltran et al. 2009 Cao et al 



2011 


Bai et al. 


2010 


Zheng et al. 


21)12 


Yu et al. 


2012) 



To define a representative set of models, assume that the WIMP couples only to standard model 
fermions. For a scalar WIMP, the interactions between the WIMP, represented by the field 4>, and 
a fermion, represented by /, are 



Scalar 
Vector 

Scalar-Pseudoscalar 
Vector-axialvector 



Fsp 
V2 



00/75/, 



^(^0)/7 M 75/- 



For fermionic WIMPs, represented by x, the corresponding interactions are 

G s 



Scalar 
Pseudoscalar 
Vector 
Axialvector 
Tensor 



xxff 



V2 

G P _ .-, 



V2 
Gy 

71 



X7 x/75/ 



^j=xi^i b xfi^f 

Gt 



V2 



X^ v xfo-u V f. 



(63) 
(64) 
(65) 
(66) 

(67) 
(68) 
(69) 
(70) 
(71) 
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Note that the various couplings that appear in these equations have different dimensions, and in 
general the couplings to fermions could be universal (independent of fermion mass), or dependent on 
the mass of the standard model fermion. It is also possible to define mixed interactions for fermionic 
dark matter, i.e. scalar-pseudoscalar interactions of the form XXflsf i m which the lagrangian is 
invariant under the combination of charge, parity, and time reversal transformations. 

From the interactions above it is possible to define the annihilation cross sections, and expand 
them as (J &an v ~ a + bv 2 . For scalar dark matter these are given by ( Beltran et al. , 2009 , Zheng 



et al. 


2012 


Yu et al. 


2012) 



(o"ann«)s - 

(o-ann'w)v - 
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m. 



ML 



2 + 



f 

ML, 



4vr Vv^y 
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(72) 
(73) 
(74) 
(75) 



and for fermonic dark matter, 
(o' aim v)p 

(0"annW)y 
(0"ann«)A 

(o"ann^)T 
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m| \ 3 / 2 



M L. 
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m. 



(M 2 dm + 2m 2 ) 

dm 



(76) 
(77) 
(78) 
(79) 

(80) 
(81) 



In the formulae above for the annihilation cross section, only leading order terms (a or b) are 
shown. The color factor is Cf = 3 for quarks and c/ = 1 for leptons. The above expressions 
show that scalar dark matter with scalar interaction and psuedo-scalar interactions with fermions 
annihilate through s-wave interactions, while scalar dark matter with vector and vector- axial vector 
interactions with fermions annihilate through p-wave interactions. Fermionic dark matter with 
pseudoscalar, vector, and tensor interactions annihilate through s-wave interactions, and fermonic 
dark matter with scalar or axial vector interactions annihilate through p-wave interactions. For 
majorana fermions, the vector current interaction vanishes. 

The above models represent only a subset of effective models because they focus on interactions 



with fermions. Models similar in spirit with couplings to gauge bosons are considered in, e.g. Cao 
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et al. (2011). A simple model that is related to those presented above is the singlet scalar model, 



in which a single dark matter particle represents the only beyond the standard model physics (Mc- 



Donald 1994; Burgess et al. 2001). Depending on its mass, the dark matter can annihilate either 
through a Higgs, gauge bosons, or fermions. For each of these different models, the annihilation 
cross section is 

£,r(fc->x t ) 



( 4M J 



(82) 



where v ew is the vacuum expectation value of the Higgs, ran is the mass of the Higgs, is the 
total Higgs decay width, and T(h —> X t ) is the partial decay width into the final state X. Note 



that, in contrast to the scalar WIMP with scalar interactions with fermions in Equation 72, in the 
singlet scalar model the annihilation cross section to fermions is proportional to the fermion mass 
because of the coupling to the Higgs. 

In supersymmetric models, neutralinos can annihilate into fermions through exchanges of Higgs, 
Z° bosons, and through sfermions. In all these cases the annihilation cross section is proportional to 
the mass of the outgoing fermion. There is a significant amount of literature devoted to computing 
neutralino annihilation cross sections; formulae for the annihilation cross section for neutralinos 



into fermions, as well as for gauge bosons and Higgs, have been compiled in Bertone et al. (2005). 
In addition to the continuum spectra, neutralinos may also annihilate to a two-photon final state or 
a final state with a photon and a Z°. In the later case, the energy of the photon is mono-chromatic, 
with an energy given by 



E. 



7 



Mdn 



rrir 



4Hn 



(83) 



The annihilation cross sections for these processes are dominated by the s-wave contributions, and 
for the neutralino parameter space maximize at approximately 10 -29 cm 3 s _1 (Bergstrom and Ullio 



1997; Ullio and Bergstrom, 1998; Bern et al. 1997), nearly three orders of magnitude below the 



thermal relic scale. 

To calculate the relic density of dark matter and the flux from particles in halos today, the 
annihilation cross sections must be averaged over the velocity distribution, (cr ann w). For the relic 



abundance calculation, this is performed via the methods outlined in, for example, Gondolo and 



Gelmini (1991), Griest and Seckel (1991), and Edsjo and Gondolo (1997). If the annihilation 



cross section does not depend on the relative velocity (a 3> b), as is the case for scalar dark 
matter with scalar and pseudoscalar interactions, and fermionic dark matter with pseudoscalar, 
vector, axialvector, and tensor interactions, then the averaged annihilation cross section in the 
early universe is the same as it is in halos today. In the opposite limit ((» > o), the annihilation 
cross section is lower in halos today than it was in the early universe, because the typical velocity 
in a halo is v/c ~ 10 -3 . In a halo the averaged annihilation cross section depends on the assumed 



halo density profile and velocity distribution, and also on the position within the halo (Campbell 



et al., 2010 Frandsen et al., 2012a) 



6.3. Photon and charged particle production 

The detectability of a model characterized by the annihilation cross sections presented above 
depends on the spectrum of photons produced through the decay or fragmentation of the fermions 
or gauge bosons in the final state. The main photon production channel is through production 
and subsequent decay of neutral pions. The photon spectra are typically determined numerically 



through programs such as Pythia 6.4.19 (Sjostrand et al. 2006) and there are fitting formulae 



that describe the output (Fornengo et al. 2004; Cembranos et al. 2011). For continuum spectra 
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the number of photons produced in the annihilation peak at an energy approximately an order of 
magnitude below the dark matter mass. 

In addition to continuum spectra, for dark matter that annihilates into final state fermions, 



photons from final state radiation (FSR) are produced (Beacom et al. 2005 Bergstrom et al. 



2005). FSR photons are produced whenever there are charged particles in the final state. The 



spectrum of photons from FSR is distinct from the line spectrum and the continuum spectrum in 
that it peaks near the mass of the dark matter ( Birkedal et al. , 2005 ) . When there is a fermion in 



the final state, the photon energy spectrum from final state radiation is (Essig et al. 2009) 



dN~ ( 
dy 



a 

7T 



y 



in 



s(i-y) 



in 



f 



(84) 



Here a = 1/137, y = Ej/Ma m , s = 4Mj m , and the maximum photon energy is M dm . Note that 
for FSR, the number of photons produced is suppressed by a factor of a relative to the continuum 
emission. 

If electrons and positrons are produced from the dark matter annihilation process, secondary 
photons that are produced by these e + e~ pairs may also be detectable. The electron-positron 
population is non-thermal, and loses energy through a variety of physical processes, including: 
synchrotron from magnetic fields, inverse Compton radiation from interaction with both starlight 
and the cosmic microwave background, and classic bremsstrahlung via interactions with interstellar 
gas. Theoretical predictions for the corresponding photon fluxes, as well as the fluxes of electrons, 



positrons, and protons, are reviewed in Cirelli et al. (2011) and Profumo and Ullio (2010) 



6.4- Dark matter decay 

In Section [7j experiments that are sensitive to photons produced from particle decays will also 
be discussed. Though there is no scale corresponding to the preferred thermal relic scale in the case 
of annihilations to preferentially set the abundance of dark matter, a large class of well-motivated 



models predict observable signals from dark matter decays ( Arvanitaki et al. , 2009, Abazajian et al 



2012) 



For decaying dark matter models, the flux corresponding to Equation 59 



is 



dn 



1 



1 



M dm T 47T 



p[r(£)]d£, 



(85) 



where r is the lifetime of the particle. As in the case of annihilation, Equation 85 nicely separates 
into a piece that depends on the properties of the dark matter particle, 1/Mrf m r, and the remainder 
which depends on the properties of the host dark matter halo. Because in this case the flux is 
proportional to the density and not the square of the density, boost factors from halo substructure 
are not relevant. The total flux within a solid angle Af2 can be then calculated in a manner similar 



to Equation 60 



The fact that the signal scales as the density and not as the density squared implies that the 
signals are more spatially-extended. The decaying dark matter signal from our Galaxy and from 
nearby galaxies is rivaled in some instances by the signal from diffuse background radiation from 
extragalactic sources, which in many instances constitute a background to the signal that is to be 



observed (Cembranos et al., 2007; Cirelli et al. 2012) 



6. 5. Neutrinos 

In addition to the above discussed signatures from high energy photons, indirect signals of dark 
matter annihilating or decaying may be detected through other channels. One such interesting 
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channel is neutrinos. By their nature, neutrinos are much more difficult to detect than photons, 
however, they represent a complementary probe to photons for two reasons. First, they are produced 
similar to photons in the fragmentation or hadronization products of annihilating dark matter. 
Second, they propagate from the source without interacting, so they are able to probe dark matter 
produced in different environments than photons. There are two primary methods utilized to search 
for neutrinos from dark matter annihilation. First, through their production in the Galaxy and 
external galaxies as discussed above. Second, they are produced by dark matter that gets captured 
and annihilates in the Earth and the Sun. Since the astrophysical phenomenology of neutrinos 
produced in the Galactic halo or external galaxies is similar to that discussed above in the context 
of gamma-rays, the remainder of this section briefly reviews the production of neutrinos from dark 
matter that is captured in the Sun or the Earth, anticipating results to be discussed in Section [9j 



The capture rate, C, in a differential volume element of the Sun can be written as (Gould 



1987) 



dC 
dV 
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. m 2 
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dE 
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(86) 



/ n a J e 

where n is the number density of dark matter near the Sun, and a is the cross section for a dark 
matter particle interacting off of a nucleus. The velocity in a given shell, uu, is related to the velocity 
at infinity u and the escape velocity, v esc {r), in the shell at radius r, as w = yt 
upper limit to the integral in Equation [86] is defined as 
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where m A represents the mass of a nuclear species. The upper and lower limits to the energy 
integral are, 
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where E ma x is the maximum energy for two-body elastic scattering, and E m i n is the minimum 
energy required for capture. The minimum energy corresponds to the capture of a WIMP into an 
orbit at any radius; requiring the WIMP to be captured within the radius of Jupiter changes the 
value of Emin ( |Kumar et al. 2012a). 

The flux depends on a competition between the capture rate of WIMPs in the core of the Sun, 
determined from Equation 86, and the annihilation rate times the relative velocity divided by the 
effective volume of the core of the Sun. For dark matter in equilibrium, the annihilation rate is 
one- half the capture rate. The time to reach equilibrium for dark matter in the Sun is 



i© ~ 8.6 x l(T 2 yr 
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where Tc is the dark matter capture rate. Assuming that dark matter annihilation and capture are 
in equilibrium, neutrinos from the Sun are primarily sensitive to the scattering cross section. The 
spectrum of neutrinos produced per annihilation is generated from Pythia; fitting formula for dark 



matter annihilation into quarks, leptons, and gauge bosons are presented in Cirelli et al. (2005) 



Low-energy neutrinos that are produced from annihilation and decay products that lose energy 



through interactions in the solar medium may also be detectable (Rott et al. 2012 Bernal et al. 
2012b. 
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From the perspective of Galactic astrophysics, the capture rate in the Sun is sensitive to the 
lowest velocity of WIMPs in the velocity distribution. Note that the situation is much different for 
capture in the Earth, which is more strongly affected by diffusion and interactions with the planets 
in the Solar System and acceleration from the Sun. This is particularly important at low veloc- 
ity (Lundberg and Edsjo, 2004). This is in contrast to direct searches, which are primarily sensitive 
to the dark matter in the high velocity tail of the distribution. So, though the neutrino rates are a 
more complicated convolution between the cross section and the dark matter velocity distribution 
(Equation 86), neutrino searches and direct dark matter searches also give complementary probes 
of the dark matter velocity distribution in the Solar neighborhood. 
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7. Indirect Detection: Experimental Results 

The previous section has reviewed the theoretical aspects of indirect dark matter detection, 
highlighting the role of both astrophysics theory and particle dark matter theory. The most exciting 
challenge confronting research in indirect dark matter detection involves understanding how to test 
theoretical models with modern experimental data. Coming at it from the other direction, it is also 
important to gain an understanding of what information can realistically be extracted from the 
modern experimental data sets that is most relevant for dark matter. With these general themes 
in mind, this section moves on to discussing the status indirect dark matter searches from a variety 
of modern astronomical data sets. 

Before going into details, it is necessary to take a step back and establish the scope of indirect 
dark matter searches. By their nature, these searches typically do not reflect "controlled" experi- 
ments, in the sense that often it is difficult to simulate experimental backgrounds and get repeated 
measurements of them. In nearly all cases, backgrounds or foregrounds to the sought after signal 
result from physical processes that are poorly understood. 

Since indirect dark matter searches are typically undertaken by experiments motivated by other 
astrophysics topics outside of the realm of dark matter, it is useful to think generally about their 
strategies when it comes to searching for dark matter. A first type of search strategy is designed 
to identify dark matter from either a known type of target or from a distinct feature in the dark 
matter spectrum. In other words, we a priori know what type of signal we are looking for, and/or 
where to find it. There are two examples that illustrate this search strategy. First, gamma-ray 
studies of dSphs, which have no intrinsic gamma-ray backgrounds and relatively well- understood 
dark matter distributions. Also in this category are gamma-ray line searches in regions of expected 
high dark matter density. Since no known astrophysical process can produce sharp lines at energies 
of hundreds of MeV and above, this is a distinct spectral feature that is difficult to mimic without 
dark matter. A defining characteristic of these types of searches is that they are based upon existing 
theoretical predictions. For example in the case of the dSph searches they are based on determi- 
nations of the dark matter distributions from stellar kinematics, combined with predictions for the 
annihilation cross section from the freeze-out arguments. Or in the case of the line searches, they 
rely on the determination of the intrinsic strength of the line signal from a theoretical calculation. 

A second type of indirect search strategy relies on more of an empirical approach. It revolves 
around analysis of astronomical data sets, identification of features in these data sets that may 
be interpreted as a dark matter signal, and building a theoretical model around it. This line of 
thinking has given birth both to surprising experimental results and creative theoretical models, 
and is an important part of the process because it ultimately highlights our theoretical ignorance as 
to the nature of particle dark matter. It is also important because, as we have learned from many 
areas of astronomy in the past, the most ground breaking discoveries are often those not at all 
anticipated. An example that defines this category is the interpretation of the recent results from 
the PAMELA satellite, which finds a rise in the local positron fraction up to energies of hundreds 
of GeV. An additional more recent example is the possible existence of a feature in the Fermi-LAT 
data near the Galactic center at approximately 130 GeV. Though this has generated a great deal 
of theoretical interest, it was certainly not predicted a priori from a theoretical model. 

The primary motivation for the above discussion is that it provides a context to interpret the 
results of indirect dark matter searches, and allows us to make better sense of the myriad of results 
that now flood the literature in this topic. As for the results that are covered in this section, 
the discussion begins by focusing on the most updated results from gamma-ray observatories, in 
particular on the Fermi-LAT and modern Air Cerenkov Telescopes (ACTs). Since we now have 
a tremendous amount of data from these instruments over all regions of the sky, we are able to 
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study essentially all classes of astro-physical sources from which a dark matter signal is expected. 
Following up on the gamma-ray analysis, the results and status of modern indirect searches that 
utilize neutrinos, cosmic-rays, and X-rays are discussed. This section closes with an interpretation 
of the most important experimental results, and the reach of future experiments in improving upon 
these results. 



7.1. Overview of gamma-ray experiments: Fermi-LAT and ACTs 



The Fermi gamma-ray space telescope was launched in June 2008 into low Earth orbit (Atwood 



et al. 2009). It was designed as the successor of the successful EGRET mission, with an order of 



magnitude better energy and angular resolution than EGRET. Fermi consists of two experiments, 
the Large Area Telescope (LAT), which is the primary instrument aboard Fermi, and the gamma-ray 
burst monitor. The LAT is sensitive to photons in the range of approximately 20 MeV — 300 GeV, 
while the gamma-ray burst monitor is sensitive to photons in the range of hundreds of keV. This 
article will focus solely on the results obtained by the LAT. 

The LAT is an imaging, wide field-of-view pair conversion telescope that measures electron 
and position tracks that result from the pair conversion of an incident high-energy gamma-ray in 



converter foils (Atwood et al. 2009). In between the converter planes are silicon strip detectors 



that are sensitive to charged particles that measure the tracks that result from the pair conversion. 
The calorimeters measure the electromagnetic energy shower produced by the electron-positron 
pair and the development of the shower. The single photon on-axis angular resolution varies from 
approximately 0.1° at energies greater than 10 GeV to approximately 3.5° at 100 MeV. The energy 
resolution, AE/E, over the energy range of interest in approximately 10%. The effective area of 
the LAT is approximately 10 4 cm 2 . 

Figure 10 shows an all-sky gamma-ray map in Galactic coordinates obtained from two years of 
Fermi-LAT data for photons in the energy range 100 MeV -300 GeV. Because of the steep nature of 
the power law spectrum of gamma-ray sources, most of the photons visible are in the energy regime 
of hundreds of MeV. As is clear from Figure 10 , the brightest source of gamma-rays in the sky is 
the diffuse emission along the Galactic plane. Diffuse gamma-rays from the Galactic plane in the 
Fermi-LAT energy bands are created by cosmic ray interactions in the interstellar medium. Over 
the energy range of hundreds of MeV to tens of GeV, the dominant component of the gamma-rays 
results from ttq decays. Above ~ 1 GeV, the next largest contribution to the diffuse gamma-ray 
flux arises from inverse Compton scattering of electrons from the interstellar radiation field. Below 
~ 1 GeV, the flux from bremsstrahlung is near the flux from inverse Compton ( Abdo et al. , 2009a). 

Bright sources in the Galactic plane, such as pulsars and supernova remnants, are clearly visible, 
as well as extragalactic sources at high latitude, in particular Active Galactic Nuclei (AGN). The 
two-year point source catalog of the Fermi-LAT contains over 1800 point sources, with over 100 
of them firmly identified in other wavelengths, and over 1100 associated with sources at other 
wavelengths (Nolan et al. 2012). Firm identification with sources at other wavelengths is typically 
done through timing and spectral information. This leaves over 30% of the point sources in the 
two-year catalog as unidentified at other wavelengths, with the fraction of unidentified sources 
increasing towards the Galactic plane, where confusion with diffuse Galactic radiation becomes an 
important issue. As discussed in more detail during the development of this section, understanding 
both the diffuse Galactic radiation and the properties of the sources in the two-year catalog will be 
important in our quest for dark matter signals. 

Complementing the Fermi-LAT searches for dark matter, there are now several ground-based 
imaging Air Cerenkov Telescopes (ACTs) that have rapidly improved in both size and sensitivity, 
giving them an important role in indirect dark matter searches. Generally, ACTs detect Cerenkov 
light from air showers that are produced in the atmosphere by cosmic rays. The shape and the 
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Figure 10 Two year all sky Fermi-LAT map of the energy flux. The bottom legend is the flux in 
units of 1CT 7 cm" 2 s _1 . FromlNolan et al.l (|2012l). 



orientation of the shower allows for a discrimination between those events generated by gamma- 
rays and those induced by other types of cosmic rays. The typical effective area for a modern ACT 
is much larger than that of the Fermi-LAT, approximately 10 5 m 2 , and the angular resolution of 
approximately 0.1 degrees over the appropriate energy regime is much better than that of the LAT. 
The lower energy thresholds are approximately hundreds of GeV, while the high energy threshold 
is much greater than that of the Fermi-LAT, up to tens of TeV. 

There are now hundreds of gamma-ray sources that have been detected and studied by modern 
ACTs, including blazars and pulsars. Most interestingly from the perspective of dark matter 
studies, ACTs are pointed instruments that are able to integrate for a long period of time on an 
individual dark matter source, as long as that source is within its field-of-view. However, relative 
to the LAT they are less sensitive to studies of the nature of the diffuse radiation over large regions 
of the sky. Further, as discussed in more detail below, given modern exposures (i.e. the integration 
time multiplied by the effective area) the sensitivity to a given dark matter cross section is typically 
a couple of orders of magnitude weaker than the sensitivity of the LAT. At present ACTs require 
tens to hundreds of hours of integration per source in order to obtain even this level of sensitivity. 
Table [5] provides some basic information on modern gamma-ray experiments that have presented 
dark matter results. 

With the above background, the remainder of this section is devoted to reviewing and inter- 
preting results from gamma-ray experiments. The results are organized according to the location 
and the classification of the source: 1) Galactic center, 2) dwarf spheroidals, 3), Galactic halo and 
extragalactic, 4) Galactic substructure, and 5) clusters of galaxies. 

7.2. Galactic center 

Due to its relative proximity, it is reasonable to assume that the largest flux of gamma-rays 
from dark matter annihilation comes from the Galactic center. However, as discussed in Section |3j 
the dark matter distribution in the Galactic center is poorly constrained from direct kinematic 
measurements, because in this region the stars in the bulge and the central nuclear star cluster 
appear to dominate the dynamics. Unfortunately this means that there is no reliable way to 
estimate the dark matter density in the central region of the Galaxy, and therefore there is no way 
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Table 5. Modern gamma-ray experiments that are active in indirect dark matter searches. 
Experiments below the horizontal line are ACTs. 



Experiment Effective Area [m 2 ] Energy Range [GeV] Location 



Reference 



Fermi-LAT 



1 



0.1 - 300 



Space 



VERITAS 
MAGIC 
HESS 



io 5 

2 x 10 4 
5 x 10 4 



50 - 50000 Arizona 
50 - 30000 Canary Islands 

10 - 10000 Namibia 



lAtwood et al.l d2009b 



A 



Holder et al.| (2006 
|20T2|) 



Aleksic et al 



mronian et al. (2006b) 



to a priori predict the precise photon flux from dark matter, and test this with the observations 
(even if we allow for the fact that we don't know the dark matter particle mass). 

To get an estimate of a theoretically-plausible prediction for the gamma-ray flux from the 
Galactic center region, consider an Einasto dark matter profile (Equation 50) with p_2 = 3.7 x 10 5 
Mq kpc -3 and r_2 = 16 kpc (Navarro et al. 2010). A dark matter halo with these parameters is 
consistent with both the measurements of the local mass density and the total Milky Way halo mass 
(Section[3j). With this set of parameters, within a region of radius of one degree around the Galactic 



center from Equation 59, J(l°) ^ 1.5 x 10 GeV cm 5 . For (a : 



3 x 10 



-26 cm 3 g 



a dark 



14. From Equation 



59 



matter mass of 100 GeV, and a bb spectrum, above energies of 1 GeV, Nj 

the photon flux above 1 GeV is F(l°) = 6.3 x 10~ 10 cm" 2 s . For comparison, the two year 
Fermi-LAT source catalog list five point sources within one square degree from the Galactic center. 
For energies greater than 1 GeV, the total flux from these point sources is ~ 1.6 x 10~ 7 cm -2 s -1 . 
So even though the predicted flux from a standard dark matter model is above threshold for the 
LAT, the point sources in this regime represent a significant background for the signal. 



At higher energies than is accessible by the LAT, there is now firm evidence from HESS (Aha- 



ronian et al] |2004[ ) and MAGIC ( |Albert et a\] |2006| for a point source, HESS J1745-290, that is 
coincident with the position of Sgr A*. The energy spectrum of this source is dN/dE ~ E~ , and 
there is a cutoff in the spectrum at approximately 15 TeV. In the two-year Fermi-LAT source cata- 
log, there is a source coincident with HESS J1745-290, with a best-fit power law of dN/dE ~ E~ 2 ' 3 
over the range of 100 MeV to 100 GeV. Above 1 GeV, the flux from this source is 7.7 x 10 -8 
cm" 2 s , making it the brightest source over these energies observed by the Fermi-LAT near the 
Galactic center. 

In addition to the above point sources, diffuse gamma-ray emission has been detected from the 
Galactic center region. HESS finds that this emission extends in both directions along the Galactic 
plane, and that it is spatially-correlated with giant molecular clouds in the central 200 pc of the 
Galaxy. The emission is attributed to the decays of neutral pions produced in the interaction of 



hadronic cosmic rays in the molecular clouds ( Aharonian et al. , 2006a ) . Characterization of the 



diffuse emission from the Galactic center measured by the Fermi-LAT is more challenging because 
of the aforementioned number of point sources. After carefully subtracting the emission associated 
with them, the diffuse emission is determined by assuming a model of the sources of cosmic rays, 
their propagation, and the properties of the interstellar medium. It is now standard to calculate 



these distributions is with the GALPROP code ( 


Strong et al. 


2000 


Moskalenko et al. 


2002 


Strong 


et al. 


2004; 


Porter et al. 


2008 


Vladimirov et al. 2011). With this modeling in place, it is then 



left to understand if any residuals remain that may be due to a dark matter signal. 

The combined effects of the point sources, diffuse emission, and uncertainty in the dark matter 
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Figure 11 Upper limits on the annihilation cross section into lines. The different sets of data points 



are for different assumptions on the dark matter density profile. From Ackermann et al. (2012a). 



distribution make extraction of a signal from dark matter difficult. In order to avoid these uncer- 
tainties, the HESS collaboration has recently searched for gamma-ray emission between angles of 
0.3 — 1.0° from the Galactic center, corresponding to a physical distance range of 45-150 pc. In ad- 
dition to avoiding contamination from gamma-ray sources in the Galactic center and the plane, this 
observational strategy is also beneficial because there is expected to be less theoretical uncertainty 
in the dark matter distribution over this region. For example, normalizing to the local density and 
the total halo mass, there is less than a factor of two difference between the corresponding densities 
of the Einasto profile and the NFW profile. However, the difference is two orders of magnitude 
when assuming a cored dark matter profile that matches both the local density and the halo mass. 
Assuming a NFW or Einasto profile, HESS finds that the upper limit on the annihilation cross 
section is (<7 ann u) — 10~ 25 cm 3 s _1 for masses greater than 1 TeV and annihilation into quark anti- 



quark pairs. Abazajian and Harding (2012) extend the limits on the HESS data into the Galactic 
center, finding an upper bound of a ann v ~ 10~ 25 cm 3 s _1 at a WIMP mass > 1 TeV for cusped 
dark matter profiles. 

The Fermi-LAT collaboration has published searches for gamma-ray lines within a 20° x 20° 



region around the Galactic center ( Abdo et al. 2010a Ackermann et al. , 2012a). To search for lines 



in the gamma-ray data, a sliding window is used over energies in order to account for the energy 
resolution. Over the different energy bins the data is fit to a power law plus signal model, where 
the index of the power law is allowed to vary and the signal represents the line feature. In the 
Galactic center, the Fermi-LAT collaboration reports no anomalous gamma-rays lines, and places 
upper limits on the annihilation cross section to lines that vary from approximately a few times 

1Q -29 



cm 3 s 



at 10 GeV to a few times 10" 



-26 cm 3 g 1 



at 200 GeV (Figure 11). There are factors 



of few uncertainty reported reflecting the uncertainty in the dark matter density. Gamma-ray lines 
from the Galactic center will be discussed again at the end of this section. 

Though the precise nature of the diffuse gamma-ray emission in the Fermi-LAT energy range 
remains to be characterized, an interesting morphological structure has come to light in recent 
analysis of the data. After subtracting out templates for the neutral pion, inverse Compton, and 
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bremsstrahlung emission, Dobler et al. (2010) first reported a detection of the "Fermi Haze," which 
extends 50° above and below the Galactic plane. Two characteristic features of this gamma-ray 
haze are that it has a long-to-short axis ratio of two, and the emission in this haze region is 
spectrally-harder than the other gamma-ray components in the Galaxy. The gamma-ray haze is 



consistent with the microwave haze discovered in the WMAP data by Finkbeiner (2004), which is 



now believed to result from synchrotron emission from the same population of relativistic electrons 
that are the source of the gamma-ray haze. A similar haze has also recently been announced by 

Further analysis has revealed that the LAT data is 



2012) 



the Planck collaboration (Ade et al. 
more consistent with a "bubble" like morphology (ISu et al., 2010) that has sharp edges. Though 



the emission is extended in a manner similar to would be expected from emission from particle 



annihilation in the dark matter halo ( Finkbeiner , 2004 ) , both the morphology and the edges seem 



to argue against a dark matter interpretation (see however Dobler et al. ( 2011[ )). 

Though there has been great progress in recent years in our understanding of gamma-ray as- 
trophysics from the Galactic center, it may leave readers of this article somewhat unsatisfied that 
a definitive limit (or detection) of dark matter has not been advocated for. While it is certainly 
true that cross sections much larger than the thermal relic scale are now clearly ruled out for 
any reasonable dark matter distribution in the Galactic center, reliably pushing down to the ther- 
mal relic scale with Galactic center observations will certainly require better understanding of the 
systematics that have been discussed above. 

7. 3. Dwarf spheroidals 

The dSph galaxies, whose properties were discussed at length in Section [4j provide excellent 
targets for indirect dark matter searches, and in particular gamma-ray searches, for several reasons. 
First, there are no gamma-ray point sources and there is no intrinsic diffuse emission associated 
with them, as is the case of the Galactic center. Indeed, the very strict upper limits on star 
formation processes and the gas content of these systems predicts a negligible gamma-ray flux 



from these sources (Grcevich and Putman, 2009). Second, they have well-understood dark matter 



distributions that are derived from the stellar kinematics, as was shown in detail in Section |4j 
Third, since the boost factor is predicted to be negligible in these systems, interpretation of a 
null-flux detection as a limit on (<7 ann t> ) is much more straightforward than it is for the cases of the 
Galactic center and larger-mass galaxies (Section [6]). 

The lack of gamma-ray emission from cosmic-ray processes or from stellar remnants is an im- 
portant point, so it is informative to be as quantitative as possible with this discussion. The most 



salient point can be seen in Figure 12 which shows the relation between the star formation rate 
and the gamma-ray luminosity for star forming galaxies (Ackermann et al. 2012b). Though it 



is not yet clear what this relation should be for galaxies with lower star formation rate than is 
shown, we can get a rough idea by extrapolating the linear relation to lower star formation rates. 
The extrapolation down to the scale of the only dSph with observed star formation activity, Leo 
IV, predicts a negligible gamma-ray flux from this object. So considering that other dSphs have 
negligible star formation, it is very likely that any gamma-rays observed from these objects arise 
from dark matter annihilation (Note that it is possible that several more star forming galaxies could 
appear as Fermi-LAT sources before the end of the mission, with the next most likely candidate 



being the M31 satellite galaxy M33 (Abdo et al. 2010c)). 



It is also insightful to compare the gamma-ray emission from local star forming galaxies to the 
emission from nearby globular clusters. The Fermi-LAT has now detected approximately ten glob- 



ular clusters in gamma-rays, with the emission most likely arising from milli-second pulsars (Abdo 



et al. 2010e). For globular clusters there is an interesting observed correlation between the density 



of the system and the gamma-ray luminosity, which may be reflective of the fact that milli-second 
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Figure 12 Gamma-ray luminosity as a function of star formation rate for galaxies detected by 



Fermi-LAT. Adapted from Ackermann et al. (2012b) 



pulsars are born in stellar captures in the dense clusters. In contrast, since the dSphs are much 
more extended than globular clusters (Figure [2]), the stellar capture rate in them is much lower, so 
the capture of stars and remnants into binary systems does not occur. Only one dSph, Sculptor, 
has a population of X-ray binaries associated with it, which likely formed from an entirely different 



physical mechanism (Maccarone et al., 2005). 



In order to precisely predict the gamma-ray signal from dSphs, we must return to understanding 
the dark matter distributions within them. It is informative to begin with some heuristic arguments 
that will be important in the context of these measurements. First, note that the dark matter halos 
of the dSphs likely do not extend beyond approximately a few kpc. This limit can be motivated 
by determining the Roche limit for a dSph for a minimum estimate of dark matter it contains, 
the known distance D to it, and the assumption that it is on a circular orbit. The Roche limit is 
specifically given by R t ~ D [Md sp h/M mw (Z))] 1/ ' 3 , where M mw (D) is the mass of the Milky Way out 
to the distance of the dSph. So in order to obtain a conservative lower estimate on Rt we must have 
a reliable estimate of the minimum mass of the dSph. This in turn depends on the largest radius 
out to which we find stars that are associated with the dSph, and understanding if these stars are 
tidally-bound to the system. For all of the classical satellites, it is certainly robust that within 



spherical radii of a few hundred parsecs the stars are faithfully tracing the local potential (Strigari 



et al., 2008) 



As discussed in Section |4j the modern kinematic data strongly constraints the integrated dark 
matter masses of the dSphs within their respective half-light radii of approximately a few hundred 
parsecs. This statement is weakly dependent on whether there is a central core or cusp in the 
dSph, so long as the log-slope is —d(logp)/dlogr < 1.5. The importance of this statement cannot 
be overestimated, because the constraints on the mass profile directly translate into constraints 



on the integral over the density-squared within the same region (Strigari et al. , 2008). Indeed, it 
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Figure 13 J- values for dSphs within a radius of 0.5° as a function of their Galactocentric distance. 
In this case NFW profiles are assumed for the dark matter density, though the results are very 
weakly dependent on the assumed central dark matter profile for dSphs with large data samples 



(compare to Figure 14). 



can be shown that in a manner similar to the calculation for the integrated mass in Section |4j the 
J-value is best constrained within an integrated physical radius that strongly correlates with the 



half-light radius (Walker et al. 2011). To better appreciate this, consider that the nearest classical 
satellites are at distances of approximately 70 — 80 kpc. For a dSph at this distance, the half-light 
radius corresponds to less than approximately one degree, which is about the angular resolution of 
the Fermi-LAT over a large energy range of interest. This is the region within which the integrated 
density and the integrated density-squared are the best constrained from the kinematic data sets. 
Thus the assumption of a core or a cusp for the density profile does not significantly affect the 
gamma-ray flux predictions for the Fermi-LAT. As discussed more below, however, for instruments 
with better angular resolution than the Fermi-LAT, the assumption of a core or the cusp is much 
more relevant. 

The theoretical developments outlined above have significantly improved the determinations of 



et al. 


2000, 


Tyler 


2002, 


Evans et al. 


2004 


and ] 


Martinez et al. 


(2009) have developed 



from stellar kinematical and photometric data using the likelihood in Equation 36 More recently 
groups have extended this analysis though in all cases the calculations are generally in good agree- 



ment (Charbonnier et al. 2011). 

For nearby dSphs that are most relevant for gamma-ray observations, the most updated de- 
terminations of the J-values are shown in Figure 13 Here an NFW profile is assumed for the 



dark matter density profile, as in Ackermann et al. (2011). However, as is shown in the proba 



bility density in Figure [14] the results are weakly dependent on whether a cored or cusped central 



density profile is assumed for the dark matter. Figure 13 clearly indicates which dSphs are the 
most interesting targets for indirect dark matter detection experiments. The two dSphs with the 
largest J-values, Segue 1 and Ursa Major II, are ultra-faint satellites with sparse samples of stars 



associated to them. Specifically, the J-values for Segue 1 and Ursa Major II in Figure 13 were 
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Figure 14 J-values for Draco within 0.1° (left histograms) and 0.5° (right histograms). The left 
panel assumes an NFW profile, while the right panel marginalizes over inner slopes of the dark 
matter density profile from a range of (core) to 1 (NFW cusp). Comparison of the peaks of the 
right-hand (blue) histograms shows that the J-value is weakly-dependent on the assumed central 
slope of the dark matter density profile within 0.5°. 



derived from samples of 66 and 20 stars, respectively. Though these dSphs have the largest mean 
flux, they also have the greatest uncertainty due to the small stellar samples. After Segue 1 and 
Ursa Major II, the dSphs with the next largest J-values are Ursa Minor and Draco, at 66 and 80 
kpc, respectively. These J-values are determined from samples of 150 and 320 stars, respectively, 
and as is seen uncertainties are much lower than the uncertainties on the J-values for Segue 1. 
In all cases, the one-sigma uncertainties are approximately log-normal; Figure 14 shows example 
probability densities from which the J-values in Figure [13] are derived. 

Using the J-values and their associated uncertainties, gamma-ray fluxes from each dSph are 
predicted using a model for the annihilation cross section and resulting gamma-ray spectrum. 
Analysis of two year Fermi-LAT data for all dSphs reveals no measured excess over the emission from 



diffuse Galactic backgrounds, extragalactic backgrounds, and nearby point sources (Ackermann 



et al. 2011 Geringer-Sameth and Koushiappas, 2012). With the null detection, a robust upper 



limit may be placed on the mass and the annihilation cross section. The procedure for determining 
the limit involves combining the posterior probability distributions of the J-value of each dSph 
with the likelihood for the gamma-ray data in the direction of each individual dSph (or a Region of 
Interest, ROI). In each ROI, the normalizations of the diffuse Galactic, diffuse extragalactic, and 
nearby points sources are separately fit to the data. While the J-values are distinct for each ROI, 
the parameters (cr ann ?;) and M^ m are shared between ROIs, thereby tying together different ROIs. 
The joint likelihood for the data, D, given N ROIs, is 



L(D\p w ,{p}) = l[L^ AT (D\p w ,p t 



'lit a. 



exp 



(log 10 (J) -log 10 (JQ) z 
2a? 



ln(10)J. 

where L^ AT is the poisson likelihood for a ROI, D is the binned gamma-ray data, p w are the pa- 
rameters that are independent of the ROI, in this case {aa, nn v) and m, and {p} are the parameters 
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Figure 15 Left: Limits on the annihilation cross section as a function of WIMP mass, assuming 
100% annihilation to bb. Upper curves are limits from individual dSphs, and the lower solid curve 
is the combined limit. Right: Limits on the annihilation cross section as a function of WIMP mass 



for several different channels. From the Fermi-LAT Collaboration (Ackermann et al. 2011). 



that depend on the ROI, in this case the background normalizations. The a l in this model are 
one sigma uncertainties on the J-values as derived from their posterior probability distributions. 



A profile likelihood method is used on Equation 90, which maximizes L for all model parameters, 
and from this profile likelihood 95% c.l. limits are defined. The above procedure for determining 
the upper limits on (caring) for a fixed mass was implemented by the LAT collaboration (Ack- 



ermann et al. 2011). Similar upper limits are obtained by using an "on-off" source analysis for 



each dSph, which involves comparing the expected source region to a nearby region in which the 
background is smooth. In this case the background components do not explicitly have to be fit 
out ( Geringer-Sameth and Koushiappas, 2011). Resulting upper limits as a function of mass are 



shown in Figure 15, for several different assumptions about the dominant annihilation channel for 
the dark matter. 

Analysis along the lines above has recently been implemented in a search for gamma-ray lines 



from dSphs (Geringer-Sameth and Koushiappas, 2012). This is probably the most robust indirect 



detection signal imaginable because it would be a gamma-ray line emanating from a region where 
the dark matter distribution is well-understood. Unfortunately the predicted signal is well-below 
that expected from standard supersymmetric models discussed in Section [6| There is no evidence 
yet for gamma-ray lines from dSphs; the upper limits on the annihilation cross section into gamma- 
ray lines are (o" ann u) ~ [10~ 27 , ~ 10 -24 ] cm 3 s _1 for masses of [10, 100] GeV, respectively. 

Similar to the case of Galactic center observations, ACTs have also studied dSphs in search of a 
signal from dark matter. Because of the high energy and thus photon-limited regime in which ACTs 



operate, the dSphs are analyzed via an on-off strategy described above. Figures 16 and 17 present 
the limits on (cr ann u) from MAGIC (|Aleksic et al] |2011[) and VERITAS (|Rico[ |2011[) observations 



of Segue 1. Figures 16 and 17 clearly show that ACTs are more sensitive to the higher dark 
matter mass regime than the Fermi-LAT analysis, though because of smaller exposures the limits 
are weaker by a couple of orders of magnitude as compared to the Fermi-LAT limits. Similar to 
the Fermi-LAT analysis, J-values were determined from the kinematic data within the appropriate 
angular regime under the assumption of an Einasto dark matter profile. The indicated limits are 
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Figure 16 Limits on the annihilation cross section as a function of WIMP mass, assuming 100% 
annihilation to bb and rf. From the MAGIC collaboration (Aleksic et al. 2011). 
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Figure 17 Limits on the annihilation cross section as a function of WIMP mass, assuming 100% 
annihilation to bb and rf, and W + W~ . From the Veritas Collaboration (Aliu et al. 2012). 



at the 95% c.l. Even though the limits from the ACTs are two to three orders of magnitude larger 
than the characteristic thermal relic scale, they do probe some of the parameter space for a dark 
matter explanation to recent cosmic ray data (see below). 

Forthcoming larger scale ACTs with angular resolution < 0.1° will be more sensitive to the dark 
matter distribution in the center of the dSphs. This is quantitatively shown in Figure [l4j which 
compares the integrated J- values within 0.1° and 0.5° for cored and cusped dark matter profiles. 
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(See also Charbonnier et al. (2011)). The uncertainties in this figure are derived only from the 
kinematic data; before the launch of next generation ACTs it will also be important to improve 
the photometric data in the center of the dSphs. As was shown in Section [4] there is a significant 
degeneracy between the central density of the stellar profile and the central density of the dark 
matter profile. 

In summary, indirect dark matter searches using dSphs are now at an advanced stage, and have 
delivered the most robust results on the dark matter mass and cross section. To the extent that 
is possible, they represent particle dark matter experiments in the sky. Gamma-ray observations 
of dSphs have been the first to achieve reliable sensitivity to thermal relic dark matter over a 
theoretically-motivated mass regime of 10 — 25 GeV, and they will continue to push down limits 
to the thermal relic scale at larger masses. The sensitivity will continue to increase because of 
improved understanding of the systematics in the gamma-ray data, and also because new Milky 
Way satellites are sure to be detected in future surveys. Also in future analyses it will be important 
to model the satellites as extended sources, improving upon the point source modeling that has 
been used to derive the limits described above. 

7.4- Diffuse Galactic and extragalactic searches 

Residual annihilation of dark matter particles throughout the Milky Way, and throughout 
external galaxies over all redshifts, provides an interesting signature that is complementary to the 
Galactic center and dSph searches. This is a truly diffuse signal that is spread out over many 
resolution elements of the detector, so by its nature it is more difficult to extract than searches 
for point sources. With its wide field-of-view, the Fermi-LAT is uniquely sensitive to a diffuse 
signal along these lines. Extracting it requires a detailed model for both the spectral and spatial 
distribution of the signal as well as the spectral and spatial resolution of the gamma-rays produced 
from cosmic ray sources. 

For diffuse searches, there are three relevant mechanisms of gamma-ray production. First, there 
is the standard prompt production that arises from hadronization, fragmentation, and subsequent 
ttq decay, which as discussed above is strongly constrained for bb and rf channels from observations 
of dSphs. Second, diffuse radiation can also be used to search for gamma-ray lines that are produced 
directly through a 77 or 7Z pair. Third, the inverse Compton component produced by the up- 
scattering of starlight by high energy electrons and positrons produced by the annihilating dark 
matter. Note that this third process is not relevant in the dSph searches because of the low density 
of starlight in the outer MW halo. 

The challenge in extracting the diffuse dark matter signal is to characterize the diffuse gamma- 
ray background from cosmic ray sources. The observed spectrum of the diffuse background radiation 
falls off like a power law over a wide energy regime, which is much different than the gamma-ray 
spectrum from dark matter. Several approaches for separating out the latter component have been 
considered in the literature. Using standard models for the Milky Way dark matter distribution 
and simply ascribing all of the background to other processes, the most conservative constraints on 



the annihilation cross section at 10 GeV are a few times 10 24 cm 3 s 1 (iPapucci and Strumia 



2010 



Cirelli et al. 2010). Improving on this simple limit requires the development of more sophisticated 



theoretical models. Baxter and Dodelson (2011) discuss an interesting technique to set a limit on 
the annihilation cross section that exploits the expected difference in the distributions between the 
signal and backgrounds in a ring of pixels around the Galactic center. They derive conservative 
upper limits of approximately 10~ 25 cm 3 s _1 on the dark matter mass of approximately 10 GeV for 
annihilation into bb and rf channels. Though this technique still is subject to uncertainties on the 
dark matter distribution it does provide an interesting method for reducing the Galactic emission 



81 



from cosmic rays that peaks in the disk, without the need to understand details of the latter signal 
other than its spatial distribution. 



Ackermann et al. (2012e) have recently performed a search for dark matter annihilation in 
the diffuse halo, modeling both the distribution of dark matter, and the properties of the cosmic 
ray source populations and their propagation through the Galaxy. A given model of the Galaxy is 
defined by the spectral index of the injection cosmic-ray energy spectrum, the scale-height of the 
diffusive cosmic ray halo, and the dust to HI ratio. After fixing a set of "non-linear" parameters to 
define a model for the Galaxy, over twenty "linear" parameters must still be determined that de- 
scribe the cosmic ray source distribution, the normalization of the isotropic gamma-ray background 
in each bin, and the amplitude of the dark matter signal. Maximizing these linear model param- 
eters allows for determination of the upper limits (or best-fitting) amplitude for the dark matter 
signal. The procedure is then repeated for a new set of non-linear parameters, and from a variety 
of these repeats different sets of constraints are placed on the linear parameters. Unfortunately, 
the large number of model parameters makes it difficult to marginalize over all of the parameters 
that describe the cosmic ray sources and their propagation. Given the complexity of the modeling, 
it is understandable that no conclusive dark matter signal can be extracted with the present data; 
at ~ 10 GeV Ackermann et al. (2012e) quote upper limits in the range cm s for 

bb and tt channels, where the uncertainty is driven by the assumed local dark matter density. 

The Fermi-LAT collaboration has reported two searches for gamma-ray lines throughout a large 
volume of the Galactic halo (Abdo et al. 2010a; Ackermann et al. 2012a). In the latest results, 



which are based upon two years of data, the annihilation cross section upper limits vary from 
approximately 10~ 28 cm 3 s -1 at 7 GeV to approximately 10 -26 cm 3 s _1 at 200 GeV. There are 
still factors of several uncertainties in these limits due to the uncertain dark matter distribution in 



the Milky Way (Figure [IT]) . 

Subtraction of the diffuse gamma-ray emission from Galactic sources leaves a residual, isotropic 
background. This background is believed to arise from a variety of unresolved sources, including 
active galactic nuclei, starburst galaxies, and shocks created during the formation of large scale 
structure. Over the energy range 200 MeV - 100 GeV, the spectrum of the isotropic emission is a 
power law that scales scales as E~ 2A (Abdo et al. 2010d). Though the above sources are expected 



to contribute to this emission, there is no well-accepted model for the relative contributions of each 



(Fields et al. 


2010 


Inoue 


2011 



energy spectrum is a pure power law, it is clear that dark matter annihilation does not explain it 
over a large range of energies. The pertinent question is then whether dark matter can explain the 
background over a smaller energy window. 

The theoretical mass function of dark matter halos is now well-understood from N-body simu- 
lations, and from this mass function it is possible to predict the contribution to the extragalactic 
background from dark matter annihilation. Because the emission is from halos and subhalos over 
the entire mass range over all redshifts, the uncertainty on the boost factor from halo substructure is 
an important systematic in determining the background component from dark matter. Further, the 
uncertainty in the contributions to the background from non-dark matter sources preclude robust 
statements from being made on the mass and cross section at this time. However, it is possible to 



rule out some exotic models of the annihilation cross section using the isotropic background (Abdo 



et al., 2010b) 



7.5. Galactic substructure 

As discussed in Section [5j N-body simulations of Galactic halos predict that the dark matter 
mass distribution in the halo is not smooth, but rather approximately 10 — 50% of the mass is 
bound up in the form of subhalos with a mass function that may extend down to Earth masses 
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Figure 18 Distance versus mass for subhalos in the Aquarius A-2 simulation (Springel et al. 2008a). 



Above the diagonal line, the subhalos have a J-value that is less than Draco, so they would not be 
visible in the Fermi-LAT, given that no signal is yet observed from Draco. Adapted from Ackermann 



et al. (2012d) 



or even below. Gamma-rays produced from dark matter annihilation in these subhalos may also 
be detectable, either from individual subhalo sources or from the cumulative effect from all such 
subhalos distributed throughout the Galaxy. 

Recall that in the LCDM model, the dSphs reside in a subset of the dark matter subhalos that 
are resolved by N-body simulations. When performing searches for dark matter annihilation from 
subhalos, the first question to ask is what fraction of these subhalos would be detectable at all given 
that no signal is yet observed from the dSphs, which have well-determined dark matter masses. 



Following Ackermann et al. (2012d), for the purpose of illustration, compare the potential signal 



from dark subhalos to the dSph Draco, which is at a distance of 80 kpc, has an integrated mass 
within 300 pc of ~ 10 Mq (Strigari et al. 2008), and a total halo mass in the range ~ 10 8 — 10 9 



Mq. Figure 18 shows the mass versus distance for subhalos from the Aquarius simulation, and 
on the diagonal line subhalos have the same J-value as Draco, assuming that the mass of Draco 
is ~ 10 8 Mq. Subhalos above the diagonal line would not be visible, while subhalos below the 
diagonal line do have higher J-values than Draco and may be detectable. It is interesting to note 
that a large majority of the dark subhalos that are detectable reside on the most massive end of 
the subhalo mass function. This behavior still holds even if the subhalo mass function continues 
as a power law dN/dM ~ M -1 ' 9 to below the resolution limit of modern simulations. In addition 
to these massive dark subhalos, if the mass function does continue down to approximately Earth 
mass scales, it still may also be possible that a dark subhalo is very local to the Sun, < 1 kpc, with 
emission detectable by the Fermi-LAT. 

Recall that the two-year Fermi-LAT point source catalog contains 1873 gamma-ray sources, 
over 30% of which are unidentified at other wavelengths. So the first place to start the search 
for subhalos is to examine these unidentified sources. Based on their energy spectrum alone, it is 



possible to show that a fraction of these can be fit to dark matter models (Buckley and Hooper 
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2010), however further tests are required for a more robust identification. 



In order to robustly identify a gamma-ray sources as a subhalo, several criteria must be met. 
First, the source must have a characteristic spatial extension. It is most likely that either a massive 
subhalo in the outer part of the Galactic halo or a much lower mass subhalo that is very nearby will 
have a spatial extension that is resolvable by the Fermi-LAT. Second, if the source is a subhalo, its 
energy spectrum will deviate from a pure power law model. If the continuum gamma-ray spectrum 
is similar to a bb spectrum, simulations find that subhalo sources can be reliably distinguished 
from pure power law energy spectra (Ackermann et al. 2012d). However, the energy spectra of 



pulsars typically have an exponential cut-off that mimics a bb spectrum, therefore based on the 
energy spectrum alone it is easy to confuse a potential subhalo source with a pulsar. Combining 
the extension and energy spectrum criteria with information on (lack of) variability, and cutting 
out sources near the Galactic plane where confusion is high, in the one year data no subhalo source 



could be uniquely identified (Ackermann et al. 2012d). For assumptions on the distribution of 
subhalos, this lack of detection corresponds to an upper limit on (cr ann f) of approximately 10~ 24 
cm 3 s" 1 at 100 GeV. 

The above analysis is based on the identification of individual subhalos as gamma-ray sources; 
since LCDM models extrapolate to predict up to 10 16 subhalos in the Milky Way, it is important to 
ask whether the diffuse emission from these subhalos may be detectable as features in the gamma-ray 
background. Due to the spatial distribution and energy spectrum of the subhalos in the Galaxy, 
the probability distribution for the gamma-ray flux has a longer power law tail than a poisson 

2009). Exploitation of the difference between this signal and background 



distribution (Lee et al 



distribution may lead to a robust identification of the dark mater mass and cross section (Baxter 



et al. 2010). It is most likely that signals such as these could be extracted if a smooth emission 



component from the Galactic center or dSphs were discovered first. 



7.6. Galaxy clusters 

Clusters of galaxies are also unique and distant targets from which a dark matter signal may 
be extracted. Gamma-ray photons from continuum emission, line emission, and inverse Compton 
emission are all produced through dark matter annihilation in clusters. Since clusters in many 
instances have well-determined mass distributions, in a manner that parallels the dSph analysis the 
gamma-ray emission from dark matter annihilation can be predicted starting from astrophysical 
data. 

The mass of a cluster is typically defined within a fixed overdensity, A, with respect to the 
critical density, p c = 3H 2 (z)/8irG, as Ma = M(< r&) and MA/(4/3)7rrA = A x p c . The radius 
that encloses the given overdensity is r&. The most robust method to derive the mass of a cluster 
involves measuring the temperature profile from X-ray spectra, and the electron gas density profile 
from the X-ray luminosity. For nearby clusters with large values of A, masses can then be robustly 
derived under the assumption of hydrostatic equilibrium, 



M(r) 



kT(r)r 

fMUpG 



dlnp(r) d In T(r) 
dlnr d In r 



(90) 



where T{r) is the temperature profile, and [i is the mean mass per particle in units of the proton 
mass. Numerical simulations find that the systematic uncertainty to due non-thermal pressure 
introduces a ~ 10% systematic in the mass determination (Nagai et al. , 2007). For clusters without 



temperature or density profiles, the total mass of the cluster is derived from several observational 
proxies, which have been calibrated using both observations of low-redshisft clusters and numerical 
simulations. Three standard observational proxies are the average temperature, the mass of the 
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hot gas, and the product of these two, which is the total thermal energy (Vikhlinin et al. 2009). 



The masses of the majority of the clusters are measured via these scaling relations; there are tens of 



clusters that have measured temperature and density profiles making it possible to use Equation 90 



Predictions for the gamma-ray signal from dark matter annihilation in clusters typically as- 
sume an NFW profile for the dark matter distribution. Recall that in the LCDM framework, the 
properties of an NFW dark matter halo are entirely specified once the total mass is determined. 
From the total mass estimate, Ma, the concentration, c v i r , is determined from Equation 51 From 
these quantities the virial radius, R v i r , is determined by standard formulae, e.g. Mailer and Bullock 
(2004). The NFW scale radius is then r s = R v i r /c v i r , and the scale density is determined from 



Ps ipvir) 



47rr? 



ln(l + c v 



1 + c v 



-i 



(91) 



Similar to the dSphs, the most important clusters for indirect dark matter detection are those 
that are the appropriate combination of the nearest and the most massive. Based solely on their 
derived mass distributions from the formalism outlined above and their distances, several recent 
studies have come to the agreement that the Fornax cluster is the brightest source of gamma- 



rays from dark matter annihilation (Ackermann et al. 2010 Ando and Nagai, 2012; Pinzke et al. 



2011, Nezri et al. 2012). Coma and Virgo are also expected to be amongst the brightest gamma- 



ray sources. Though formalism above for determining cluster masses is well-motivated from a 
theoretical point of view, systematic uncertainties will be introduced if the mass profile deviates 
from an NFW model. This is particularly relevant for the gamma-ray analysis because a few of the 
most nearby clusters are predicted to be the brightest sources of gamma-rays. 

Unlike the case of dSphs, there is expected to be a significant increase in the gamma-ray signal 
from clusters due to dark matter substructure. This is pleasing from a detection perspective because 
it provides a well-motivated reason for clusters to have brighter emission than is deduced from the 
smooth dark matter component. However, in the case of a null-detection, it introduces a theoretical 
systematic that makes it difficult to interpret the meaning of an upper limit on the annihilation 
cross section. 

Numerical simulations provide the most reliable method for determining the gamma-ray emis- 
sion from subhalos in clusters. From a sample of nine cluster-mass dark matter halos, which have 
individual particle masses of ~ 10 6 M & and identify subhalos down to a mass scale of ~ 10 7 M , 
Gao et al. (2012) directly determine the overall boost factor and the surface density profile of the 
substructure component for clusters as a function of the cluster mass M2oo- Due to the resolution 
limit of their simulations, only the gamma-ray luminosity near the cluster virial radius increases 
due to subhalos. Starting from a mass scale above which subhalos are resolved and extrapolating 

" 6 Mr: 



down to 10 



'0, 



Gao et al. (2012) determine a boost factor of 
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7rm(17) rloo + 16^ 2 ' 

where R is the projected radius and C is the luminosity of the smooth component of the cluster. 
Though it relies on an extrapolation, the increase in the luminosity over the smooth component in 
clusters is much greater than is expected in dSphs. 

Also unlike the case of dSphs, in clusters there is expected to be a significant flux of gamma- 
rays due to cosmic ray processes. For clusters of interest, the emission separates into component 

85 



functions that depend on radius and energy, so that the source function decomposes as 



dN 



dtdVdE 



A(r)s(E). 



(94) 



The spatial factor, A(r), is proportional to the square of the gas density in the cluster, and the en- 
ergy component, s(E), varies as dn/dE ~ E~ 2 ' 3 for energies up to approximately 100 GeV. Because 
the gamma-ray luminosity in clusters traces the gas density, it is more centrally-concentrated than 
the dark matter-induced signal, which is more extended because of the emission from subhalos. 
Theoretical calculations have shown that both the cosmic ray and the dark matter annihilation 



induced gamma-rays may be detectable by the Fermi-LAT (Pinzke et al. 2011 Nezri et al. 2012) 



However, no gamma-ray emission has been conclusively associated with clusters in the Fermi-LAT 



one year data, assuming that they are point sources (Ackermann et al. 2010). Modeling clusters 



as extended sources, no gamma-ray emission was yet found in two year data (Han et al. 2012). 



At a dark matter mass of 10 GeV, the limits on the annihilation cross section deduced from these 
observations vary anywhere from 10~ 24 cm 3 s~ 1 to below 10 -26 cm 3 s , depending on what is 
assumed for the emission from subhalos. 

HESS has made a pointed observation of the Fornax cluster, obtaining a total exposure time 
of 14.5 hours (Abramowski et al. 2012). From a null detection, it has placed limits on a ann v 



for different assumptions of the dominant annihilation channel and the slope of the central dark 
matter profile. The limits obtained are very sensitive to both of these assumptions, in the latter 
case because the angular resolution of HESS probes the central region of the dark matter profile of 
the cluster. Depending on the assumption for the amount of substructure in the cluster, the limits 
on (<7 ann t> ) are in the range 10 -21 cm 3 s _1 to 10" 



-23 cm 3 g 



7. 7. Summary of gamma-ray results 

Table [6] provides a summary of modern sensitivity limits for dark matter of mass 10 and 100 
GeV from the various gamma-ray searches described above. 

Because no clear dark matter signal has been detected yet in gamma-ray data, the above 
analysis has primarily focused on methodology for deriving limits on the annihilation cross section. 
Recalling the discussion in the beginning of this section it is always important to keep in mind that 
the nature of the signal we are searching for is not precisely predicted, so it is necessary to pay 
careful attention to details and residuals in data sets where a signal may be lurking. It must also be 
kept in mind that our knowledge of astrophysical systems is incomplete, and also our understanding 
of experimental and detector related issues is often incomplete. It is often tempting to associate 
our lack of understanding of conventional astrophysical processes or detector systematics with a 
dark matter signal. 

Though as discussed above the published Fermi-LAT results report null searches for gamma-ray 
lines, recently a feature has come to light in the data near the Galactic center where an excess over 



a pure power law model for the background is found (Bringmann et al. 2012 Weniger, 2012) 



Even though the hint of a line-like signal is intriguing, there are certainly reasons to proceed with 
caution. First, the signal does not directly coincide with the Galactic center, so a mechanism must 
be invoked to displace the center of the dark matter halo from that of the baryons. Second, if 
the signal is in fact a result of dark matter, it should be much more spatially extended above 
and below the Galactic plane. Third, the power law parameterization of the background may not 
be appropriate for the small spatial regions that are being examined. Fourth, from a theoretical 
perspective, the line emission should generically be accompanied by continuum photons produced 
via hadronization and fragmentation. Since no continuum emission is found above the power law 



background (Cohen et al. 2012; Cholis et al. 2012), the ratio of the continuum to line signal must 
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Table 6. Summary of limits on the annihilation cross section times the relative velocity, cr ann z;, in 
units cm 3 s" 1 for dark matter of mass 10 and 100 GeV for the different sources. For mass of 10 
GeV, the limits are given for 100% annihilation into bb and r + r~. For a mass of 100 GeV, the 
limits are given for 100% annihilation into bb, t + t~ and W + W~ . 
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Ackermann et al. (2011); Geringer-Sameth and Koushiappas (2012) 



b Results for NFW profile. Subject to the uncertainty in the local dark matter density. 
From Ackermann et al. (2012e) 



Abdo et al. (2010b) 



d Assumes just the emission from the smooth component and not subhalos. From Han 



et al. (2012); Ando and Nagai (2012) 



be much less than is generally found in particle dark matter models. A larger sample of data and 
cleaned data sets should be able to shed more light on this feature over the next few years. 

7.8. Neutrino searches 

Large-scale neutrino experiments, in particular Ice Cube and Super-Kamiokande, are now gain- 
ing sensitivity to study dark matter through indirect detection. Both of these experiments detect 
"up-going" muon neutrinos that interact with nucleons in the Earth to produce muons that make 
Cerenkov light in the respective detectors. This mode of observation is used primarily to reduce 
the background from atmospheric muon neutrinos, which produce a background from "down-going" 
muons above the horizon. The lower energy threshold for Ice Cube observations is 10 GeV, while 
Super-Kamiokande extends to a lower energy threshold allowing it to study lower mass dark matter. 
The angular resolution for up-going muon events is approximately one degree for these experiments. 
This section reports results from searches for neutrinos from the various sources discussed above 
in the context of gamma-ray searches. Searches for neutrinos from dark matter that gets captured 
by scattering off of nuclei in the Sun will be discussed in the section on direct dark matter lim- 
its (Section [9]); these observations are most sensitive to the spin-dependent WIMP-nucleon cross 
section. 

Ice Cube has now reported limits on the annihilation cross section into several channels through 
observations of the diffuse Galactic halo with an exposure of 276 days and a 22-string configura- 



tion ( Abbasi et al. , 2011 ). For dark matter mass greater than 20 GeV, the best-constrained channel 
is the direct annihilation of dark matter to neutrinos, for which the limits are cr a nnW ^ 10 -22 cm 3 s _1 . 
For a bb annihilation channel, for a dark matter mass of 30 GeV the upper limit is o~ ann v ~ 10~ 19 
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-1Q S — 1 

a cur s . 



For 



, while for a dark matter mass of 10 TeV the upper limit is o" ann t> 
the other channels the limits are between the bb and direct neutrino limits. 

Ice Cube has very recently presented limits on the annihilation cross section into neutrinos 



from direct observations of the Galactic center (Abbasi et al. 2012b). For direct annihilation into 
neutrinos, the limit on the annihilation cross section at approximately 200 GeV is ~ 10~ 23 cm 3 
s _1 , which improves upon the aforementioned diffuse limits by about an order of magnitude. For 
the remaining channels, the limits are typically improved in comparison to the diffuse limits for 
WIMP masses less than a few hundred GeV. 



7.9. Antimatter searches 

Positrons and anti-protons produced in dark matter annihilation will propagate to the Earth, 
and their detection is a complementary probe of dark matter in the Galaxy. There are now de- 
tailed theoretical treatments that analyze the propagation of positrons and anti-protons from their 
production sites, either from a dark matter source or astrophysical source such as pulsars, through 
the Galaxy. For non-dark matter sources, secondary production of positrons and anti-protons 
is probably the best established mechanism for describing their propagation. In this model cos- 
mic ray nuclei interact inelastically with interstellar gas, producing charged pions that decay to 



positrons (Moskalenko and Strong 1998). This process typically results in a positron fraction that 



declines with energy. 

Though not specifically designed for indirect dark matter detection, within the past few years 
it has been realized that antimatter experiments do have a unique position in Galactic searches 
for dark matter. The PAMELA experiment has found an interesting rise in the positron fraction 



at energies up to 100 GeV (Adriani et al., 2009b), extending the previous results from the HEAT 



experiment (Barwick et al. 1995). This result has been corroborated by the measurements of 



the positron spectrum from the Fermi-LAT collaboration (Abdo et al. 2009b Ackermann et al 



2012c) and HESS (Aharonian et al., 2009). Measurements of the anti-proton flux by PAMELA 



on the other hand are in agreement with the predictions for local astrophysical sources (Adriani 



et al. 2009a). Though the rise of the positron spectrum may be connected to interesting local 



astrophysical sources, it is unlikely to be due to dark matter because the required cross section is 
approximately two orders of magnitude larger than the thermal relic scale. Forthcoming results 



from the AMS-02 spectrometer are expected to shed more light on the PAMELA results (e.g. Hooper 



and Xue (2012)). 



7.10. Decaying dark matter and X-ray searches 

Up to this point the analysis in this section has focused on stable dark matter particles that 
annihilate into detectable standard model final states, restricting to energies in the gamma-ray 
band and above. A variety of models also predict that the dark matter is unstable, and a small 
fraction of the particles are decaying into detectable final states at the present epoch. For these 
searches the energy of the radiation produced can extend down to the X-ray energy band, and even 
below. 

Sterile neutrinos at the mass scale of approximately a few keV that are stable over the lifetime 



of the universe provide good candidates for dark matter (Abazajian et al. 2012). Because the 



sterile neutrino mixes with standard model neutrinos, it radiatively decays into a photon and a 
lighter neutrino. For a majorana neutrino, this decay rate is given by 



r 



1.4 x 10" 30 s _1 



sin 2 26 



10 



-7 



m. 



IkeV 



(95) 



SS 



where m s is the mass of the sterile neutrino and sin 2 26 is the mixing angle. This mixing angle is 
constrained to be < 1CP 6 to avoid overproduction in the early universe. 

Beyond approximately 20 keV, the null detection of M31 from XMM-Newton observations 



provides the best limit on the mass and mixing angle of sterile neutrinos ( Watson et al. 2006 , 2012 ) . 



Sterile neutrino masses greater than about 30 keV can be ruled out from the non-observation of 



lines from the INTEGRAL satellite (Yuksel et al. , 2008). Null reports have also been reported from 



searches for an X-ray signal from decaying dark matter from Galactic satellites (Loewenstein and 



Kusenko, 2012). 



7.11. Summary of results and interpretations 

As has been outlined through the course of this section, over the past several years indirect dark 
matter searches have made an extraordinary amount of progress. We have gone from wondering 
what contributions the Fermi-LAT could make to dark matter studies, to for the first time probing 
thermal relic scale cross sections. It is now appropriate to take a step back and consider some 
examples of what these results tell us about some theoretical models. 

because of the lack of velocity suppres- 



In terms of the models presented in Equations 72 ■ 80 



sion in the annihilation cross section the indirect detection results are most sensitive to fermions 
with psuedoscalar, vector, and tensor interactions, and scalar dark matter with scalar, scalar- 
pseudoscalar, and vector- axial vector interactions. In order to see how the results presented here 
constrain these dark matter models, it is informative to consider a specific example. For scalar 
dark matter consider the case of scalar interactions, with an annihilation cross section given by 
Equation 72 For universal couplings to fermions, the relic abundance is set by annihilation into 



all of the fermions. The relic abundance serves to set the couplings F s . For the sake of clarity, say 
that the relic abundance is set by the annihilation into bb, for which case the coupling is then fixed 
to F s ~ 3 x 10~ 4 GeV -1 . Including all fermionic annihilation channels decreases this by a factor of 
a few. It is clear that the Fermi-LAT dSph results now probe scalar couplings at this level in the 
mass range 10 — 25 GeV. 

A second interesting example is fermionic dark matter particles with pseudoscalar couplings. 
This is particularly interesting because, as discussed in detail in the following two sections, these 
are most probably unconstrained by direct detection experiments (in the absence of light mediating 
particles), so indirect methods are the only ones that can probe dark matter with dominant pseu- 
doscalar interactions. In contrast to the scalar dark matter case discussed above, the annihilation 
cross section in this case depends on the square of the dark matter mass (Equation 77). From 
the same arguments as above, for full annihilation into bb, the coupling that the Fermi-LAT dSph 
results are probing in the mass range 10 — 25 GeV is ~ 5 x 10~ 6 GeV -2 . This analysis clearly indi- 
cates that there is a complementarity between indirect and direct dark matter searches in probing 
different dark matter models. 
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8. Direct Detection: Theoretical Developments 



Goodman and Witten (1985) first suggested that if weakly interacting massive particles are 



the dominant component of mass of the Galactic halo, they could be detected by interactions with 
nuclei in low background underground detectors. Though the guiding theoretical framework is very 
much similar to what was originally discussed in the Goodman and Witten ( 1985 ) paper, in the now 
almost thirty years subsequent to this work, a tremendous body of literature has been established 
that has refined and expanded upon direct dark matter detection theory. The theoretical work has 
led to a steady understanding that the theory of direct dark matter detection is very rich, both 
from the perspective of the dependence on forefront topics in modern Galactic astrophysics and on 
the dependence on particle dark matter models. 

Since the initial papers on direct detection, numerical simulations have dramatically-improved, 
and it has been realized that they play an important role in advancing our understanding of the 
dark matter distribution in galaxies like our Milky Way. Although they are still nowhere near the 
resolution limit that is needed to study the detailed phase space structure of the dark matter on 
scales that are relevant for direct detection, they do allow us to extend beyond the most naive 
models and provide clues as to how different phase space models affect direct detection. From the 
particle physics perspective, the flurry of experimental results, and hints of possible signals, has 
motivated a wide variety of theoretical models designed to accommodate them. Even if in the 
long run the experimental hints are unable to withstand further scrutiny, these models have led 
researchers down interesting theoretical paths that likely would not have been taken otherwise. 
It is certainly true that, in order to eventually extract necessary information from experimental 
data, we must have an in depth understanding of how to separate the dependence of the signal on 
Galactic model parameters and those parameters related to the nature of the interactions of the 
dark matter particle. 

There are several excellent reviews that lay out both the astrophysical and the particle physics 



aspects of the theoretical framework for direct detection (Jungman et al. 1996, Lewin and Smith 



1996). This section builds upon these reviews by highlighting the important theoretical progress 



that has occurred during the past several years in this field. This section begins by reviewing the 
ingredients that go into determining the flux and the detection rate of WIMPs at direct detection 
experiments. Following this is a discussion of basic models for the WIMP-nucleus elastic scattering 
cross section; modern experimental constraints on these models will be discussed in Section |9j The 
focus then moves to a discussion of how the detected rate of WIMPs depends on the Galactic 
halo model, accounting for the smooth dark matter distribution, the distribution of dark matter 
substructure, and modern estimates for the WIMP velocity distribution. Finally, recent theoretical 
work is discussed that attempts to understand how the properties of the WIMP can be extracted 
from experimental data sets by rigorous modeling of the dark matter distribution in the Galaxy. 



8.1. WIMP detection methods 

The local flux of WIMPs is given by their number density in the Solar neighborhood times their 
velocity, nv . For target particles in an underground detector, the total cross sectional area seen 
by the WIMP is the number of targets, Nt, times the cross section per target, a. The WIMP 
interaction rate is then R = Ntcuv. 

The first part of the discussion in this section will take the observable signature of a WIMP 
to be the energy deposited to a nucleus in an underground detector, E r . Define the local velocity 
distribution of WIMPs as f(v,VE), where the notation is now so that v represents the WIMP 
velocity relative to the Earth and ve is the velocity of the Earth relative to the Galactic rest frame. 
The recoil momentum of the nucleus, in the lab frame, is q = 2m r vcos8, where 9 is the angle 
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between the incoming WIMP and the scattered nucleus, rriN is the mass of the nucleus, is the 
WIMP mass, and m r = m^Mdml '(tun + M^m) is the WIMP-nucleus reduced mass. The recoiling 
energy of the nucleus is then E r = q 2 /(2m^) = 2m 2 i> 2 cos 2 O/tun, and the minimum speed to 
scatter a nucleus to a recoil energy E r is 



2m 2 



(96) 



Note that Equation [96] is strictly valid only for elastic scattering interactions; for inelastic scatter- 
ings, a term & j \/2mjqE r is added to Equation 96 (Tucker-Smith and Weiner, 2001), with 5 being 



the mass splitting for particles in the spectrum. 

The velocity of a WIMP in the Galaxy is of order v ~ 10 _3 c. For WIMP masses and nuclei 
that we are considering, the reduced mass is typically well approximated by m r ~ 100 GeV. For 
the above velocity and m r , the momentum transfer is q ~ 100 MeV. The energy imparted to the 
nucleon is then E r = \q\ 2 /2m N ~ (100 MeV) 2 /(2 A GeV) = 10 4 keV/(2^), where A is the mass 
number of the nucleus. Since we will be focusing on nuclei with A ~ 100, the recoil energy of the 
nucleus is in the range of tens to hundreds of keV. 

With these definitions, the differential event rate as a function of nuclear recoil energy is then 
given as 



dR 



N T 



Pe 
M dm 



da(E r ,v) 
dE r 



vf{v, v E )d ,v, 



(97) 

p(Rq) is the local mass 



where v esc is the escape velocity in the Galactic reference frame, and p Q 
density of WIMPs. 

The velocity of the Earth with respect to the Galactic rest frame, ve, is a sum of three separate 
components: the velocity of the LSR, the peculiar motion of the Sun with respect to the LSR, and 
the orbital motion of the Earth around the Sun. To calculate each of these components, consider a 
coordinate system so that U points in the direction of the Galactic center, V points in the direction 
of Galactic rotation, and W towards the Galactic north pole. In these coordinates, the LSR velocity 
is vlsr — (0,220,0) km/s. The most recent determination of the mean velocity of the Sun with 
respect to the LSR gives v e = (U, V, W) = (11.1, 12.24, 7.25) km/s flSchonrich etaLj |201~0[ ) ■ 

Because of the motion of the Earth about the Sun, the flux of WIMPs in the Earth rest frame 



varies periodically during the course of the year (jFreese et al.[ 1988[ 2012|). This effect manifests 

Including both the 



itself in a time dependence in ve in the velocity distribution in Equation [97 
motion of the Sun through the Galaxy and the Earth motion about the Sun, 

v E (t) = v LS R + ^0 + ^0 [cosuj(t - ti)ei + sinuj(t - ii)e 2 ] • 



(98) 



where V® = 29.8 km/s is the orbital speed of the Earth around the Sun, u = 2-7r/year, e*x is the 
direction of the velocity of the Earth at a time t, and e% is the direction of the velocity of the Earth 
at a time t\ + 0.25 yr. 

The amplitude of the annual modulation signal is typically approximately 1-10% of the averaged 
total event rate. It is a function of v m i n , and depends strongly on the anisotropy of the velocity 
distribution. Due to the strong dependence on v m i n , from a precise determination of this phase 



reversal of the signal it may be possible to strongly constrain the WIMP mass ( Lewis and Freese 



2004) 



The above analysis has focused on determining the event rate as a function of energy deposited 
in a detector from a WIMP-nucleus scattering event. It has ignored the direction of the recoiling 
nucleus, by integrating over all possible recoil directions. Understanding, and measuring, the 
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distribution of nuclear recoil directions is also important, because it contains information both on 
the direction of the Earth through the Galactic halo and on the velocity distribution. Because of 
the former dependence on the direction of the Earth motion, it provides a clear confirmation that a 
signal determined from the energy spectrum is due to WIMPs and not from a various background 
component. Because of the latter dependence on the velocity distribution, the directional signal 
may contain information on the phase space properties of the dark matter halo (note that neither 
the outgoing WIMP momentum or its direction are detectable after a scattering event). 

The simplest derivation of the directionally-dependent event rate spectrum neglects both the 
motion of the Sun with respect to the LSR, and the motion of the Earth around the Sun. Define 
d£l = d(f)d(cos 9) as the solid angle associated with the direction of the recoiling nucleus, where as 
above 9 is the angle between the direction of the incoming WIMP and the scattered nucleus, and 
vr as the direction of the recoiling nucleus. From the dark matter velocity distribution f(v), the 



recoil rate as a function of both the energy and the direction is ( Gondolo , 2002 ) 



dfi 1 n f v esc fJ n 

" : N Ti^ d 3 vv 2 —f(v)5[(v-v Q )-v R -v mm (E r )}. (99) 



dE r dQ R 2vr M dm J Vmm dE, 

Here the explicit dependence of the cross section on the recoil energy and incoming velocity is 
indicated; simplified models for this cross section will be considered in more detail below. Note 
that for an isotropic velocity distribution, the event rate will be a maximum in the direction opposite 
to the motion of the LSR. 

8.2. Scattering cross section models 

The WIMP-nucleon scattering cross section fundamentally derives from the WIMP interactions 
with quarks. A theoretical model predicts specific interactions, which in general may reduce to 
functions of the WIMP spin, nuclear spin, the WIMP velocity, and the momentum transfer to the 
nucleus. Given the scale of the velocity and momentum transfer for WIMP-nucleon scattering, 
interactions that are not suppressed by powers of them will typically be the most relevant for 
direct detection. Spin-independent (SI) interactions are defined as those that do not depend on the 
nuclear spin, while spin-dependent (SD) interactions do depend on the nuclear spin. 

It is standard to define the total cross section in terms of the zero momentum transfer cross 
section, uq, as 



a 



7 2 

imax 



J dq 2 F 2 {q) (100) 



where qlL ax = 4m^w 2 is the maximum value of the momentum transfer. This assumes that the 



scattering is described by a contact interaction. The form factor is (Lewin and Smith, 1996), 

2 9 [sm(qRn) - qR n cos{qR n )] 2 

F iq) = ' (101) 

where R n is the approximate nuclear radius, 



Rr, 



0.9L4 1/3 + 0.3 



x 10~ 13 cm. (102) 



All of the dependence on the momentum transfer is then contained in the form factor. 

Determining o"o is then a matter of specifying the nature of the WIMP and its interactions. 
Effective models in which a scalar WIMP couples to a fermion through scalar or vector interactions 



(Equations 63 and 64) and have the right relic abundance of dark matter are already strongly 
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constrained by direct detection experiments (Beltran et al. 2009 Yu et al. 2012). For Dirac and 
Majorana fermions that have scalar interactions with quarks (Equation 67), the cross section is 



(o"o)5 
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(103) 



where C = 1(4) for a Dirac (Majorana) particle. Here Z is the number of protons, and A is the 
mass number of the nucleus. This model gives SI interactions. The factors Gg^ p ^ translate the 
coupling of the quarks to those of the protons and neutrons, 
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(104) 



The various factors in Equation 



0.003, 



104 



0.036 ± 0.008, and f^ n) 



are = 0.020 ± 0.004, = 0.026 ± 0.005, /£J = 0.014 ± 



0.118 ± 0.062 ( Jungman et al. 



1996). 



For a Dirac fermion with vector interactions to quarks (Equation 69 ) , the cross section is 
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where 



Gv, P 
Gv,n 



2Gv, u + Gy,d 
Gy,u + 2Gv,d 



(106) 
(107) 



This model also gives SI interactions. Note that the vector current vanishes for Majorana fermions. 



For a fermion with axial vector interactions to quarks (Equation 70 ) , the cross section is 
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(108) 



where C = 4 (16) for Dirac (Majorana) fermions, J a is the nuclear spin, and Sp and Sj{ are 
the expec tatio n values for the spins of the proton and neutron, respectively. The couplings in 
are given by G A ,( p , n ) = T,n=- ll .A.s G A,qA q p ' n \ and A p u 

\P An _ aP An _ 
^d' ^d — ^u, ^ s ~ 



Equation 
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-0.085, A! 7 



0.842. 



interaction has a similar form to Equation 108 with C = 16 



-0.427, A? 



Ag (Jungman et al. 1996). For a Dirac fermion, the tensor 



From Equations 103, 105, and 105, and assuming universal coupling to protons and neutrons, 
the cross section per nucleon is 



<?N = — A 0-0, (109) 

where m r> jv is the reduced mass between the WIMP and the nucleus and m r ^ is the reduced masses 
between the WIMP and the nucleon. 

The interactions that lead to Equations 103 105, and |108| are typically the most important 
for direct detection. The remaining interaction terms are suppressed by at least one power of the 



momentum or the velocity (e.g. Kurylov and Kamionkowski ([2004]) and Fan et al. (2010)). When 
the interaction is not suppressed at high momentum transfer, the SI interactions are enhanced 
relative to the SD interactions by the factors in the brackets proportional to A 2 . 

On the most general grounds, a Dirac fermion that couples to the Z° through vector interactions 



(Equation 105) has a significant SI cross section that was large enough to be ruled out by the 
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initial direct detection experiments nearly twenty years ago. This bound of course can be evaded 
for Majorana fermions since they do not have vector currents. It is straightforward to construct a 
Majorana fermion that has axial vector interactions and couples to the Z° and a scalar coupling 
to the Higgs; from the cross sections above these lead to SD and SI coupling respectively. For the 
later couplings, the scale for the scattering cross section is approximately a few times 10 -44 cm 2 ; 
this scale is now being probed by direct detection experiments sensitive to SI couplings (Section [9]). 
It is informative to compare to the predictions of different theoretical models. For example, 



in the scalar singlet model (McDonald, 1994 Burgess et al. 2001), the scalar WIMP couples to 



the nucleons through a coupling to the Higgs. In this model the cross section for elastic scattering 
is a function of the coupling of the Higgs to the nucleus and the coupling of the dark matter to 
the Higgs. In supersymmetric models, there are two primary interactions that survive in the low 
energy limit. The dominant contribution to SI interactions between neutralinos and quarks arises 
from Higgs and squark exchanges. Neutralino SD interactions on the other hand have contributions 



from both Z and squark exchanges (Jungman et al. 1996). 



8.3. Dependence on local dark matter density 

The dependence of the WIMP rate on the local dark matter density is straightforward- it is 
simply a normalizing factor that scales Equation 97 Recall from Section [3] that there has been 
a variety of recent determinations of the local dark matter density. Although it is typical for 
experiments to quote results assuming a canonical local dark matter density of 0.3 GeV cm -3 , in 
Section [3] it was seen that the local dark matter density may be up to three times larger than the 
canonical value of 0.3 GeV cm' 3 , though these results still systematically depend on the inputs of 
the analysis and the specific stellar population that is utilized (Garbari et al. 2011, 2012). Note 



that even these estimates are still below the local baryonic material by up to a factor of several, and 
they still carry both significant systematic and statistical uncertainties. Various other analyses that 
add in constraints from the total Galactic potential find that this uncertainty is reduced, and the 



mean central value for the dark matter density is slightly larger (Catena and Ullio, 2010 McMillan 



2011). 



The local kinematic measurements described above are primarily sensitive to the "smooth" 
distribution of dark matter in the Solar neighborhood. However, as discussed in Section [5] from 
theoretical predictions of dark matter halo formation in cold dark matter cosmology, the distribution 
of dark matter in the Galactic halo is not smooth, so in principle it is possible that the Sun resides 
in either a significant local over or under density of dark matter, which may affect the implied 



constraint on the WIMP elastic scattering cross section. Numerical simulations (Vogelsberger et al 



2009, Kuhlen et al. 2010), as well as analytic models (Kamionkowski and Koushiappas 2008), have 



begun to address this issue, finding that the probability for the Sun to reside in a significant over 
density or under density is small, of order 10 _4 %. Further, there are predictions of a dark matter 
disk in the Galaxy, that may have its origin in the accretion of a massive satellite galaxy that 



was dragged into the disk by dynamical friction (Read et al. 2008). However, analysis of stellar 



kinematics that extend out beyond a few kpc place strong limits on a dark matter disk component 



in the Galaxy (Bidin et al. 2010). So in sum, our modern theoretical expectation is that the event 



rates in WIMP direct detection experiments will be dominated by the smooth component, i.e. the 
component that is not locked up into subhalos or streams. 

8.4- Dependence on the velocity distribution 

The dependence of the WIMP event rate on the velocity distribution is more phenomenologically- 
interesting than the dependence on the local dark matter density. However, because it cannot be 



94 



directly inferred from stellar and gas kinematics like the local density is inferred, the velocity dis- 
tribution is essentially impossible to deduce from astronomical observations. Indeed, probably the 
only method to measure the local velocity distribution is through a direct detection experiment 
itself. This means that our understanding of the impact of the velocity distribution on the WIMP 
event rate relies entirely on theoretical modeling. 

As was shown in Section [5j numerical simulations do provide a theoretical starting point for 
understanding the dark matter velocity distribution. From a theoretical perspective, the simulations 
provide a Monte Carlo realization of the velocity distribution in a given volume. However, the 
resolution of the distribution in velocity space is affected by both the finite number of particles and 
the finite size of the volume within which the velocity distribution is determined. To understand 
the effect of the finite number of particles, define N' as the number of particles within the volume 
of consideration, and a as the one-dimensional characteristic width of the distribution of velocities 
in the halo. The characteristic separation in velocity space then scales as Av ~ a / '(N') 1 ^ 3 . From 
the definition of the velocity distribution function, the effect of the finite volume size scales as 
Av ~ a(N') 1 / 3 . So in order to accurately measure the velocity distribution from simulations, a 
compromise must be made. In position space, the highest resolution simulations of Milky Way mass 



2009 



Kuhlen et al, 2010), 



halos have up to 10 5 particles in a 2 3 kpc 3 box ( Vogelsberger et al. 
which implies a mean inter-particle separation of ~ 50 parsecs. For comparison, a terrestrial dark 
matter detector covers a distance of approximately 220 km/s xl yr ~ 2 x 10~ 4 parsecs per year. 
So clearly we are a long way from resolving the distribution function on the scales that are relevant 
for direct dark matter detectors, with even the highest resolution cosmological simulations. 

From Equation |97[ the WIMP event rate depends on the integral of the velocity distribution as 
g(vmin) = C, esc \f(v, ve)/v] d?v. This scaling can be understood by noting that the scattering cross 
section depends on the velocity as da/dE r ~ l/v 2 , and that the WIMP event rate is proportional 
to the mean WIMP speed. For the general case of anisotropic distributions, the function g(v m i n ) 
can be expressed as 



g{Vmin) 



dv'v' l dtt'f(v',n'). 



;no) 



For a given model of the velocity distribution, the uncertainty in the WIMP event rate due to the 
unknown velocity distribution depends on this function. 

Because direct detect experiments at present still operate in the "discovery" phase, and because 
the lack of resolution in simulations outlined above, in order to appropriately compare results 
experiments utilize a standard model for the velocity distribution. The velocity distribution function 
that the experiments use to interpret their reuslts is the standard halo model (SHM), which is an 
isotropic maxwellian distribution with a cut-off at the Galactic escape velocity, 

(v + v E ) 2 ' 
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where the normalization is 
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In the limit that v esc — > oo, the SHM reduces to the maxwell-boltzmann distribution in Equation [52j 
From the form of Equations 111 and 112[ there is a link between vq and the velocity dispersion a 
associated to any of the components of the velocity. This link can be make more clear by examining 
the circular speed vo as derived from the spherical and isotropic jeans equation, 
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Figure 19 Integral in Equation 110, normalized to the case with v m i n = 0. The left panel is for 
a germanium target, and the right panel is for a xenon target. The solid red curve assumes the 
standard halo model, and the dashed blue curve assumes an isotropic NFW model with a velocity 
distribution computed from Eddington's formula. For the NFW profile, the scale radius is 20 kpc. 
For the inner curves the WIMP mass is 10 GeV, and for the outer curves the WIMP mass is 50 
GeV. 

where the circular velocity is v c = \J GM/r, vq = \^2a, and p is the density. If the term in the 
brackets is Equation 1 13| is unity, the one-dimensional velocity dispersion that appears in the SHM 



velocity distribution is related to the cirular velocity as v 2 = 2a 2 . 

Although the SHM is a practical and versatile theoretical model for the dark matter velocity 
distribution, as was shown in Section[5]the velocity distribution from numerical simulations deviates 
from the SHM, and may not ultimately be an adequate description of the velocity distribution in our 
Galactic halo. Recall from Section [5] that simulations of dark matter halos (without the inclusion 
of baryonic physics) are converging on a form for the velocity distribution over a wide range of halo 
masses. In particular, the velocity distribution from simulations has a very high velocity tail that 
is steeper than a corresponding maxwell-boltzman distribution, a peak that is "flattened" relative 
to the more peaked maxwell-boltzman distribution, and a larger fraction of particles at velocities 



less than the peak. This behavior is quantified by Equation 53 



The difference between different velocity distributions, and their impact on direct detection 



rates, is quantified by Equation 110 (Fox et al. 2011). Figure 19 shows this function for two 



different models, for a WIMP mass of 10 and 100 GeV for a germanium and xenon target. The 
first model is a standard halo model with vq = 220 km/s and v esc = 553 km/s, and the second 
is an isotropic NFW model with a velocity distribution computed from Eddington's formula. For 
the NFW profile, the scale density is 8 x 10 6 M Q kpc -3 and a scale radius of 20 kpc. Converted 
to nuclear recoil energy, the most significant difference between the two models is at high nuclear 
recoil energy. 

8. 5. Determining WIMP properties 

Having developed an understanding of how the WIMP event rate depends on the scattering 
cross section, the local dark matter distribution, and the velocity distribution, it is now appro- 
priate to address the important question of what WIMP properties and Galactic halo properties 
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can ultimately be extracted from direct detection experiments. Because at this stage we only have 
upper limits on the dark matter cross section (as will be seen from the discussion in Section [9]), at 
this point in time this question will be addressed, most optimistically, by next generation detec- 
tors. However, there is some sense of urgency to improve the theoretical projections, because this 



information will both help to set the program for future direct dark matter searches (Section 10) 
and help to clarify the need for complementarity between direct searches, indirect searches, and 
colliders. The three main signals in direct detection experiments, specifically the mean event rate, 
the annual modulation signal, and the direction of a nuclear recoil, all will play an important and 
complementary role in determining WIMP properties. 

To gain some intuition on how well the WIMP mass can be determined, it will be useful to start 
by appealing to some simple scalings. Consider the spin-independent zero momentum transfer cross 
section. In this case, for no form factor suppression the event rate scales simply as the inverse of 
the WIMP mass, so lighter WIMP have a larger flux than heavy WIMPs. However, the lighter the 
WIMP mass, the larger the value of v m i n because the reduced WIMP-nucleus mass factor. This 
latter point implies that the range of velocities that is integrated over in Equation [97] is smaller 
for lighter WIMP masses, and that the flux is more sensitive to the shape of the high velocity tail 
of the WIMP distribution. It is easy to see in this case why the event rate is a sensitive function 
of the WIMP mass. For increasing WIMP masses, the suppression due to the nuclear form factor 
is more significant, so that the cross section is reduced from the zero momentum transfer value. 
In this case, for larger WIMP masses, (> 100 GeV), the event rate is less sensitive to the WIMP 
mass. Understanding more precisely how the above competing factors play off of one another will 
help better understand how well a WIMP of a given mass may be constrained from an observed 
sample of events. 

Initial detailed studies that focused on reconstructing the WIMP mass and cross section from 
direct detection data were based on fixed models of the WIMP velocity distribution and local dark 



matter distribution (Green, 2008). These results showed that for an exposure of 3 x 10 5 kg-day 
(~ 0.8 ton-yr), a 50 GeV mass WIMP with a spin-independent cross section of 10 -44 cm 2 could be 
determined with an accuracy of ~ 10%. The simple scaling arguments above are thus confirmed, 
in that for WIMPs > 100 GeV, the constraints weaken so that only a lower bound can be placed 
on the WIMP mass. 

In order to understand the origin of these results and to make more rigorous predictions, it will 
be useful to develop a simple model to determine the constraints on the WIMP mass and cross 
section as a function of the most relevant experimental parameters. Begin first by considering the 
SHM velocity distribution. If future experiments have a large enough sample of events over their 
range of recoil energies, it will be possible to bin the events to obtain spectral information. From 
this type of data set, the probability to obtain a given number of counts in an energy bin can be 
approximated as a product of poisson distributions, 

N N b 
b 7Vfe - 

where JVj is the number of observed counts in each bin, and Nr is the number of expected counts 
in the bin. The number of counts in each bin are determined from the theoretical model with 
parameters represented by 0. 

When binning the experimental data and comparing to a theoretical model, it is important to 
account for the finite energy resolution of the detector. This requires smearing the true event rate 
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spectrum according to 
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where dR/dE r is the true event rate spectrum as determined from Equation 97 and e is the energy- 
dependent energy resolution. The total number of events within a recoil energy range (E\,E2) is 
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where e e t / is the effective exposure as a function of energy. 

For the assumptions above from the poisson likelihood, we can define the elements of the inverse 
covariance matrix as 

p _ N exp dRi dR t 
1=1 



R t dO a de b ' 



;ii7) 



where N exp is the exposure of the experiment, i.e. the mass of the detector times the exposure 
time. From Equation 
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the one-sigma uncertainty on parameter a is F~ , evaluated at the 
fiducial values for the parameters. Here R t is the predicted rate in the i th energy bin. Though for 
the analysis in this section backgrounds are not included in the analysis, these can be added as 



necessary (see Section 10). 



Figure 20 shows the projected constraints on the WIMP mass and cross section, assuming the 
SHM velocity distribution with fixed parameters V c = 220 km/s and p@ = 0.3 GeV cm -3 . The 
results in Figure [20| provide a simple way of assessing how well, in an ideal case, the WIMP mass and 
cross section can be determined. Because the SHM is probably only an approximation to the WIMP 
velocity distribution, and because a range of uncertain Galactic model parameters contribute to 
the determination of the WIMP event rate, when a detection is confirmed a more rigorous model 
for the Galactic dark matter distribution will be required. Using a discrete binned model for the 



WIMP velocity distribution, Peter (2011) and Kavanagh and Green (2012) still show that there 



is a significant degeneracy in the reconstruction of the WIMP mass and cross section. Combining 
data from different nuclear targets does improve the constraining power on these quantities. 



Strigari and Trotta (2009) introduced a Bayesian method to determine the WIMP mass and 



cross section using the general parameterization for the Milky Way dark matter density profile in 



Equation 31 Combining the uncertainties in the five parameters of the profile in Equation 31 with 
uncertainties in the local dark matter density, the circular velocity, the dark matter halo mass, 
as well as the properties of the Milky Way disk, for a one ton-year exposure with liquid xenon, 



Strigari and Trotta (2009) find that the mass of a fiducial 50 GeV WIMP can be determined to a 



precision of less than approximately 25%. For WIMPs greater than 100 GeV, as above the mean 
event rate distribution will only be able to provide a lower bound on the WIMP mass. 



The Strigari and Trotta (2009) methodology developed a more robust connection between Galac- 
tic model parameters and the signal in direct dark matter detectors. This methodology can be 
improved upon by establishing a more theoretically-motivated connection between the dark matter 
density profile and the VDF. For isotropic models this can be readily accomplished through the 



use of Eddington's formula in Equation 13 In order to implement this method, the density profile 



parameters in Equation 31 are constrained from the astrophysical data, and the constraints on 
these parameters are then propagated through to determine the isotropic velocity distribution, and 
from this the WIMP event rate distribution. This method can be extended to a few anisotropic 
models for which it is possible to obtain an analytic VDF from a density profile. 
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Figure 20 Projected constraint on the WIMP mass and spin-independent cross section assuming a 
1 ton-yr exposure. The left panel shows the projections for a germanium target, and the right panel 
shows the projections for a xenon target. From left to right in each panel, the contours assume 
fiducial WIMP masses of 50, 100, and 200 GeV. Both the local density of WIMPs and the velocity 
distribution are fixed. 



Using a general parameterization for the isotropic WIMP VDF motivated by double power 
law density profile models, Pato et al. (2011) first examined the complementarity between ton- 
scale xenon, germanium, and argon targets. For assumed ton scale and greater exposures, the 
constraining power of argon, germanium, and xenon experiments were found to be able to reduce 
the uncertainty on the WIMP mass by a factor of two relative to utilizing just germanium and 
argon alone. Also, Pato et al. (2011) show that future experiments with different nuclear targets 
will be able to "self-calibrate," so that uncertainties on astrophysical parameters will be able to 
largely normalize out. 

It is also interesting to considering what can be learned about particle physics, beyond sim- 
ply the mass and the zero-momentum transfer cross section. With these motivations a useful 
general theoretical parameterization of the WIMP-nucleus cross section are effective field theory 
models (Fan et al. 2010; Fitzpatrick et al. 2012). Phrasing the cross section in these terms clearly 
shows how different components of the cross section depend on the nuclear recoil energy. 

The aforementioned methods outlined for determining WIMP properties using the mean event 
rate spectrum are complemented by using information on the annual modulation. The phase of 
the annual modulation signal is a function of recoil energy, such that for typical WIMP masses for 
lower recoil energy, there is a complete phase reversal of the signal. The recoil energy at which 
the phase shifts depends on the WIMP mass. This fact has been utilized to determine how well 



detectors are able to constrain the WIMP mass from the annual modulation signal (Lewis and 



Freese, 2004). More detailed determination of the potential to extract the WIMP mass from the 



annual modulation signal using the methodology of Strigari and Trotta (2009) and Pato et al. 



(2011) represents the next step forward in these studies. Finally, using distribution function based 
models and future directional dark matter detectors, Alves et al. (2012) discuss the potential to 
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reconstruct the WIMP distribution. These authors find that probably a sample of thousands of 
events is required in order to detect the presence of anisotropies in the dark matter halo. 
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9. Direct Detection: Experimental Results 



The first direct dark matter detection experiments that reported significant scientific results 
utilized germanium crystals, and were largely dedicated to other scientific purposes such as double- 
beta decay. However at the time when Dirac neutrinos were strong candidates for particle dark 
matter, there was some sentiment that the detectors would be able to detect these particles in 
abundance. Null results were obtained, though in the process these experiments did rule out 



the first theoretically well- motivated dark matter candidate (Ahlen et al., 1987, Caldwell et al. 



1988). As the technology has developed since these initial results, the field has been characterized 



by a continual stream of impressive results, and along with them a good deal of controversy. 
The first modern limits on spin-independent WIMP-nucleon cross section were reported by the 
Cryogenic Dark Matter Search (CDMS), which at the time operated at a shallow underground site 
at Stanford University. CDMS was able to rule out a cross section of a few times 10 -42 cm 2 as 100 
GeV ( Abusaidi et al. , 2000[ ). The improvement upon the initial germanium results was in large part 
due to improved background discrimination methods, in particular discrimination between electron 
and nuclear recoils. These limits were strong enough to essentially rule out an annular modulation 
signal that the DAMA experiment, which utilized sodium iodine crystals, claimed was due to dark 
matter (Bernabei et al. 2000). 

The past decade has been an even more extraordinary time for direct dark matter detection 
experiments. In addition to the development of larger and more sensitive germanium detectors, a 
variety of different technologies have been developed. Each of these technologies relies on different 
methods to eliminate backgrounds to extract the signal of a WIMP scattering off of a nucleus. 

This section reviews the results from modern direct dark matter detection experiments, focusing 
on those experiments that are now operating and have reported science results. The presentation of 
the results is broken up according to whether they were obtained from spin-independent searches 
or spin-dependent searches. Following this main line of discussion will be a short discussion of 
modern directional dark matter detection limits. For each of these three categories, the discussion 
is organized according to which experiment has the best limit at the time of the writing of this 
article. For much larger incarnations of these experiments, or for those that use technology that are 
currently not competitive with modern experiments, discussion is reserved until Section 10 After 
reviewing the experimental results, a connection is made to the both the theoretical models that 
were discussed in Section [8] and the results from indirect searches discussed in Section 



9.1. Spin-independent limits 

The organization of the spin-independent limits done according to which of the experiments 
have the best sensitivity. Note that some of the experiments discussed here also have reported 
competitive spin-dependent limits; where this is the case discussion of these results is reserved until 
the following subsection on spin dependent limits. 

• XENON100- The XENON100 experiment has been in operation in the Gran Sasso under- 
ground laboratory since 2009. It is the predecessor of the XENON10 experiment which operated 
in Gran Sasso from 2005-2007. It is a dual phase noble liquid detector that detects both primary 
scintillation light and secondary light proportional to the amount of ionization from particles that 
scatter in the detector volume. A recoiling particle induces ionization and excitation of Xe atoms. 
For a given amount of energy deposited into the detector, recoiling electrons produce a much larger 
amount of ionization than recoiling nuclei. In XENON100, an electric field is employed in the 
detector that drifts the ionized electrons into a region where the noble is in a gas phase, and there 
the electrons produce a signal proportional to the amount of ionization (this is the so-called "SV 
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Figure 21 Limits on the spin- independent WIMP-nucleon cross section as a function of mass from 
the experiments indicated. Plot produced by DM tools (http://dmtools.brown.edu:8080/). 



signal). The primary scintillation signal (the so-called "Si" signal) results from the production 
of photons from the de-excitation of Xe atoms that have not been absorbed within the detector. 
Electron and nuclear recoils are then discriminated based upon the ratio of S2/S1. 

The low energy threshold of approximately 8 keV for XENON100 is determined by photoelectron 
yield. The precise value of the WIMP mass that this corresponds to depends on the local WIMP 
velocity distribution (Sections [5] and [8]). Because of the large atomic mass of Xe, XENON100 reaches 
its best sensitivity at a lower WIMP mass than detectors that rely on smaller atomic masses. The 



energy resolution of XENON100 is < 1 keV at 1 keV and ~ 3 keV at 30 keV (Aprile et al., 2012). 



In July 2012, the XENON100 collaboration announced their results for 13 months of data 
between 2011-2012, with a total exposure of approximately 7665 kg-days. In the nuclear recoil 
energy range of 6.6 — 43.3 keV there were two events that were observed, at recoil energies of 
7.1 and 7.8 keV, which was consistent with the background expectation over this energy interval. 
From the lack of detected WIMP events over this energy interval, the deduced upper limit on 
the SI WIMP-nucleon cross section as a function of mass is given in Figure |2~Tj The XENON100 
limits are the strongest limits on the SI cross section for masses greater than approximately 7 GeV. 
XENON100 reaches its maximal sensitivity at approximately 40 GeV, at which the limit on the 
cross section is 10~ 45 cm 2 . 



• CDMS-II- The Cryogenic Dark Matter Search (CDMS-II) operated in the Soudan Mine from 
2001-2008. CDMS-II uses germanium and silicon detectors that reject electron recoils using a 
combination of phonon light and ionization. In a cryogenic detector, the struck ion dissipates 
energy via scatterings with electrons and ions on the lattice. The dissipated energy can be detected 
by either phonon production, ionization, or scintillation. Phonons are produced as the kinetic 
energy of the recoil is converted into a collective excitation of the crystal. The shape and timing of 
the phonon pulses from ejection of electron events near the surface and the ionization yield allows 
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for discrimination of bulk electron events from nuclear events within the detector. 

In December 2009, the CDMS-II collaboration announced the results of their data run from 
2007-2008 with an exposure of 612 kg-days. There were two candidate events detected in the nuclear 
recoil energy window of 10 — 100 keV, at energies of 12.3 and 15.5 keV, which was statistically- 
consistent with the expected background over this energy interval. At this time of announcement 
this was the best limit for WIMP masses greater than approximately 40 GeV, though now these 
results have been superseded by the aforementioned XENON100 results at all WIMP masses. 

In order to improve the sensitivity to WIMPs in the lowest mass regime that is probed by 
CDMS-II, the collaboration has presented an analysis with its threshold energy lowered to 2 keV. 
Since the discrimination between nuclear and electron recoils significantly degrades at these lower 
energies, the limits on WIMP masses less than approximately 10 GeV is degraded relative to 
WIMPs of larger mass. Even with the larger background acceptance, with an exposure of 241 
kg-days CDMS-II was able to improve the limits on the SI WIMP-nucleon cross section below a 
mass of approximatley 9 GeV (Ahmed et al. 2011b), and exclude cross section > 3 x 10 -40 cm 2 
for WIMP masses greater than 5 GeV. 



• EDELWEISS-II- EDELWEISS-II is a cryogenic germanium detector located in Laboratoire 
Souterrain de Modane. EDELWEISS-II has a different strategy for rejecting surface backrgound 
events relative to the method employed by CDMS-II. From an exposure of 384 kg-day in 2009-2010, 
EDELWEISS-II detected five events above a recoil energy of 20 keV, four of which have energies 
less than 25 GeV, and one has an energy of 172 GeV. This is also consistent with the expected 
number of background events (Armengaud et al. 2011). For WIMP masses larger than about 85 



GeV, EDELWEISS-II achieves a sensitivity that is similar to that of CDMS-II. EDELWEISS-II 
and CDMS-II have recently combined their data sets, achieving a total of 614 kg-days of effective 



exposure, improving the respective individual limits by a factor of approximately 1.6 (Ahmed et al. 



2011a). 



EDELWEISS-II has also presented a low threshold analysis with a goal of improving the sen- 
sitivity to WIMPs with mass below 10 GeV (Armengaud et al. 2012). From an exposure of 113 
kg-days and the observation of one event in the signal region, at 5 GeV the limit on the cross 
section is approximately three times weaker than the corresponding limits from CDMS-II. 



• ZEPLIN-III- ZEPLIN-III is a liquid xenon experiment that operates in the Boulby Under- 
ground Laboratory. Similar to XENON100, it is a dual phase detector that relies on the detection 
of both ionization and scintillation light. From an exposure of 1344 kg-days between 2010-2011, 
ZEPLIN-III observed a total of eight events in the nuclear recoil energy range of 7 — 29 keV, again 
consistent with its expected backgrounds. Below approximately 70 GeV, ZEPLIN-III achieves 
better sensitivity to the Si-cross section than CDMS-II or EDELWEISS-II do individually. 



• CRESST-II- The Cryogenic Rare Event Search with Superconducting Thermometers (CRESST) 
is located in Gran Sasso National Laboratory, using a scintillating material of Calcium Tungstate 
(CaW04). CRESST-II uses a combination of phonon and scintillation light to reduce backgrounds, 



and now has an exposure of 730 kg-days (Angloher et al. 2012). It has observed 67 events in its 



acceptance nuclear recoil energy range of 10 — 40 keV. The indication from the analysis is that 
the known backgrounds are unable to account for these events observed in the signal region; from 
a maximum likelihood analysis, the CREST-II results can be made consistent with WIMPs with 
mass less than 25 GeV. However, the cross sections which are able to explain the CRESST-II data 
are inconsistent with the CDMS-II and XENON100 limits described above. 
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• DAMA- The DArk MAtter (DAMA) experiment is located in Gran Sasso and uses radiopure 
Nal (TI) scintillators as target detectors. DAMA does not distinguish between a WIMP signal and 
backgrounds on an individual basis, rather it infers the existence of WIMPs from single scatter 
interactions at low energy. 

The DAMA/Nal experiment has collected seven years of data between 1995-2002 with a mass of 
100 kg and a total exposure of 0.29 ton-yr. Including the second phase of DAMA/LIBRA, the total 
mass is 250 kg and the total exposure is 0.53 ton-yr. In the 2 — 6 keV energy interval, DAMA/Nal 
and DAMA/LiBRA report a 8.2a detection of an annual modulation signal with both a phase and 
amplitude that is consistent with what is expected from a dark matter signal. No modulation is 
seen in the multiple scatter events, which can be interpreted as evidence that the signal is not from 
a modulating background. 

There have been numerous efforts to reconcile the observed DAMA signal with other null results, 



and windows have remained open for low- mass spin- independent interpretations (Gondolo and 



Gelmini, 2005) until the recent announcement of the recent XENON100 results. There is manifest 



tension between these results, and it is not yet clear how to explain the DAMA results in light of 
all of the other dark matter limits. There is still the possibility that the reported results arise from 



an unknown experimental background; such possibilities are nicely laid out in Schnee (2011). 



• CoGeNT- CoGeNT (short for Coherent Germanium Neutrino Technology) is a germanium 
detector cooled to liquid nitrogen temperatures of 77 K, with sensitivities to sub-keV nuclear re- 



coils (Aalseth et al. , 2011b ) . These experiments only measure deposits of energy through ionization, 



and do not discriminate between electron and nuclear recoil backgrounds. 

In 2010, CoGeNT identified an excess of events at energies less than 3 keVee, where here the 
energy units stand for equivalent electron energy. At the time it was determined that, if interpreted 
as a WIMP-induced spectrum, this excess is compatible with the DAMA signal of a 7 GeV WIMP. 
As is the case of the DAMA signal, there is significant tension between this result and the recent 
results presented by the XENON100 collaboration. In 2011, CoGeNT presented 2.8a evidence for 



an annual modulation signal that is consistent with a WIMP of mass 7 GeV (Aalseth et al. , 2011a). 



CDMS-II searched for an annual modulation signal in its nuclear recoil energy range of 5 — 11.9 keV 
and reported null results that appear in conflict with the CoGeNT and DAMA annual modulation 



signals (Ahmed et al. 2012). 



TEXONO- TEXONO is a germanium detector similar to CoGeNT, which is located in the 



Kuo-Sheng Laboratory ( |Lin et al. 2009). From an exposure of 0.338 kg-day, TEXONO has reported 
limits on the spin-independent cross section down to 2 GeV, constraining the cross section to 



< 10 38 cm 2 at this mass. 



• XENON 10- XENON10 has published a limit on low mass WIMPs that does not include 
the primary scintillation "SI" signal (Angle et al. 2011). Without this signal, it is unable to 



discriminate nuclear from electron recoils, so the sensitivity is reduced, however it does allow for 
nuclear recoil energy to be measured by the electron signal. This provides sensitivity down to 
approximately 1 keV nuclear recoil energy. As shown in Figure 21, this limit is in significant 
tension with the results that are interpreted as a WIMP detection at a mass near 10 GeV. 



9.2. Spin- dependent limits 

This subsection now moves on to discuss spin-dependent limits. As above the results are or- 
ganized based upon which experiments have the best limits as of the writing of this article. The 



updated limits discussed below are shown in Figures 22 and 23 
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Figure 22 Limits on the spin-dependent WIMP-neutron cross section as a function of mass, 
from the XENON10 Collaboration ( |Angle et ~aT| |2008[ ). Plot produced by DM tools 
(http://dmtools.brown.edu:8080/). 



9.2.1. Direct detection 

Begin by considering the so-called threshold detectors (SIMPLE, PICASSO, and COUPP); 
these operate by tuning the temperature and pressure of a superheated liquid to a point at which 
they are insensitive to electron recoils. Though these detectors provide limited information on the 
energy of an event, they do hold the promise of being readily scalable to larger volumes. 



SIMPLE- The Superheated Instrument for Massive ParticLe Experiments (SIMPLE) (Fe- 



lizardo et al. 2010) operates in the Low Noise Underground Laboratory in southern France. It has 
sensitivity to spin-dependent WIMP-proton interactions through scattering on X 9F, 3 5C1, 3 7C1, and 
1 3C. Combining the exposures for Stage 1 (13.47 kg-days) and Stage 2 (6.71 kg-days) of the SIM- 
PLE experiments, it is able to exclude a WIMP-proton spin-dependent cross section > 5.7 x 10~ 39 
at 35 GeV. 



PICASSO- The Project in Canada to Search for Supersymmetric Objects (PICASSO) (Ar- 



chambault et al. 2012) is located at SNOLAB in Canada; it is sensitive to spin-dependent WIMP- 
proton interactions primarily through scattering on 19 F. In 2012 PICASSO reported results from a 
total exposure of 114 kg-days, and from these results constrains the cross section to < 3.2 x 10~ 38 

2 

cm . 



COUPP- The Chicagoland Observatory for Underground Particle Physics (COUPP) (Behnke 



et al. 2011) operates in the MINOS tunnel at the Fermi National Accelerator Laboratory. It is 
sensitive to spin-dependent WIMP-proton interactions primarily through scattering on flourine. 
From a 28.1 kg-day exposure in 2009, COUPP achieved a maximal sensitivity to the cross section 
at a WIMP mass of approximately 50 GeV. 
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Figure 23 Limits on the spin-dependent WIMP-proton cross section as a function of WIMP mass. 
Plot produced by DM tools (http://dmtools.brown.edu:8080/). 



• KIMS- The Korea Invisible Matter Search (KIMS) is located in the Yangyang Underground 
Laboratory. KIMS uses Csl detectors and a pulse shape discrimination method to statistically- 
discriminate between nuclear and electron recoil events. From an exposure of > 24524. kg-days 
between 2009-2010, KIMS has recently placed an upper limit on the spin-dependent WIMP-proton 



cross section that is more strict than any other experiment for WIMP masses > 20 GeV (Kim 



et al. 2012). 



XENON10- The most stringent limits on the WIMP-neutron spin-dependent cross section 

From an exposure of 136 kg-days, 



come from the XENON10 experiment (Angle et al. 



2008[ ) 

XENON10 achieves a maximal sensitivity of < 10~ 38 cm 2 at approximately 25 GeV (Figure 22). 



9.2.2. Neutrinos 

As discussed in Section [6j WIMPs get captured into the Sun and annihilate, and the only 
particles that can effectively escape to be detected are neutrinos. The neutrinos then interact in 
the Earth to create muon neutrinos which are then detectable via Cerenkov radiation. As a result 
neutrinos from the Sun are particularly effective at probing the spin-dependent WIMP-nucleon 
cross section. 

Both Super-Kamiokande (Tanaka et al. 2011[ ) and Ice-Cube ( Abbasi et al. , 2012a) have searched 
for neutrinos that come from dark matter capture and annihilation in the Sun. As in the case of 
diffuse neutrino searches discussed in Section [7J neutrino experiments quote their limits for hard 
(W + VF~) and soft (bb) spectra for the sources of neutrinos. For Super-Kamiokande, the total 
exposure time for the search was approximately 3110 days, while for Ice Cube the total exposure 
was 1065 days. For both searches, no excess of events above backgrounds were reported, leading 
to flux upper limits of muon neutrinos from the Sun. For Super-Kamionkande, assuming a dark 
matter mass greater than 100 GeV, the flux upper limits are approximately 6 x 10~ 15 cm~ 2 s -1 
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Figure 24 Limits on the spin-dependent WIMP-proton cross section as a function of WIMP mass, 
including neutrino limits from Super-Kamiokande and Ice Cube. 

for the soft channel, and approximately 4 x 10~ 15 cm~ 2 s^ 1 for the hard channel. For the soft 
channel, the flux limits extend down to 20 GeV, and are a factor of two larger than the limits 
above 100 GeV. For Ice Cube, the reported limits cover a dark matter mass range of 50 to 50000 
GeV, with the corresponding muon neutrino flux upper limits of approximately 10 4 km -2 yr _1 
to 10 2 km -2 yr 1 over this range. The limits on the WIMP-proton spin-dependent cross section 



from Super-Kamiokande and Ice Cube are shown in Figure 24 in comparison to the most stringent 



direct detection limits from Figure 23 



9.3. Directional searches 

Having surveyed the status of spin-independent and spin-dependent searches, now briefly con- 
sider the limits from dark matter detectors with directional sensitivity. The best limit on dark mat- 
ter from this class of detector is from the Dark Matter Time Projection Chamber (DMTPC) ( |Ahlen 



et al. 2011). The DMTPC is a CF4 gaseous detector, with a sensitivity to spin-dependent interac- 



tion of WIMPs off of flourine. The first results from the DMTPC ( [Allien et al.[|201ip utilize a 35.7 
kg-day exposure, which is several orders of magnitude less than the exposure of the experiments 
reported in this section. They obtain an upper limit on the spin-dependent cross section of ~ 10~ 32 
cm~ 2 . The threshold energy of the DMTPC is much higher than the experiments described above, 
~ 80 keV, so that the best sensitivity is obtained for WIMP masses ~ 100 GeV. 

9-4- Summary of experimental results and interpretations 

Taking a step back now and surveying the results presented in this section, it is clear that 
the direct dark matter searches are moving at a rapid pace. DAMA, CoGeNT, and CRESST 
have shown that their results cannot be explained by identified backgrounds. While each of these 
individually may be reconciled with a WIMP mass in the ballpark of 10 GeV, there is no hint of a 
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signal at this mass scale from either XENON100 or CDMS-II. While there have been models that 



have been put forward that attempt to reconcile these results (Feng et al. 2011), it is not clear 



that these satisfy all of the most recent constraints, even taking to account uncertainties in the 
dark matter local density and velocity distribution (Fox et al. 2011; Frandsen et al. 2012b) 

The limits presented in this section have significantly improved constraints on dark matter 



with effective interactions (Beltran et al. 



, 2009 


Zheng et al. 2012 


Yu et al. 


2012 



Kumar et al. 



2012b). Scalar and fermionic WIMPs with mass < 1 TeV with scalar interactions and universal 



couplings with fermions are now strongly excluded by direct detection results. Note that fermionic 
WIMPs with scalar interactions have velocity-suppressed annihilation cross sections, so they are 
not detectable with Fermi or other indirect detection experiments, so direct detection is the only 
means to probe these models. Fermionic WIMPs with vector interactions (such as the heavy 
Dirac neutrino) are also strongly excluded by direct detection; these models do not have velocity- 
suppressed annihilation cross sections so they are also constrained by indirect detection results. 
Fermionic WIMPs with scalar interactions and universal couplings to fermions are now also strongly 
excluded; for mass-dependent couplings to fermions the limit on the WIMP mass is < 200 GeV. Note 
again that direct detection experiments are not sensitive to WIMPs with pseudo-scalar interactions. 

Direct detect searches, in particular the new XENON100 limits, are now reaching into the well- 
motivated focus point regime of the supersymmetric model parameter space. Though searching 
the entire supersymmetric parameter space for models that match relic abundance constraints and 
experimental constraints is challenging ( Akrami et al. , 2010 Buchmueller et al. , 2012 ), from general 
calculations of the neutralino-proton scattering cross section the spin-independent cross section is 



right where the XENON100 constraints are sensitive to for neutralino masses > 100 GeV (Feng and 



Sanford, 2011 ). This conclusion is further strengthened if the Higgs mass is ~ 126 GeV (Akula et al. 



2012), and tells us that in the coming years it is very plausible that direct detection experiments 



will see hints of neutralino dark matter from supersymmetry. 
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10. Future Directions and Complementarity with Other Searches 



This section reviews expected developments in dark matter searches over the next decade. 
Several issues are addressed, broken up into the following components: future developments in 
direct and indirect searches, complementarity with recent results from collider and cosmological 
searches, and finally forthcoming astronomical observations 

10.1. Direct detection 

There are a variety of direct dark matter searches that will improve on the limits presented in 
Section [9] during the next decade. Several of these experiments are extensions of modern experi- 
ments, while others represent new technological directions. 

• Xenon- As discussed in Section [9j liquid xenon experiments are readily scalable, so that an 
order of magnitude increase in the fiducial volume is expected within the next few years. The next 
phase of the XENON100 experiment is Xenon 1 ton, which is fully funded and will be located 
in a tunnel in Gran Sasso near the XENON100 experiment. As its name indicates, the fiducial 
active volume of this experiment will be 1 ton. In addition to the increased volume, additional 
improvements are expected to increase the experimental sensitivity. These improvements include 
a water shield to be used as a muon veto to reduced backgrounds. Further there is expected to 
be a reduction of Kr and Rn backgrounds via the development of removal towers. Xenon 1 ton is 
currently in the planning and constructing phase, with first data expected in 2015. If the projected 
sensitivity is achieved, it is expected to cover a large portion of the so-called focus point region of 



supersymmetric parameter space (Feng et al. 2000). 



In addition to Xenon 1 ton, the Large Underground Xenon (LUX) experiment, which is located 



in the Homestake mine in South Dakota, is expected come online soon (McKinsey et al. 2010). 
Because it is liquid xenon, the principles are similar to those of the XENON100 experiment. LUX 
will have a 100 kg fiducial volume, but a background that is suppressed relative to XENON100 
because of 200 kg of LXe that acts as self-shielding. Additionally, water shielding around the 
detectors act as a veto against muons. Data from its first science runs is expected in 2013. 

• SuperCDMS- SuperCDMS represents the next generation of the CDMS-II experiment. Su- 
perCDMS will also be located in the Soudan mine, with a goal of eventually moving to the deeper 
SNOLAB facility in Sudbury, Ontario. The individual detectors of SuperCDMS will be increased in 
mass by greater than about a factor of two and, and it will have improved background discrimina- 
tion relative to those of CDMS-II. The total target mass of the detectors is expected to be 100-200 

• Argon and Neon- Like liquid xenon detectors, liquid argon detectors are excellent scintillators 
and ionizers. Liquid argon detectors are easily scalable, though there is a significant contaminating 
radioactive background from 39 Ar. There are several collaborations pursing a ton-scale liquid argon 
detector including Argon Dark Matter (ArDM) |^] and the WIMP Argon Program (WArP) 

Like both xenon and argon, neon offers potential for good background rejection through pulse 
shape discrimination. Neon is being implemented as part of the DEAP/CLEAN program, with 
a ~ 360 kg prototype, mini-CLEAN, currently being built ( Hime| |2011 ). Liquid neon has the 
additional benefit of being scalable to approximately 100 ton. 



2 http:/ /cdms. berkeley.edu/ 
3 neutrino.ethz.ch/ArDM/ 
4 warp .lngs.infn.it/ 
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As a final note on direct dark matter detectors, recall that they were first conceived as Solar 



neutrino experiments (Drukier and Stodolsky 1984). This is because if a neutrino strikes a nucleus 



a coherent recoil happens if the three-momentum imparted to the nucleus, q, obeys the inequality 
qR n < 1, where R n is the radius of the nucleus. For typical nuclear radii, coherent scattering 
leads to nuclear recoils in the range of a few keV for incoming neutrinos with energies in the range 
~ 1 — 100 MeV. Though this neutrino-nucleus coherent scattering remains a fundamental prediction 
of the standard model, this process has yet to be detected. 

There are several sources of neutrinos that are expected to appear in next generation under- 



ground detectors that will serve as backgrounds to direct dark matter detection signals (Monroe 



and Fisher 2007 Strigari, 2009). These include Solar neutrinos, atmospheric neutrinos, and diffuse 
supernova neutrinos (Horiuchi et al. 2009; Beacom, 2010). The neutrinos produced from the 8 B 



decay process are the most-easily detected component of the Solar neutrino flux. In the near term, 
direct detection experiments with sensitivities to nuclear recoil energies of a few keV or less will be 
most sensitive to this signal. 



Figure 25 shows the expected neutrino fluxes in xenon and germanium detectors as a function 
of nuclear energy. Over the entire nuclear recoil energy range, atmospheric and diffuse supernova 
neutrinos rival the WIMP event rate for a WIMP of mass 100 GeV and cross section of < 10~ 48 cm 2 . 
Though the neutrino signal mimics the WIMP signal in many ways, some differences between the 
signals can be exploited to distinguish the neutrino signal from the WIMP signal. The most obvious 
method is to use directional dark matter searches, which would clearly indicate the presence of a 
WIMP signal due to the the motion of the Earth through the Galactic rest frame. By comparison 
both the diffuse supernova neutrino and atmospheric signal are completely isotropic. Even in 
the absence of a directional signal, it may be possible to perform a multi-component spectral fit 
to the neutrino spectra, since the Solar and atmospheric neutrino spectra are well understood 
from previous neutrino experiments. It is also interesting to consider the possibility for extracting 
new physics from the coherent neutrino signal, for example as a probe of non-standard neutrino 



interactions (Pospelov, 2011 Harnik et al. , 2012 ). Note that in their low threshold analysis, CDMS- 



II has reduced the energy threshold to ~ 2 keV, making them at this stage closest to measuring 
the coherent scattering from 8 B Solar neutrinos. If the exposure in Ahmed et al. (2011b) were 
extended down to 1 keV, the 8 B event rate in CDMS-II would be one per ten years. 

Next generation dark matter detectors will also be sensitive to scattering of Solar neutrinos with 



electrons. Though the rate for this interaction is lower than the nuclear recoil rates in Figure 25 



this scattering on electrons is dominated by the neutrinos created directly in the proton-proton 
chain deep inside the Solar interior. Though in terms of number of neutrinos over all energies, the 
proton-proton component of the Solar neutrino flux dominates, neutrinos of this energy are too low 
to be detected with modern detectors. The neutrinos from the proton-proton chain would make a 
flat spectra of electron recoils that extends up to electron energies ~ MeV (Mailing et al. 2011). 



10.2. Indirect detection 

Indirect dark matter detection experiments will also significantly improve their sensitivity to 
WIMPs over the course of the next several years. Focusing on gamma-ray searches, the progress 
can be characterized along two different fronts: more and improved data from the Fermi-LAT, and 
larger and more sensitive ground based ACTs. 

• Fermi-LAT— The Fermi-LAT will improve the results presented in section[7]for several reasons. 
For dSph searches along the lines of the current analysis, the limits will clearly improve linearly 
with exposure time in the high mass WIMP regime of a few hundred GeV. With a total of ten years 
of data, extremely sensitive limits or detections of dark matter will be obtained for the entire mass 
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Figure 25 Neutrino and WIMP coherent scattering rates in xenon (left) and germanium (right). 
The neutrino rates are determined as in Strigari ( 2009 ) . The two WIMP spectra assume a mass of 
100 GeV and cross section per nucleon as indicated. 



regime of 10 — 700 GeV, and there is even reason to believe that better sensitivity will be achieved 
than can be currently projected. For example, galaxy surveys that are now coming online, such 
as Pan-STARRS, the Dark Energy Survey (DES), and farther into the future the Large Synoptic 
Survey Telescopes (LSST), are certain to discover more faint satellite galaxies. As new objects 
are detected, kinematic analysis of their constituent stars will provide their dark matter masses, 
and adding these objects to the Fermi-LAT analysis will certainly improve on the current limits. 
In addition, there will be improvements in sensitivities on all sources, and concurrent theoretical 
advances that will improve the understanding of their dark matter and astrophysical gamma-ray 
emission. Even extending beyond Fermi-LAT, Gamma-400, a Russian space-based mission that is 



designed to detect gamma rays in the energy range of 100 MeV to 3 TeV (Galper et al. , 2012), will 
complement the above results. 

• ACTs- There is room for significant improvement in dark matter searches with ACTs. Modern 
experiments such as MAGIC, VERITAS, and HESS will continue to take data and improve the 
results presented in section [7j The HESS experiment is being expanded with the addition of a 
fifth 600 m 2 telescope (HESS-II), which will both expand the effective area and reduce the energy 
threshold. Extending farther into the future, the proposed Cherenkov Telescope Array (CTA) 
which is expected to begin operation near 2017, will extend WIMP dark matter limits to much 
higher masses, beyond the TeV scale, with better angular resolution than Fermi-LAT. CTA is 
expected to improve present limits from dSphs by at least an order of magnitude, and extend the 



Fermi-LAT studies of the Galactic center to a much higher energy regime (Doro et al. 2012) 



Combining the satellite searches with the results from the Galactic center, line searches, and 
the other sources presented here through multiple wavelengths, the analysis will only improve over 
the next several years. In addition, we of course cannot anticipate surprising discoveries that will 
be made with more Fermi-LAT data, as well as forthcoming data from larger ACTs and neutrino 



www.cta-observatory.org/ 
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Figure 26 CMS and ATLAS limits on the WIMP-nucleon spin-dependent (left) and spin- 
independent cross section (right) from monophoton searches. Data from Chatrchyan et al. (2012) 
andlAad et all fl2012b|). 



experiments. Finally with regards to indirect detection experiments, new anti-matter detectors, 
in particular AMS-2, will begin reporting results, and test the increase in the positron fraction 
observed by PAMELA. 

10.3. Complementarity with initial LHC searches 

The effective model of dark matter interactions discussed in sections [6] and [8] predicts interesting 
signals at colliders, which are complementary to indirect and direct searches. For couplings of dark 



matter to fermions, it is possible to translate the interactions in Equations 66 67 to annihilation of 
fermions into dark matter and constrain the values of the respective couplings. Events in which a 
pair of dark matter particles are produced can be detected in colliders if a photon or a jet emanates 
from an incoming fermion. 



Recently the CMS (Chatrchyan et al, 2012) and ATLAS (Aad et al., 2012b) detectors have 



reported results on searches for so-called monophoton events. These events are characterized by a 
photon and missing transverse momentum in the final state; the missing transverse momentum in 
this case is the probe of new physics such as the production of dark matter (without the accom- 
panying photon, the missing energy signal is unobservable at a hadron collider). CMS observed 73 
of these monophoton events, which was determined to be consistent with the background expected 
from standard model processes. ATLAS observed 116 events, again consistent with its background 
expectations. The consistency with the standard model background implies an upper bound on 
the spin-independent and spin-dependent scattering elastic scattering cross sections. The limit on 
the former cross section is determined from scalar and vector couplings, while the later is deduced 
from axial- vector and tensor interactions. The upper limits are shown in Figure 26, in comparison 
to the limits from direct detection discussed in section [HJ Over the entire WIMP mass range, the 
limits from ATLAS and CMS are now more stringent than the spin-dependent limits deduced from 
direct detection. Above about 1 TeV, the only more stringent constraint is the Ice Cube W + W~ 
limit. For WIMP masses around 1 GeV, the ATLAS and CMS limits are more strict than the 
spin-independent limits from direct detection. 
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Limits on monophoton searches are also complementary to the Fermi-LAT limits from dSphs 
(Figure 15). When translated into a limit on the annihilation cross section, the ATLAS and CMS 
limits on WIMPs with mass of 10 GeV with vector current interactions described by the operator 
in Equation 69 are just as stringent as the Fermi-LAT limits. Further, the implied ATLAS and 
CMS upper limits on 10 GeV mass WIMPs with axial vector interactions described by the operator 



in Equation 70 are an order of magnitude more stringent than the Fermi-LAT bounds (Aad et al. 



2012a). Above a WIMP mass of 100 GeV the Fermi-LAT bound is more stringent than the ATLAS 



and CMS bounds. 

The aforementioned limits were obtained with 5 fb _1 of data for CMS and 4.6 fb _1 of data for 
ATLAS. With the increase of the luminosity at the LHC, more stringent limits on monophoton 
events will be set, and also the search will continue for supersymmetric dark matter through the 
detection of colored particles, quark jets, hard leptons, and missing transverse energy that signifies 
the presence of physics beyond the standard model. 



10. 4- Cosmological searches 

In the early universe, WIMPs annihilate and inject energy into the photon-baryon plasma, 
ionize the gas and increase the number of free electrons. This process leaves an imprint on the 
CMB anisotropies, is most sensitive to WIMPs that do not annihilate largely into neutrinos, and is 
particularly sensitive to WIMPs that annihilate directly to electrons. Several authors have recently 



developed the theory for energy injection by WIMPs (Galli et al. 2009; Slatyer et al. 2009 Huetsi 



et al. 2009; Finkbeiner et al. 2012); these authors show how the constraints on the WIMP mass 



and cross section are particularly sensitive to the fraction of the energy that is absorbed by the 
plasma. Using WMAP7 data, for direct annihilation into electrons the limits on the annihilation 
cross section reach the thermal relic scale at a WIMP mass of approximately 10 GeV, and are a 
factor of a few weaker for direct annihilation into muons (Galli et al. , 2011 ). Data from the Planck 
satellite mission, which will be released in 2013, is expected to further improve these limits by 
about an order of magnitude. 



10.5. Astronomical observations 

If WIMPs are not detected with the next decade, it is likely that there will be a paradigm shift 
in how dark matter is viewed in astrophysics, cosmology, and particle physics. In this circumstance 
of events, it will becomes more important to understand what astronomical surveys can teach us 
about the nature of dark matter. As a general point, over the course of the next several years, the 
improvement in the astronomical observations will be important from two perspectives. On the 
one hand, we will have both better measurements of astrophysical parameters as input to direct 
and indirect detection experiments. In addition, observations that are sensitive to deviations from 
the standard CDM model will become more relevant. 

• Milky Way surveys- The GAIA satellite is expected to launch in 2013. GAIA is the next 
generation large scale astrometric survey of the Milky Way which will significantly improve the 
corresponding measurements from the Hipparcos satellite. For stars brighter at 20th magnitude 
GAIA will have an astrometric accuracy of approximately 1 mas and a photometric accuracy of 60 
mmag. From the GAIA data, it will be possible to re-analyze the local distribution of dark and 
luminous matter via the methods discussed in section |3] and to obtain better constraints on the 
local dark matter density. Further, GAIA will be able to probe the shape of the Milky Way dark 
matter halo from kinematic information on the population of stars in the outer halo, improving the 
measurements from current disk populations of stars discussed in section [3] It will also be able to 
identify streams and other features in the phase space of stars. 
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• Kinematics of Satellites- Larger samples of resolved stars in dSph satellite galaxies with 
kinematic information will continue to improve the measurements of their dark matter distributions. 
More line-of-sight velocities will be important to further test the initial results that have been 
obtained from multiple populations of stars in a small number of dSphs. Also, new ultra-faint 
satellites that will be discovered in future surveys will require kinematic follow up. It is very possible 
that a new survey will detect an object that is as massive and nearby as the SDSS discovery of 
Segue 1, so follow up on all satellites that are discovered will be important. Beyond measurements 
of line-of-sight velocities, in the future it will also be possible to obtain transverse velocities of 
bright stars in several dSphs. This will go a long way towards breaking the degeneracy between the 
velocity anisotropy and the dark matter mass profiles, and distinguishing between cold and warm 



dark matter models (Strigari et al. 2007a) 



• Surveys for Satellite Galaxies- A variety of large-scale galaxy surveys are beginning to take 
data, or will take data, in the very near future: these include Pan-STARRS, Dark Energy Survey 
(DES), and even farther into the future the Large Synoptic Survey Telescope (LSST). Pan-STARRS 
will have a survey region that largely overlaps with the SDSS sky region, however, DES will cover 
about one-half of the Southern Sky that was not surveyed by SDSS, and LSST is scheduled to 
survey the entire southern sky down to a magnitude limit of g = 24. These surveys are each 
expected to make significant contributions to Galactic astrophysics and the studies of the dark 
matter distribution in the Milky Way. Though it is impossible at this stage to predict how many 
ultra-faint satellites will be detected in these surveys, only a very small fraction of the total Milky 
Way volume has been surveyed to the SDSS depth, so there is certainly a vast amount of discovery 



space available (Koposov et al. 2008 Tollerud et al., 2008 Willman, 2010). 



As discussed in section [5j surveys for faint satellites around external galaxies are now reaching a 
significant level of maturity. The results obtained to date will be improved by a combination of the 
aforementioned surveys. For example, it is likely that these surveys will identify more than 4000 
host galaxies with satellites to Fornax magnitude differences, approximately 1600 galaxies with 
satellites as faint as Leo I, and approximately 100 host galaxies with satellites as dim relative to 
their primary as Sculptor is to the Milky Way. A large statistical sample along these lines will also 
go a long way towards probing the subhalo population around other galaxies, and distinguishing 



between cold and warm dark matter models (Strigari and Wechsler 2012) 



• Gravitational Lensing- As discussed in section [5j strong gravitational lensing is a powerful 
probe of the mass distribution of galaxies. Strong lensing is particularly interesting as a probe 
of dark matter substructure, providing the only method for detecting substructure outside of the 
Local Group. When a background quasar is lensed by a foreground galaxy, at optical and X- 
ray wavelengths the observed flux ratios are affected not only by dark matter substructure, but 
also by stars in the lensing galaxy. Because the contribution of the stars to the lensing signal 
varies, the contribution to the flux measurement from the stars is typically subtracted out by long 
term monitoring of the lightcurves. At radio wavelengths however, microlensing by stars is not 
important because the emission comes from a larger region within the source. Subtracting out this 
stellar emission through multi- wavelength studies of strong lensing systems is a promising technique 



for measuring the properties of dark matter substructure in these systems (Moustakas et al., 2009). 



Time delay measurements of strong lensing systems may probe the substructure mass function 
down to mass scales < 10 9 M (iKeeton and Moustakasl |2009|) . 



One can additionally ask whether dark subhalos in our own Milky Way Galaxy may cause a de- 
tectable gravitational lensing signal. The primary problem with lensing signatures in the Milky Way 
is that subhalos are predicted to be too diffuse- only if subhalos happen to have inner slopes that 
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are near that of a singular isothermal sphere will they generate a photometric microlensing signal. 
The situation is slightly better when considering an astrometric signal, however this measurement 



also requires slopes for subhalos steeper than are indicated in numerical simulations (Erickcek and 
Lawl 12011). 
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11. Summary 



Particle dark matter theories have in many ways become a victim of their own success. A myriad 
of models have been developed that explain a variety of observational phenomena, and because of 
this it is easy to see why researchers not directly tied to this field may have difficulty understanding 
the extent of the progress that has been made. The reality is of course that a significant amount of 
work remains to be done to solve the dark matter problem. As has been discussed in this article, the 
phenomenon of dark matter provides us with the clearest evidence for physics beyond the standard 
model, touching several different areas of modern scientific research, and many lingering questions 
remain as we try to bring together the various sets of observations from astronomy and particle 
physics. 

The process of obtaining a solution to the dark matter problem will likely continue to involve 
many more twists and turns. To appreciate the timescale that is required to fully understand 
the dark matter problem in an ideal scenario, we may compare to related problems in particle 
astrophysics. For example there was a twenty year interval between the time Pauli proposed the 
existence of the neutrino and its confirmed detection by Reines in 1954. Then it was not until the 
turn of the century that experiments were able to measure theoretical properties of the neutrinos 
such as their masses and mixing angles. Further, there was also a thirty year time window in 
between the initial discovery of the CMB to the time when the primary anisotropy was detected. 
Given its similarly broad implications for astrophysics and particle physics, the detection of dark 
matter will certainly generate an equivalent amount of interest, stoking an extended period of both 
experimental and theoretical research. 

It is clear that the phenomena of dark matter is one of the most profound mysteries in modern 
science. Solving this mystery will without a doubt led us down the path to new physics, and be 
one of the most significant scientific accomplishments in generations. One of the more exciting 
prospects in modern science is that, one way or another, in the very near future our prevailing view 
of dark matter will be affirmed, or there will be a paradigm shift in our view of it. 
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